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PREI’ACE 

Pi^PXBiTiJfiS of dieleotrio materials are of interest to scientists 
in various branches: physiciste, chemists, electrical engineers, 
biologists. Their interests concern different aspects; thus an 
electrical engineer Trill require the dependence of dielectric loss 
on frequency and temperature in order to find a substance 
which is nearly loss fcee in a certain range. The chemist on the 
other hand can use this knowledge to draw conclusions on the 
properties of molecules. For these, and for many other pur¬ 
poses, it is imperative to have a theory of dielectrics. 

The present book is intended to give a systematic account of 
the theory of the dieleotrio constant and of dielectric loss. This 
it is hoped will satisfy the requirements of the various branches 
of research interested in dielectrics. In writing thia account I 
found that the subject deserves interest afeo from a methodical 
point of view as an application of classical statistical mechanics. 
That this application is far from trivial is shown by some 
of the controversies in the literatme which have lasted until 
very recent years; also the general theorems derived in § 7 seem 
to be novel. 

It was my intention to write this monograph for the use of 
applied scientists. I hope, however, that the sections dealing 
with the general theory wall also be of value to students. The 
required mathematical technique only occasionally exceeds 
acquaintance with calculus; even so I have been told that its 
extensive use might be too heavy for biologists. The reader is 
assumed to have a certain elementary knowledge of atomic and 
molecular physios, statistical mechanics, and electrostatics. 
Quantum mechanics wUl not be required; its relation to the 
theory of dielectrics is discussed in van Vleck’s book ( VI). 

Units unless stated otherwise refer to the electrostatic c.g.s. 
system. Vectors are represented by bold type. Unfortunately 
it was not always possible to avoid repetition of symbols. The 
meaning of the symbols oc, is respectively propor¬ 

tional to, approximately equal to, order of magnitude of. As 
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usual, h is Boltzmann’s constant, h is Planck’s constant, 
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CHAPTER I 


MACROSCOPIC THEORY 


1. Static fields 


Tke> electric properties of dielectrio substances are usually 
described in terms of the dielectric constant. For most materials 
this quantity is independent of the strength of the electric field 
over a wide range of the latter, but in the case of alternating 
fields depends on the firequenoy. It also depends on parameters, 
such as the temperature, which define the state of the material. 
In the macroscopic (phenomenological) theory, which will be 
summarized in the present chapter, the dielectric constant is 
supposed to be known empirically. The purpose of the rest of 
the book will be to derive the dielectric constant (and its varia¬ 
tion with temperature, frequency, etc.) from the atomic structure 
of the material. 

Throughout this book we shall be interested in homo¬ 
geneous substances only, and in electric fields which are inde¬ 
pendent of the space coordinates, although they may depend 
on time. 


Consider now a condenser consisting of two parallel plates in 
vacutun whose distance apairt d is small compared with their 
linear dimensions, and suppose that they are charged electri¬ 
cally with charges -\-aA and —a A respectively. A is the surface 
area of a plate, and hence -j-a is the charge per unit area which 
will be denoted as surface density. The charges give rise to an 
electric field which inside the condenser is practically homo¬ 
geneous and directed perpendicular to the surface. Its amount is 


given by 


E = 4770-, (vacuum). 


( 1 . 1 ) 


In this equation the factor 47r is due to the particular way in 
which the unit of electric charge is defined, leading to the electro¬ 
static c.g.s. system of units usually chosen in atomic physics. 
Equation 1.1 should be completed by the definition of the electric 
field as the force (in dynes) acting on the unit of electric charge 
which then leads to a unique definition of this unit. 

49S0.il n 
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Between the two plates of the oohdenser a voltage can be 
measured whose absolute value ^ is given by 

0 = \E\d. (1.2) 

Now suppose the space between the plates to be filled with a 
homogeneous dielectiio material while the charge on the plates 
remains unaltered. This will cause the voltage to drop to a 
smaller value, and the ratio of its former value to the present 
one is denoted by the static dielectric constantSince equation 
1.2 still holds, the electric field strength has also decreased, its 
present value being ^ 


From equations 1.1 and 1.3 it follows that the drop in the field 
strength accompanying the insertion of the dieleotrio might also 
be achieved by reducing the surface-density a of the electric 
charge by the amount 



(1.4) 


Therefore the infiuence of the electric field on the dieleotrio 
is equivalent to charging the two surfaces of the dieleotrio with, 
charges of opposite sign in such a way that the positive condenser 
plate is faced by the negatively charged surface of the dielectric, 
and vice versa. The surface density of the charge is constant, 
and amounts to P. This behaviour of the dielectric is to be 
expected from the atomic point of view according to which any 
substance containing no net charge consists of an equal number 
of positive and negative elementary charges. In a dieleotrio in 
particular these charges cannot move &eely through the medium 
(as in a conductor), but they can be displaced. Clearly negative 
charges will be displaced towards the positive plate and con¬ 
versely. The total charge passing through any unit of area 
within the dielectric, which is parallel to the condenser plates, 
is the same and its amormt is equal to P. P, therefore, is called 
the polarization of the dielectric. 

In macroscopic physics the introduction of the polarization 
P through the displacement of charges is rather fictitious 
because these charges cannot be removed from the dielectric. 
They compensate charges of opposite sign but equal absolute 
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value of the condenser plates. These latter charges are some¬ 
times denoted as ‘bound’ charges, 'whereas the original charges 
at the condenser plates are called ‘true’ charges (cf. !Fig. 1). It 
is cus'toinary -then to introduce a new field-quan'tity which is 
desciibed in terms of the true charge. This is the electric displace¬ 


ment, defined by 


D = 4770'. 


( 1 . 6 ) 


In vacurun, therefore, D = E. 


1.3 and 1.5 



In dielectrics, however, using 

( 1 . 6 ) 


Clearly &om equations 1.4, 1.5, and 1.6 


D = E-\-4:ttP. 


(1.7) 


+ + + + l±)_£t2_to_+ + + + + + 



VC? a '=5' — _ — • 

la 1b 


Fio. 1. Dielectrio material with, dielectric constant e, 2 between condenser 
plates, (a) Macroscopic description. The left-hand aide t>f the fig;ure ^ows 
true charges + > — as sources of D. The right-hand side shows that the true 
charges oeua be considered as composed of bound charges 0, 0 and of £ree 

charges | + | . being sources of 47rP and £1 respectively, and leading 'to 

JD = jBr-|-4wi’. 

(b) Atomic description. Only true charges exist, and the field is described 
"by JEj only; the polarization of the dielectric (indicated by -1-1-..* —)» how¬ 

ever, lea^ to surface charges which compensate some of the charges on the 

condenser plates. 

Thus in maoroscopio physics the electric field in a dielectrio 
must be described by two field quantities. The electric field- 
strength E and the electric displacement D are usually chosen, 
and the polarization P can then be derived with the help of 
equation 1.7. D is defined by the (true) charge according to 
equation 1.5, and E can be derived from D wdth the help of a 
special relation 1.6 which is characteristic for the particular 
dielectric material. 

To link up this description with atomic physios it should be 
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noted that the surface charges i PA give rise to an electric dipole 
moment M of the dielectric given by 

M = PAd = PF, (1.8) 

•where F = Ad is its volume. On the other hand, as sho'vm in § 4, 
this electric moment M can be oalcula'ted from the configuration 
of the positive and negative elementary charges which form the 
constituents of the substance. Thus from 1.8 with the use of 
equations 1.6 and 1.7 it follows that 

6,-1 = 4^MIVE == 4^PIE. (1.9) 

This equation provides "the required link between macroscopic 
and atomic theory. 

2. Time-dependent fields 

Now consider that the charges on •the condenser plates, and 
hence the electric field, depend on time. As in the static case, 
a dielectric which is placed between the plates 'will be polarized 
by •the field. The displacement of charges connected 'with this 
polarization usually shows some inertia. Thus if a constant 
field is suddenly applied the polarization •will not reach its static 
value immediately, but •will approach it gradually (cf. lig. 2). 

As in the static case, two field quantities are required to de¬ 
scribe 'the electric field inside •the dielectric and one usually 
chooses the electric field-strength E and the electric displace¬ 
ment Z>. The latter is still defined by equation 1.5, and the 
correlation between E, D, and P given by equation 1.7 still 
holds. Equation 1.6, however, correlating E and D, is no longer 
valid in the present case, but has to be replaced by a more general 
relation. 

Consider first the important case of a periodic field, e.g. 

E = PqCoscu^, (^.1) 

where E^ is independent of time and w/ 27 r is the frequency in 
cycles per second. If a field of this type has persisted for a 
s'officient length of •time, JD too must be periodic in time. In 
general, however, I) will not necessarily be in phase •with E, but 
■will show a phase-shift i.e. 

E = Z>oOos(<o^—^) = Z>]^cos<u<-t-.^2sin<oiS, 


( 2 . 2 ) 
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TIME-DEPENDENT FIEDDS 

according to elementary trigonometry 

Di = X>o cos Da = I>o sin (2.a) 

most dielectrics Do will be proportional to JSg, but the 
-Z>o/-®o usually depends on frequency. Therefore, two 





Tiine-depend0noe of ths polflcization P of a diolootiio wben a constant 
eleotrio field E is suddenly applied to it. 


^^rf^x^nt dielectric constants, ei(n)) and both firequenoy 

^xxdent, can be introduced by 


Di = £iJ?o and Dg = J7, 


(2.4) 


from 2.3 and 2.4 

tan^ = —. (2.6) 

€i 

j-fc will be shown in § 3 that £3 is proportional to the energy loss 
in. ddeleotrios. 

J^s the frequency approaches zero the present description 
im-usb become identical with that given in § 1. Thus (assuming 
•fcliafr ■fchere exists no dielectric loss in static fields) 

£3(0)) —>• 0, *i(t*') as ^ 0. (2.6) 

"W© shall add the further relation (of. also § 10 ) 

(2.7) 

"vvfrioli is to be understood in such a way that e*, is the value which 
e-j. (<o ) a-pproaches at the highest frequencies contemplated in the 
-px:*es©nt book. They correspond to wave-lengths in the infra-red 
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The above equations can be written in a condensed form by 
introducing a complex dielectric constant 

e = ci+ica (2-8) 

and replacing equation 2.1 by 

E = (2.9) 

considering, however, only the real part of this equation (which 
is identical with 2.1). Then the real part of the equation 

.D .= eE (2.10) 

is identical with equations 2.2 and 2.4. 

The two dielectric constants and if considered as functions 
of the fiequency o> are not entirely independent if the relation¬ 
ship between E and is a linear one.f This linear relation¬ 
ship is usually expressed by the principle of superposition and is 
best explained with the help of a more general time dependence 
of the held than has been considered above. Assume that during 
the time interval between u and u-\-du an electric field of strength 
E{^u) has been applied to the dielectric and that the electric 
field vanishes outside this time interval. A displacement D will 
result which in view of the inertia of the polarization P will 
persist at times t > it-j-dw, but which will gradually vanish. J 
D is thus a function of t — u, i.e. 

D{t — u) — E{v,)ot{t — u) du if i >■ u-\-du, 

where a(<— u) is the decay function describing the gradual 
decrease of D, in particular 

oL{t—u)^0 if «->-oo. (2.11) 

The displacement D contains a part which can follow the field 
practically immediately, and which in view of the meaning 
given to €„ will be assumed to be equal to Coo E{u). 

Thus 

D{t — u) = €a 3 E{u)-]r E{u)a{(i) du if u Ct d u-\-du, 

t Of. B. Gross [G4] and S. "Wliitehead [_WS], where further references are 
given. 

j; Note that according to 1.7 D — 4fnP If E = 0. 
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'vehere <x may be oonsidered to remain at the value a(0) during 
the short interval du. 

Suppose now that at a later time interval between u' and 
u'+du' another field E(u*) is applied. Then by the principle of 
supeiposition it will be assumed that the corresponding dis¬ 
placement D{t — u') is superposed linearly on the former one. 
This principle of superposition if applied to a continuous tipae- 
dependent field E{u) initiated at the time u — 0 requires that 
the displacement D{t) at the time t is given by 

D{t) = e„ E{t)-\- \ E{u)oi{t — u) du. (2.12) 

0 

I 

This equation will now be applied to periodic fields. Intro¬ 
ducing JE from equation 2.1 into 2.12 thus leads to 

t 

cosojf =: J a(i—u)oos cju du 

0 

t 

= JIq J a(a;)cos to{t —*) dx 
0 

if a! = t—u is introduced. It should be noted that in all integra¬ 
tions t is to be considered as a parameter. Again it will be 
assumed that the field has persisted sufficiently long to make D 
a periodic function of time. This means that t is larger than the 
time iff at which C3i(f) practically vanishes. Then in view of 2.11 
the above integration over the coordinate x can be extended to 
infinity without appreciably altering the value of the integral, i.e. 

oo 

D{t )—ffoo JS^Q cos cof = Eq j tx{x)<iqa oi{t — x)dx, 

0 

or applying a simple trigonometric formula, 

00 

— E^COBiot = EoCOBOit \ a.{x)<iOB UiX dx 

0 

cX> 

-l-^oSinwt j (x(x'j&Xkaixdx. (2.13) 
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Comparing this expression with equations 2.2 and 2.4 it follows 


that 


00 


«i(a>)—e* == f ai(x)ooaQiX dx 


CO 


«2(<o) = J cx(a;)sin oxb dx. 


(2.14) 

(2.16) 


Both functions ei(cu)—and € 2 (a>) can thus he derived &om the 
same function a(a;) and, therefore, cannot be independent. 

A calculation carded out in the appendix (A, 1 .iii) shows that 

00 

6i(a))—eoo —- f Ca(/x>)- 2 — (2-16) 

. TT J jjr — 

0 


and €a(a»)=- f {gi(/*)—€«} — ^ dfjL, (2.17) 

IT J Ol*- 

0 

where /x is a variable of integration. Both integrals are principal 
values. 

Bquation 2.16 can be used to calculate the static dielectric 
constant &om e 2 (<o). In this case clearly 

00 

«« = «i(0) = Ceo+- f ga(A^) —, (2.18) 

TT J fl 

0 

which indicates that substances for which —€„, is very small 
cannot show appreciable dielectric losses (which are proportional 
to € 2 )- 

I^inally, it should be mentioned that the correlation between 
maecoscopio and atomic theory is provided by an equation 
similar in nature to equation 1.9. It is based on the fact that 
equation 1.7 holds iu the time-dependent case as well as in the 
static one. With a similar argument to that in § 1 it follows that 
the polarization P is equal to the electric moment Jf per unit 
volume, leading to equation 1.8. Introducing the complex 
dielectric constant e with the help of 2.10, using 1.7, it follows 

(€-1)^ = 4^MIV, (2.19) 
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in which only the real part of the expression on the left-hand side 
is considered. Thus if 

M = (2.20) 

then —1 = 4:7T3fjJVEQ, = 4!TT3i^VjEQ. (2.21) 

m 

3. Energy and Entropy 

A. Static jidds'f 

In many text-books the expression 

is stated to represent the electric energy per unit volume of a 
dielectric material with a static dielectric constant «« in the 
presence of an electric field E. This statement is misleading 
whenever depends on temperature. In fact it suggests that 
the energy difference per unit volume of the dielectric, first in 
the presence and then in the absence of an electric field, is always 
given by the above expression. This energy difference ought, 
however, to depend on the state in which the dielectric is kept 
while the electric field is applied; this might, for instance, be 
done isothermally or a<diabatically. A more accurate discussion 
given below shows that the above expression actually is the 
change of the free energy of the dielectric. 

Before this discussion is commenced the reader will briefly 
be reminded of the two fundamental laws of thermodynamics. 
Consider as a simple example a gas of volume v, pressure p, 
and temperature T. By expanding it by the small volume dv, 
work amounting to pdv will be done. Conservation of energy, 
therefore, requires that 

dU — dQ—pdv (3-1) 

is the change in the energy content of the gas if dQ is the influx 
of heat during the expansion. Equation 3.1 represents the first 
law of thermodynamics in this simple case. An analysis shows 
that the quantity dQ is not a total differential, i.e. that no 
iinique function Q of the variables exists such that dQ is the 


t Cf. Abraliam-Beoker [Al, Chapter XI]. 
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differenoe between two neighbouring values and ^ 2 * ^ 

reversible process, the expression 

(3.2) 

does, however, represent a total differential of a function 8^ 
the entropy. ^ is of fundamental importance in connexion with 
the second law of thermodynamics (cf. text-books). 

With the help of 8 the Helmholtz free energy F can be derived 

F^U^TS. (3.3) 

It represents the maximum amount of work which the system 
can be made to do in an isothermal (i.e. constant temperature) 
process. 

. In the case of dielectric material in an electric field, electro¬ 
magnetic theory (of. Appendix A l.i) shows that the quantity 

~EdD (3.4) 

47r 

represents the influx of energy into the dielectric (per unit 
volume) if the displacement D is increased by the small amount 
dD. 

Assume now that the volume of the dielectric is always kept 
constant, and that the temperature T is the only parameter 
besides the electric field F considered to be varied. Then the 
iuCTease dU of the energy U per unit volume of the dielectric 
in a process in which either T or E, or both, are varied slightly 
is given by 

dU = dQ+^dD (3.6) 

477 


if dQ is the influx of heat per unit volume. 

This equation is similar in structure to equation 3.1 for a gas 
if JEf and E are replaced by —p and v. The relation existing in 
gases between p, v, and T (the equation of state) is, however, 
difiEerent from the relation between E, D, and T. For the latter, 
equation 1.6 will be supposed to hold with a dielectric constant 
which may depend on T but is independent of E. Thus 

dD = d{€,E) = e^dE-^Ede^ = e^dE+E^dT, 
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which means that a variation of D may be due to a change in 
the field-strength E at constant temperature, and to a change 
in temperature at constant E. For the following it will be useful 
to take T and E^ as the independent variables. Squation 3.6 
representing the first law of thermod 3 na.amios then becomes 

dQ-\-^d{E^)-\-—^dT = d.U = diE*) -4-^ dT. 

8{E^) ^ ’^dT 

(3.6) 


A further relation •will now be obtained from the entropy law 
according to which dS, given by equation 3.2, is a total differen¬ 
tial. This means that a unique function 8{T,E^) must exist 
such that 

d(J5*), (3.7) 


38 jff, 08 
dS — — dT+—^ 


Thus if it is found that 


d8 = A{T, E^) dT-\-BiT, E^) d{E^), (3.8) 


where A and B are both functions of the two variables T and 
JSr®, the condition that d8 is a total differential requires that 


dB _ dA 
dT ~ d{E^) 


(3.9) 


because both sides of the equation 3.9 are equal to 8^8ldTd(E^). 
Now according to 3.2, inserting dQ from 3.6, 


d8 




8U 

8T 


E^8e„ 
47r 8T 



8U 

d{E^) 



(3.10) 


This equation is of the type 3.8, and equation 3.9, therefore, 
becomes 


8 (1 I 8U 8 (I fdU .©’‘aeAl 

8T\T\8{E^) Stt// 8{E^)[t\8T 4^ ST)]' 

Carrying out the differentiations one finds 

_—/g -I- T 

8(E^) ~ 877\ 8Ty 

In'tegrating with respect to E^ yields the energy density 

£/=0i(T) + (..+ 2’U)^, 


(3.11) 
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where Vq{T) is mdependent of but depends on T a,xid. thus 
represents the energy of the dielectric in the absence of ^ field* 
The entropy 8 can now easily be calculated because on oom- 
paring equation 3.10 with 3.7, both dSfdT and d3/d(^^) R^e 
known if 27 is introduced fiom equation 3.11. Thus 

88 _1 aZ7o E^ 8\ 88 _ 1 de^ 

8T TaT "‘" SttST*’ 8{E^) 9m 8T' 

or integrating ^ (3.12) 

where 8q{T) is the entropy in the absence of a field.. JFrom 
equation 3.3 one finally finds for the free energy 

= (3.13) 

where Eq{T) is the ficee energy in the absence of a field.. This 
proves oiur original contention. 

The above expressions for 27, 8, and F are very instmotive. 
Thus equation 3.13 for the free energy shows (in analogy^ to the 
meaning of J' in gases) that the amount of electric energy avail¬ 
able in an isothermal reversible process is E^/Sir. 

from equation 3.11 for the energy it follows tlxat for 
substances such as dilute dipolar gases for which (cq is inde¬ 
pendent of T) _ g^_^QQjigtajit/T, (3.14) 

the energy change due to the field is given by E^fSir. Tlius in 
this case the temperature-dependent part of e, does not make 
any contribution to the energy. The remaining onergy 

(cg— €0)E^/S7t available besides eQE^fSir is thus entirely duo to a 
change in entropy, f 

!FinaJly, equation 3.12 shows that the entropy is increased by 
the field if 8€g[8T is positive, and decreased if this quantity is 
negative. Since the entropy is a measure of the molecular dis¬ 
order, an external field creates order in dipolar liquids and gases 
for which decreases with increasing T. This may be expected 

l^ear tlie absolute zero of temperature this woulci oontradiob blxe tiliird 
law of thermodynamics; a temperature dependeuoe of the type 3.14= can, 
therefore, not be valid near T « 0, 
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because the field will orientate some of the dipoles which in the 
absence of a field are at random. In some dipolar solids, on the 
other hand, increases with T, which means that an external 
field increases the disorder. This too is understandable if one 
assumes that in the absence of a field the dipoles are in a well- 
ordered state as may be expected in solids. The field by turning 
some of the dipoles into different directions can thus only 
decrease the existing oxd.er. 



Eio. 3. Sobematio temperatiire-depeiidenoe of the dieleotrio oonstaat t, and 
of the entropy change S oc BtfIdT due to the polarization by a field. TfS >>■ 0 the 
field creates disorder, if iS ■< 0 it creates order. Near the absolute zero of 
temperature the substance is already perfectly ordered. Henoe dejBT cannot 

be negative near T => 0. 


B. Periodic fields 

The corresponding calculations in the general case of a time- 
dependent field become very involved. It is easy, however, to 
consider the isothermal case in the presence of a periodic field, 
and to calculate on an average over one period the amount of 
electric energy which is transformed into heat. 

On an average over a period the energy U of the dieleotrio 
cannot change because the temperature is kept constant, and 
the field E is periodic. Thus in equation 3.5 dU = 0, and henoe 
integrating over one period 

[ dQ = —^ ( EdD = —^ I E^dt. 

J 47rJ 47r J at 

o 

The heat produced per second and per unit volume, i.e. the 
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rate of loss X ofem^ from tke electric field, is therefore given 

liy a»/» 

r 0 f «Si) 


8ff» J U 

0 

Inserting here B from equation 2.1 and D from 2.2 and 2.4, 
one finds 

r (oilS r J dm0t , fisincoft 
! = _• oosorfL-— 


or after integration, L ^ 


^ € 2 -^ 0 ^ 


8tr 


(3.16) 


TUs can be expressed in terms of the phase-shift ^ (equation 
2.6) by P 2 

i = ll^tanf (3-16) 

oTT 


For this reason is usually described as the loss angle. 

Another derivation of equation 3.16 is given in the appendix 
(A l.ii). The connexion of and “with the optical constants 
of the material is also discussed in the appendix (A. l.iv). 


4 



CHAPTER II 

STATIC DIELECTRIC CONSTANT 

4. Survey 

In the present chapter the intention is to calculate the electric 
dipole moment induced by an external held in a dielectric from 
its atomic and molecular structure. The dielectric constant e, 
can then be obtained with the help of equation 1.9. 

In § 7 a formula of general validity will be derived which 
connects the static dielectric constant c, with structural proper¬ 
ties of the substance. The calculation of explicit values of e, 
and its temperature dependence is usually, however, beset with 
great formal dihiculties and two types of approximations are, 
therefore, introduced. 

Firstly, a simple model is in general chosen to represent a 
material of much greater complexity. Secondly, mathematical 
approximations are often introduced which hold only within a 
certain range of some parameter such as the temperature. 
Consequently, in order to decide whether an approximate 
formula is applicable for a given material, one has to judge (a) 
whether the basic model can be chosen to represent the actual 
material, and (6) whether the mathematical approximation 
holds for the given range of the parameter concerned. 

A dielectric substance can be considered as consisting of 
elementary charges e^, and 

2e, = 0 (4.1) 

if it contains no net charge. The electiio dipole momentf of a 
charge relative to a fixed point is defined by the vector if 
is the radius vector from the fixed point to . The total moment 
of the whole system is the vector sum of all the individual dipoles, 
2 If This quantity is independent of the position of the fixed 

point if the sum of the charges is zero. For using this fact (cf. 4.1) 
the dipole moment relative to a point at a distance b from the 
original one is 

2 = 2 2 = 2 

’f Also referred to as eleotrio moment or dipole moment or moment. 
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It will now be asstimed that in the lowest energy state (ground 
level) of the substance its dipole moment vanishes. Then if is 
the positional vector of the charge in the ground level. 



(4.2) 


Therefore, if is the displacement of the charge from its 
equilibrium position in the ground state, 1< = l<o+*‘<* Then 

M(X) = |; 1< = r< (4.3) 

is the electric dipole moment of the substance for a given set of 


displacements 


i — (Tj,, 1*2,..., r^,...) 


(4.4) 


which has been e^ressed in abbreviated form by X . Many such 
sets X may lead, of course, to the same moment M. 

Very often it is useful to collect some of the elementary charges 
into a group forming an atom, a molecule, a unit cell of a crystal, 
ox some larger unit. Let the jth unit of this type contain the 
s elementary charges 

ajy (x^,x^g,..., x^jj.,..., x^^) (^•®) 


be an abbreviation for the set of all their displacements x^^,..., T^g. 


Then m(a;j) = 2 ej^x^,. (4.6) 

is the electrio moment of this ^'th group of charges, and 

M{X) = y xnixj) (4.7) 

where the sum extends over all the groups. The vector sum of 
their individual moments in(a;^) thus forms the total moment 
M(X). Our task is to hnd the average displacements, and hence 
the average electric moment tmder the influence of an external 
electric field. 

In order to obtain a preliminary idea about the average con¬ 
tributions of certain displacements to the electric moment we 
shall consider two cases, each of which has its characteristic 
type of displacement: 

Case (i) the displaced charge is bound elastically to an equili¬ 
brium position; 
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Case (ii) a charge has several equilibrium positions, each of 
which it occupies with a probability which depends on the 
strength of an external field. 

The interpretation of case (i) is that on displacing the charge 
s, carried by a particle of mass m a distance r, a restoring force 
proportional to —r acts on the particle in a direction opposite 
to the displacement (hence the — sign). Thus, if a constant 
external field f is applied (t = time) 

_ = -a,$r+-f, (4.8) 


where denotes the frequency of oscillation, and — 

is the restoring fprce. Equation 4.8 can be written 

S(r—r) = —cu§(r—r), (4.9) 


where r *= —(4.10) 

i.e. dHtJdt = 0. The change e, therefore, carries out harmonic 
oscillations about the position r which thus represents the time 
average of its displacement, i.e. if G and S are constant 

r 3r“|“ G cos(cu0 S). 


The average electric moment is, therefore. 


>2 


er 




(4,11) 


As an example of case (ii) consider a particle with charge e 
possessing two equilibrium positions A and B, separated by a 
distance b. In the absence of an electric field the particle has 
the same energy in each position. Thus it may be assumed to 
move in a potential field of the type shown in Eig. 4. If in 
equilibrium with its surroundings it will oscillate with an energy 
of order kT about either of the equilibrium positions, say about 
A. Occasionally, however, through a fluctuation, it will acquire 
sufficient energy to jump over the potential wall separating it 
from B. On a time average, therefore, it will stay in A as long 
as in B, i.e. the probability of finding it in either A or B is 
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Th© presexLoe of a field f "will afiect tliis in tTwo "ways, firstly, 
as in oa^ (i), the equilibrium positions will be shifted by an 
^ount r which for simplicity will be assumed to be th© same 
in A ^d B. Secondly, th© potential energies of the par¬ 

ticle in the two equilibrium positions will be altered because 



^rgy of a charged particle with two eqTiilibrium positions. 
The clotted curve holds in the presence of an external held f. 


its interaction energy with th© external field differs by e(bf), i.e. 

^ = e(bf). (4.12) 

Th© particle should, therefore, on the average spend more 

tune ne^ B than near A. Actually, since according to statistical 

mechanics, th© probabiKty of finding a particle with an energy 
V is proportional to 

^ Pb g-K/Arj I (4-13) 

are the probabilities for positions A and B respectively. They 
ave been normalized in such a way as to make 

Pa+Pb = 1 (4.14) 

m agreement with the physical condition that th© particle must 
be in on© of the two positions. Thus from 4.13 and 4.12 

^ ^ (4.16) 
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It follows from the definition of the probabilities jp^ and ^ 
that if the condition of the system over a long time is con¬ 
sidered, the particle will spend a time (use 4.14) 

in position A, and a time PsH = Ir position 

B. It has thus been displaced by the distance b from A to B 
during the fraction ^{Ps — Pa) of f^e time The average 
moment induced by the field is thus 

^^MPb—Pa)- (^- 16 ) 

Hence if ^ is the angle between b and f, the projection of the 
induced moment into the field direction is, using 4.16 and 4.16, 


given by 


\eb cosd 


^ebfoosOlkT _ 


(4.17) 


gc6/ooe0/fca?_|_i* 

In most cases it is permissible to assume 

e6/ < kT, (4.18) 

for putting e — electronic charge,/= 300 volts/cm. = 1 e.s.u., 
b = 10~®cm. distance between neighbouring atoms in a 
molecule, and T = 300° (= room temperature) one fin ds 

ebf 4-8 X 10-10x10-8x1 


kT 


10 - 4 . 


(4.19) 


1-4 X 10-1® X 300 

Developing 4.17 in terms of ebfJkT, the average induced 
moment in the field direction is found to be 

{^d?)*coa^dflkT-\-ef, (4.20) 

where ef is a term similar to those considered in case (i) which 
has been added to account for the elastic displacement. 

Often two charges -f- e and — e are strongly bound, forming an 
electric dipole p. — ed where d is the distance between them. 
The above case (ii) then leads to the same result as that of a 
dipole p having two equilibrium positions with opposite dipole 
direction, but with equal energy in the absence of a field. In 
a field f the energy of interaction between field and dipole is 

”5' -(Irf). (4.21) 

SO that 2p/oo8 6 is the energy difference between the two positions 
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if 0 is the angle between (i and f. This is equivalent to equation 

jjL = ^eb. . (4.22) 


4.12 if 


Actually putting an immobile charge —e half-way betvit^een A 
and B turns case (ii) into the present case. Clearly the in^np^d 
moment must be the same for both cases because the charge 
—e is immobile, and its distance from A and B is leading to a 
dipole moment /x, 4.22. Introducing this into 4.20 yields folrthe 
induced moment in the field direction 


( 4 . 23 ) 

In contrast to case (i) the electric moment now depends on tem¬ 
perature. In view of equation 1.9, a substance consisting of a 
great number of such dipoles will have a temperature-dependent 
dielectric constant in contrast to a substance in wbioh all charges 
are bound elastically. According to equation 3.12 this means 
that in the dipolar case (ii) the entropy of the substance is 
decreased by the field. This is evident because the field causes 
the fiaction of dipoles with components in the field direction 
to be larger than the fraction of dipoles with components in 
the opposite direction, thus leading to a state which is disordered 
to a smaller degree (i.e. having a lower entropy) than the state of 
complete disorder in which pj^ = p^. 

The difference between the action of the field in the two oases 
(i) and (ii) should be well understood because this is essential for 
the whole theory of dielectric constant. In case (i) the field 
exerts a force on an elastically bound charge, thus shifting its 
equihbrium position. In case (ii) this force of the field on the 
charge again leads to contributions of type (i) denoted by cr in^ 
equations 4.20 and 4.23. It would be wrong, however, to assume 
that the field by this force turns a dipole from one equilibrium 
position into another. This is effected in a more indirect way 
because the field slightly alters the probabilities of a jump of a 
dipole from one equilibrium position to another. This will be 
described in greater detail in § 9, where the dynamic properties 
of the present model will be investigated. 

It should also be realized that though every charge is displaced 
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elastically (case (i)) the fraction of dipoles turned by a field of 
reasonable strength (case (ii)) is very small. This fraction is 
given by —ik) which, according tp 4.16, 4.18, and 4.19, is of 

the order 10~* in a field of 300 volts/cm. Even in a field of 
100,000 volts/cm. only about 2 per cent, of all dipoles are 
orientated. 

5. Dipolar interaction 

I 

In a dielectric two essentially different types of interaction 
forces should be distinguished. Forces due to chemical bonds, 
van der Waals attraction, repulsion forces, and others have all 
such short ranges that usually interaction between nearest 
neighbours only need be considered. Compared with these 
forces dipolar interaction forces have a very long range. This 
can be readily shown as follows. 

As indicated previously (cf. § 4) a polarized dielectric can be 
considered as composed of small regions each having a certain 
dipole moment, and the total dipole moment of the body is the 
vector sum of the moments of these regions. Now it is well 
known from macroscopic theory that the energy per unit volume 
of a macroscopic specimen depends on its shape (cf. Appendix 
A 2.iii). This implies that interaction between dipoles must be 
taken into account even at macroscopic distances and illustrates 
the great importance of the dipolar interaction forces. 

Due to the long range of the dipolar forces an accurate cal¬ 
culation of the interaction of a particular dipole with all the 
other dipoles of a specimen would be very complicated. How¬ 
ever, very good approximation can be made by considering that 
the dipoles beyond a certain distance, sayt ciyn> o®-*! replaced 
by a continuous medium, having the macroscopic dielectric 
properties of the specimen. Thus the dipole whose interaction 
with the rest of the specimen we are calculating may be con¬ 
sidered as surrounded by a sphere of radius cont ain ing a 
discrete number of particles, beyond which there is a continuous 
medium. To make this a good approximation the dielectric 
properties of the whole region within the sphere should be equal 

f The suffix m stands for xnacrosoopic. 
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to those of a maorosoopio speoimen, i.e. it should coutain 
a sufficient number of molecules to make fluctuations veiy 
small. 

jProm this point of view we are thus led to H. A. Lorentz’s 
[Zr5] method for the tmatment of dipolar interaction: fiom a 
maorosoopio specimen select a microscopic spherical region 
which is sufficiently large to have the same dielectric properties 
as a macroscopic specimen. The interaction between the dipoles 
inside the spherical region will then be calculated in an exact way, 
but for the calculation of their interaction with the rest of the 
specimen the latter is considered as a continuous medium. To 
demonstrate this method we shall make use of a very simple 
model. It consists of a cubic lattice of * atoms’ whose linear 
dimensions are very small compared with the lattice distance. 
Sack atom consists of a positive charge -\-e which is rigidly 
bound to the lattice point, and of a negative charge —e which is 
bound elastically to it. The force acting on a negative charge 
when it is displaced while aU others are at rest at their respective 
equilibrium positions will be denoted as the restoring force. When 
other charges are displaced as well, an additional force will act 
on the charge due to a change in the interaction between the 
charges. It will be assumed that electrostatic interaction only 
exists, and in particular that there are no short-range forces. 
iFurthermore, the temperature will be assumed to be so close to 
the absolute zero, T = 0, that thermal oscillations can be dis¬ 
regarded. 

On these assumptions when a maorosoopio electric field E is 
set up in the substance all negative charges are displaced by the 
same amount, say r, each forming a dipole 

m = (—c)r. (5*1) 

The field f acting on one charge and displacing it against the 
restoring force is often referred to as the local field, or the inner 
field; it has to be distinguished from the macroscopic field £. 
Thus if the restoring force is denoted by —c®r. 
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and hence according to 6.1 

m=^f. (6.2) 

To calculate f according to Lorentz’s method, separate the 
sources of f into those inside the spherical region and those out¬ 
side it, making contributions and respectively, 

f = (6.3) 

Since aU induced dipoles are identical, the field f must be the 
same at every lattice point. To find consider therefore the 
dipole at the centre of the spherical region and calculate its 
interaction energy I with all the other dipoles of the region. 
This will depend on f , and 

=—grad7(f). (5.4) 

To calculate this interaction energy I assume the field to be 
sufficiently weak to make r very small compared with the lattice 
distance a^. Sach dipole xn can then be considered as a point 
dipole. The electrostatic interaction energy between two parallel 
dipoles of equal moment is given by 

^(l-3cosV), (6.6) 


where I is the distance between them and tp is the angle between 
1 and m. For simplicity let m be parallel to a crystal axis, say 
in the z-direction. Then the three components of 1 are 

X = nag, Y = pUq, Z = qa^ 

where n,p, q are positive or negative integers. Since cos ^ = Zjl, 
the total energy of interaction follows from 6.6 by summation 
over all lattice points, i.e. 


m* 


2 






0 n.3>,<r 


(TO®-|-p®-l-g®)** 


Now to each set of three values n, p, g, say n n^, p = p^, 
g = go, two others can be coordinated by cyclic permutation, 
namely n = p^, p = q^, g = and % = go, p = g == po- 
These three terms just cancel in the above sum, and hence 7=0. 
It thus follows, using 6.4, that 

/i = 0. 


( 6 . 6 ) 
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The external contribution has to be calculated maorosoopi- 
oally. It represents the electric held inside the spherical region 
due to aU sources except the polarization inside this region. To 
understand this clearly it should be remembered that the 
macroscopic electric held £ is partly due to true charges outside 
the specimen (or at its surface) and partly due to the polariza¬ 
tion P of the dielectric which acts in the opposite direction. 
To obtain the contribution to £ due to the spherical region 
should be omitted. Therefore if is this contribution, 

f, = E-E.. (6.7) 

By dehnition, Eg, the self-held, is the held inside a spherical 
specimen with a permanent polarization P. Its direction is 
opposite to P. A simple electrostatic calculation (Appendix 
A 2.21) shows that 

E.- ^P. (6.8) 

Therefore, using 6.3, 6.6, 5.7, 6.8, and 1.9, 

f = fg = E+^P = f^E. (6.9) 


It should be noted that this expression is independent of the size 
of the spherical region, which means that the contribution of any 
homogeneously polarized spherical shell to the field inside it 
vanishes. In other words, the interaction between such a shell 
and a dipole inside it vanishes. This is in agreement with the 
above result 6.6 obtained by considering dipole-dipole inter¬ 
action in detail. The present model, therefore, has been chosen 
in such a way as to make macroscopic and mioroscopio treat¬ 
ment identical. Such a result, of course, could only be obtained 
by assuming short-range forces to be entirely absent. The exis¬ 
tence of such a force would modify but it would not alter our 
result for f 


The fact that the interaction energy of the dipoles of a spherical 
region vanishes must mean that if a spherical specimen is brought 
into a homogeneous field Eq in vacuum this field represents the 
local field acting on each dipole. Actually from electrostatics 
it follows that the field inside a spherical specimen of dielectric 
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constant is B — 3Eo/(c,+2) (of. Appendix A 2.16), and hence 
with equation 6.9, f = E^. 

The above restdt 6.9 can be derived by an alternative method 
due to Onsager [01] which is capable of a generalization and 
which will be of use at a later stage of the development of a 
general theory. 

The reader should realize, however, that in Onsager’s paper 
the spherical region contains one single dipolar molecule only. 
This is, of course, open to the objection that the spherical region 
is certainly not large enough to have the same properties as a 
macroscopic specimen. In our use of his analysis, however, the 
spherical region is sufficiently large to make it certain that it is 
correct to treat the outside macrosoopioally. Now consider, 
in the absence of a macroscopic field, a spherical region which is 
polarized homogeneously, thus having a dipole moment M = PF^ 
where F is its volume. If considered separately without the rest 
of the specimen the electric field inside the spherical region is the 
self-field Eg, equation 6.8. If, on the other hand, the spherical 
region is considered inside the specimen there will be a certain 
interaction between the polarized sphere and the surro undin gs, 
giving rise to an altered field inside the spherical region. Its 
difference from Eg is denoted as the reaction-field R and is given 
by (Appendix A 2.18; = radius of the spherical region) 

P 2(gg-l) M 2 eg-1 47rM ^ 2 eg-1 

2eg+l al, 3 2eg+l F 3 2fg+l 

( 6 . 10 ) 


Thus R is the field produced inside the spherical region by the 
surroundings if polarized by the former. If a macroscopic field 
E is now produced in the specimen without altering the moment 
M of the sphere, the field ins ide it will be increased by an addi¬ 
tional field, the cavity-field G (cf. Appendix A 2.16), 



( 6 . 11 ) 


Thus the total field inside the sphere, due to outside soturces, is 
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This expression holds whatever the value of P. In our particular 
case if 47rP = (c,— 1 )£ (of. § 1) is introduced, the above expression 
becomes identical with 5.9. 

Pinally to calculate insert f firom 6.9 into 6.2. The polariza¬ 

tion P is then given by 

P = JVim = 5 s± 2 f isr.E, 


where is the number of particles per unit volume. 

e,-l-2 3 


Hence 

(6.13) 


which is usually known as Clausius—Mossotti formula (of. 
Clausius, • G1 ; Mossotti, MS). The above derivation of this 
formula is exact except for the assumption that f is a small 
quantity. It will be noted, however, that the change in e, as f 
increases is essentially a manifestation of the held dependence 
of the dielectric constant (since the change in magnitude of f 
depends only on held strength) and this does not concern us here. 
That is, formula 5.13 is correct for the limiting case in which the 
dielectric constant is independent of the held. The reader should 
realize, however, that this derivation holds for the above model 
only; there are few substances which it could claim to represent, 
though for simple non-polar substances it will often be a useful 
approximation. 


6. Dipolar molecules in gases and dilute solutions 

Molecules can be divided into two classes, polar and non-polar, 
according to whether or not they possess an electric dipole 
moment in their lowest energy-level (ground state). In general 
it is easy to recognize the class to which a molecule belongs 
because a non-polar molecule must have a point of symmetry 
dehned in such a way that the distribution of charges along (or 
near) any straight line passing through it must be symmetrical 
with respect to this point. Thus a diatomic molecule is polar 
unless its two atoms are equal (e.g. the polar molecules HCl, 
CO, and the non-polar molecules H 2 , O 2 ). Triatomic molecules 
of the type AB^, where A and B represent different atoms, are 
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polar unless their nuclei lie on a straight Hne with A half-way 
between the two B atoms. Examples are the triangular polar 
HgO and the straight non-polar COg molecules. CoH, is a more 
complicated non-polar molecule forming a plane hexagon with 
the centre as point of symmetry. On replacing one H atom by 
another type of atom, say by Cl, the resultant molecule CeHsCl 
becomes polar. 

The magnitude of molecular dipoles is usually of the order 
of one electronic charge (<^4'8xl0-^® e.s.u.) displaced by 
Jx 10“® cm., i.e. about 10““ o.g.s. imits. Often dipole moments 
are expressed in Debye units, one such unit being 10““ c.g.s. 
units. Debye was the first to recognize the importance of the 
investigation of dipole moments for a study of the constitution 
of molecules (cf. his book, reference D2). A detailed discussion 
of molecular structure is, however, not intended in the present 
book. Our purpose in this section is to determine the properties 
of dielectric substances consisting of molecules having per¬ 
manent dipole moments and our model for such a molecule (in 
vacuum and ftee from perturbing influences) will, therefore, 
have a dipole moment (i.«. 

In addition to its translational motion, a free molecule can 
carry out oscillations and rotations. However, unless stated 
otherwise, it will be assumed that these do not alter the average 
value of the dipole moment. In a complicated molecule intra¬ 
molecular rotation of dipolar groups may lead to considerable 
changes in the dipole moment of the molecule as a whole (of. 
the book by Le F^vre, L2). Such molecules wiU not be considered 
at present. 

A constant electric field f will influence the molecule in two 
ways. Firstly, it will perturb the free rotation- of the dipole, 
and secondly, it will induce a further dipole moment, say of, 
by elastic displacement of the atomic electrons relative to their 
respective nuclei, and to a smaller extent by elastic displacement 
of the nuclei relative to each other. The total moment of the 
molecule is thus 

m = lA^+of. (6.1) 

The quantity a. has the dimensions of a volume and is called the 
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polarizability of the molecule. In anisotropic molecules (that is, 
molecules having different polarizabilities along different axes) 
the induced moment need not always have the same direction 
as the field f; the polarizability a is in this case not a scalar but 
a tensor quantity. Its average value, a, obtained by edlowing 
the molecule to have all possible directions relative to the field, 
is, however, a scalar. Assuming that the polarizability of a 
■ molecule is entirely due to electronic contributions, a is deter¬ 
mined by the optical refractive index n. For, according to Max¬ 
well’s relation, n* is the dielectric constant at optical frequencies 
at which there arC no dipolar contributions (because the time 
required by the dipoles to attain equilibrium with the field is 
much longer than the period of the field). Assuming the mole¬ 
cules to be isotropic, it will bo shown below that the Clausius— 
Mossotti formula holds approximately in this case, i.e. 


a — a 


a 


1 


where a is the radius of a sphere which on an average contains 
one molecule. 


Gases 

We shall now proceed to calculate the dielectric constant e* 
of a gas of dipolar molecules. In order to simplify calculations, 
it will be assumed that the density of the gas is so low that the 
dipolar interaction energy is small compared -with the thermal 
energy (£^ kT per molecule) and can therefore be neglected. 

According to § 6 the dipolar interaction energy is of the order 
/LiJ/P if is number of molecules per unit volume. 

It will thus be assumed that 

(6.3) 

This simplifies the calculation of considerably because it 
implies that the local field f, acting oh a dipole, is entirely due to 
external sources, i.e. f = D == e,E (cf. § 1). Also if AJ, is small 
c,— 1 must be small as well, so that the further assumption 

(6.4) 

leads to f = E. (6.6) 

It will be found below that 6.4 follows from 6.3. 
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In order to calculate the electric moment M of the gas we make 
use of the fact (§ 4) that M is equal to the vector sum of the 
moments m of all- N molecules. Since, on the other hand, the 
time average m of m is the same for all molecules. 


M = iVm. 

(6.6) 

Thus from 1.9 e,— 1 = 4TTNf^mjE, 

1 

(6.7) 

1 

where according to 6.1 and 6.5 

- 

m = {i^+oiE. 

(6.8) 


In this expression the first term represents the average value of 
the intrinsic moment, and the second one is the average induced 
moment per molecule. The calculation of m is similar to that 
carried out in case (ii) of § 4. In the present case, however, 
instead of having only two possible directions as described in 
§ 4, a dipole has a whole continuum of possible directions. As 
in § 4, the behaviour of the dipole may be considered on a 
statistical basis, without inquiring into its dynamic properties 
under the influence of the field £!. Thus since —cos 6 is 
the energy of the dipole in the field, the probability of finding 
jjL^ in a direction forming an angle between d and d-{-d0 with E 
is given, according to statistical mechanics, by {2Train 6 d6 is the 
element of the solid angle between 6 and d-\-dd) 


OikT gJn d dd j j gBlfi^fioaOlkT gin 6 d6, (®-9) 

provided that thermal equilibrium has been attained. At 
present in considering static properties this will always be 
assumed to be the case. 

In evaluating 6.9 we shall assume, as in the case of equation 
4.18, that the field is sufficiently weak to make 


kT 


< 1 - 


( 6 . 10 ) 


Then developing 6.9 in terms of EjkT, and keeping terms up 
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to the first order only, one obtains for the average value off 
cos^, 

~ ~ E 


•JT 


OOS0 


J 


008 sin 6 dd 


\I j ^Ett,ooaeikT sin d • 


( 0 . 11 ) 

Similarly the components of |i.^ per^ndicular to E are found 
to vanish, so that the average moment of a molecul e is a vector 
with the same direction as E, amounting to p.^ = cos 0, i.o. 

( 6 . 12 ) 


2 


a = Jfi-E. 
3JfcT 

On using 6.7 and 6.8, therefore, 


*a 


1 = i^^« + 4,rSy«, if ..-I < 1 


(6.13) 


or introducing €«> = l+4«-a^ (6*1^) 

as dielectric constant at frequencies so high that the dipoles have 
no time to attain equilibrium. 


_ _ 4^i4No if€ _1 1 

€oo ^JcT * ^ ^ 


(6.16) 


In agreement with the conclusions drawn in § 4 the dipolar 
contribution to e, is temperature-dependent in contrast to the 
non-polar contribution Measurement of the temperature- 
dependence of thus makes possible the separation of the dipolar 
contribution, — e^o, and hence the determination of (examples 
in § 16). 


Dilute solutions 

The above derivation suggests that a formula similar to 6.16 
should hold for dilute solutions of dipolar molecules in a non¬ 
polar liquid having the temperature-independent dielectric 
constant cg. For at sufdciently low densities the interaction 
between dipoles can again be neglected. Also, replacing 6.4 by 

e*—(6.16) 


f For with x = cos 6 and y = JSjhT 1, 



ev^dxcn 


1 

J (as+yas*) 


—1 



dx 


0“f- 

2 


y/3. 
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ir. f 6 


equation. 6.5 holds as before, leading to 6.11 because in a liquid, 
in the absence of an external field, all directions of the dipole 
are equally probable. Furthermore, if too is again the non-polar 
contribution to then —l)^?/47r represents the non-dipolar 
contribution to the electric moment of the solution per unit 
volume. Since —l)j©/47r is the total moment, (e^—eoo)-®/4w» 
in view of the additivity of moments, is,* as before, the dipolar 
contribution. Assuming that the dipole is rigid, i.e. that a = 0, 
equation 6.16 would follow with e^o — eo* Actually a molecular 
dipole is not rigid, but can be polarized. In a solution this leads 
to a change in the effective dipole moment of the molecule, as 
was recognized by various authors (cf. Weigle, WJ; Frank, Fi; 
Higasi, 112; Frank and Sutton, I'd; these authors also consider 
other effects which may lead to an alteration in the effective 
dipole moment when the molecule is non-spherioal or has a 
large quadripole moment). The dipole polarizes its surroundings 
which in turn produce a reaction field at the position of the 
dipole. This field polarizes the molecule and thus alters its dipole 
moment. Thus if we define the resultant moment which the 
molecule has in the solution as the ‘internal moment* p.^, the re¬ 
action field will be proportional to p.^, say g'p.^. In the absence 
of an external field, therefore, from 6.1 (using f == and 


m = p^) 


Pi — P'»~l~0^P’i> 


(6.17) 


or Pi — — . (6.18) 

i—ag 

A quantitative calculation of the reaction field, and hence 
of p^, is very difficult, however, except in the case of spherical 
molecules. 

It follows from the above that, when a field E, produoed*by 
external sources, exists in the solution the resultant field acting 
upon a molecule in the solution will be different from E. Also 
the moment pj of the molecule in solution will be different from 
that it has in vacuum (namely p„). In order to take account of 
these facts in calculating the dielectric constant of the solution 
we may consider the dipole in one of two ways. In the first of 
these we calculate the field acting upon the ‘internal moment* 
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itf. This will, of course, depend on the shape chosen for our model 
of the molecule. In the second method we calculate the inter¬ 
action between the held of the dipolar molecule in the solution 
around it and the applied field E. The same dipolar field might 
alternatively be produced by a rigid dipole having in vacuum a 
moment |j,g which will be denoted as the ‘external moment* of 
the molecxile.f A solution of rigid dipoles fJLg should, therefore, 
have the same dielectric constant as the solution which wo 
actually consider. 

The former method will be used below in deriving the Onsager 
formula. At present the latter metho^ will be chosen because it 
enables us to make use of our previous derivation for gases to 
determine the dieleotiio constant of the solution, i.e. to employ 
equation 6.16 after replacing jji.,, by [t^. Thus with the use of 6.18 

___ Y 

* “ SkT ~~BkT [T^l * 

^ ««—«<» <1 and e<o —€0 •< 1. (6.19) 

Here the ratio of external to internal moment depends on 
the shape of the molecule. For spherical molecules, according 
to the defimtions of fig and /u,j, the internal moment is identical 

with the moment of a sphere inside the dielectric containing 
at its centre. An ordinary calculation in electrostatics shows 
that in this case (of. Appendix A 2.31) 






2€,-hl 




( 6 . 20 ) 


and if a is the molecular radius (of. Appendix A 2.19) 


2e,+ l a» 


( 0 . 21 ) 

to disnl^mp^^^r polarizability is isotropic and is mainly due 
refi»ive index n pS.id'rfiTi'* 

With 6.18, using 6.21 6 2 the molecule. Then 

molecule is ^ Jiitemal moment of a spherical 

lit =B { W’^4-2\ 

V2e,+na~r~;**^’ (6.22) 
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e^(w.«+2) 


(6.23) 


The dielectric constant in this case is (of. 6.19) given by (replacing 
eg by cq on the right-hand side according to 6.16) 

477^4-2'^2)'® 


'OO 


/go(^*+2)Y 

3ArT \ 2eo+n® ) 

H^)‘[ 


(6.24) 


_ 4,7/^i^e /fo±2\*/i 2 (co-1)(6o-w®) \® 

(2€0H-n®)(e0-i-2)/ 


SkT 


for spherical molecules and €g —!• 


In some cases we might expect the approximation of spherical 
molecules to give a correct order of magnitude for the deviation 
from unity of the ratio of the dipole moments in vacuum and in 
solution, but we should not expect to obtain a more accurate 
result than that. 


Onsager formula 

In the case of spherical molecules Onsager [02] has shown 
that it is possible to go one step farther in the approximate 
calculation of the dielectric constant. The interaction between 
molecules wiU no longer be entirely neglected, but one com¬ 
ponent of it, namely the long-range dipolar interaction, will be 
taken into account. The following assumptions will, therefore, 
be made: 

(а) A molecule occupies a sphere of radius a and its polariza¬ 
bility is isotropic; 

(б) The short range interaction energy is negligible (i.e. 

kT per molecule). 

Assumption (6) means that the surroundings of a molecule 
will be treated as a macroscopic continuum with dielectric 
constant e^ because long-range forces only will be considered. 
Estimates on the range of validity of this method can be made 
(cf. reference F12), but will be postponed until the more general 
Kirkwood formula is derived (§8). 

On the assumptions made above, the contribution of a single 
molecule to the dielectric constant €g can be calculated in the 
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same way for a piue liquid or for a mixture. In both cases the 
looctl field f acting on the molecule is the field inside a cavity of 
radius a within a continuous medium of dielectric constant e^, 
and is composed of (i) the cavity field G due to external sources, 
and (ii) of the reaction field R due to the moment m of the mole^ 
oule itself. Both G and R were discussed in § 5, although there 
the cavity of radius d contained a great number of molecules. 
To use the same values for G and R in the present case is only 
possible because of the assumptions (a) and (6). Thus, according 
to 6.10, 6.11, and 6.12 (M/a^i of § 6 corresponds now to m/a®) 
and using 6.21 

f = G+R = (6.26) 

26 ,+ 1 

Inserting this value into 6.1 the moment m becomes 


m 




36, 


8 




osE+o^in, 


or solving vsdth respect to m, 


m 




cxE 


1—otg 26 ,+ ! 1—aff’ 

With this value for m the internal field f (6.25) becomes 


(6.26) 


f = 


36 


8 


£ 


-f 


ff 




(6.27) 


26,+1 1 —ag * 1— ag^' 

To calculate the average polarization m we shall require 
Considering that 



36, cos 0 
26,+1 1 —ag 


Q 

1 — acflr 



(6.28) 


is the energy of a dipole in the field f, the probability of finding 
it in a direction forming an angle d with £ is (similar to 6.9) given 
by 

^ / w 

eMftT smede / J sin 6 dd, (6.29) 


where should be introduced from 6.28. It wiU be seen that 
the second term in 6.28 is independent of 6 and therefore in 6.29 
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only the first term remains. Following the same procedure as 
in the derivation of 6.11 one finds 


and hence 




3g8 _ fJ'yE 

2e.+ l ZkT(l-agy 

3^8 _ 

26^+1 3kT{l—ag) ‘ 


(6.30) 


Thus from 6.26, using 6.30 (and 5c = oc for isotropic polarization), 
the average moment m of the molecule is given by 


m = 


3e 


a 


( 




2e,+ l \3;fcr(l—agr)a 


+ 


l~ocg) 


£. 


(6.31) 


Let us consider first the application of 6.31 to the simple case of a 
pure liquid of non-polar molecules for which — 0. Since by 
definition a molecule occupies an average volume of 


47r 

"s' 




(6.32) 


equations 6.7 and 6.31 yield 


a 3 
2 €g-f-l 1 —ag a* 


(6.33) 


Introducing g from 6.21 and solving with respect to <x leads to 
the Clausius—Mossotti equation 


-1 _ Q£ 

g«+2 “■ 


(6.34) 


This equation can be considered to prove 6.2, which, together 
with 6.21, transforms 6.31 into 


m = 


Sgg E f 2gg-|-1 



(6.35) 


The Onsager formula for a pure dipolar liquid is then obtained 
by inserting 6.35 into 6.7, making use of 6.32: 



—- 3ea 47rjMgi\^/n^4-2 \^ 

26,+71» SkT \ 3 J ’ 


(6.36) 


For very small densities this expression becomes identical 
with 6.15, as should be expected (because for Cg —1 1 and 
w*—1 < 1, (ri2-j_2)/3 1 and 3e,/(2e„H-l) 1). 
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Consider next a mixtuire containing N^, . Ng,..., mole- 

oules per o.o. of z different oompoxinds. Then 2 is 

8 — 1 

total moment per o.o. where iHg is the average moment of a 
molecule of the type s obtained &om 6.36 by inserting for (jt^, 
n, and a-the respective values and a, which these quantities 

have for the sth kind of molecule. From 1.9 the dielectric con¬ 
stant of the mixture thus follows: 

eg —1 = 47rNg rhgfE. (6.37) 

As a particular case take a mixture containing per c.c. 
polar (/Api, n^, ciy) and non-polar (/i „2 — ^a> ^a) molecules. 

Inserting and from 6.36, equation 6.37 leads to 

3^/49t/4j^i / 2€,+1 »ii+2\s 


’8 


1 = 


\ SJcT 
2e.+ l 




3 


■)' 


+ 




}• 


(6.38) 

For a dilute solution of polar molecules, i.e. when N^, 

this expression should become identical with 6.24 if n = n^, 
= ^o, and N-^ — Nq. Actually since 47ra|.^/3 represents the 
volume occupied by the non-polar molecules in 1 c.c., 

means that 4e7Ta\N^Z 1. Since the second term in 6.38 can be 
neglected, we find with the above substitutions for Wg, and N^, 


'8 


1 = 


3e, 


/- _|V I 3 ea 477^1^^ / 263+1 

' ® ^“•■2€^-1-1 BkT \2ca-fn* 3 / ’ 


2e,-t- Cq 

which, after a simple transformation, becomes identical with 
6.24 (since e, £=i Cq cri e«,). 


7. Greneral theorems [Fid] 

The foirmulae derived in the previous section for the dielectric 
constant e, have the great advantage of representing e. in a very 
simple way with the help of a few parameters. It should not be 
forgotten, however, that they only hold subject to certain con¬ 
ditions which in many practically important oases are not ful¬ 
filled. In the present section, therefore, expressions for the 
static dielectric constant e, will be derived which hold in a very 
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general way for any dieleotiio substance which is not per¬ 
manently polarized. Clearly the derivation of such expressions 
requires some mathematical abstractions, but in view of their 
importance it was thought desirable to give a detailed aocotint 
of their deduction. 

As in § 6, select from an infinite homogeneous specimen a 
macroscopic spherical region of volume F which is large com¬ 
pared with a region just big enough to have the same dielectric 
properties as a macroscopic specimen. The surface of the 
spherical region need not be an exact geometrical sphere; 
deviations of a molecular magnitude are permitted and the 
exact surface will be laid in such a way that no moleotiles are out 
by the surface. This has no influence on the field at distances 
which are large compared with atomic dimensions. We shall 
calculate the projection of the average electric moment of 
the sphere in the direction of the macroscopic field £. For this 
purpose all particles inside the sphere will be treated according 
to the rules of classical statistical mechanics. The outside, 
however, will be considered as a continuous dielectric described 
by the macroscopic dielectric constant e^. It will be assumed 
throughout that the macroscopic field E is sufficiently weak to 
prevent saturation, so that e, is independent of E. 

The spherical region consists of a number of elementary 
charges each of which can be described in terms of its displace¬ 
ment from the position it would have in the lowest energy-level 
(ground state) of the whole system. This displacement is a vector 
quantity and is denoted by r^. A set of all the displacement 
vectors will be collectively denoted by X, according to equation 
4.4. Except at the absolute zero of temperature, such a system 
of particles does not stay in the same configuration, even if 
macroscopically it is in equilibrium. Owing to thermal fluctua¬ 
tions there is a probability 

e~u{X,micTdXH e-mx.micTdX (7.1) 

of fintling it with any set of displacements in a space element 
between 

X — (r^, Tg,..., and X-\-dX = (ri-f-dri,..., 
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Here, ?7(X, E) is the potential energy of the system in the con¬ 
figuration X, in the presence of a fi.eld E, and 

dX = dTidr^..,,dr^... (7.2) 

is the product of the volume elements of all the displacements 

each being again a product of its components dr^dr^yd/ti^. 
The integration must be carried out over all possible values of all 
the displacements. 

Sach set of displacements X leads to a dipole moment M(X) 
as sho^m in equation 4.3. Therefore, if 0 is the angle between 
hl(X) and £, JIf(X)oos^ represents the projection of M(X) in 
the direction of the macroscopic electric field £, and is 
the average value of JIf(X)oosd. Thus using 7.1, 

= j M(X)oob e dX / J e "uiXMkx ax. (7.3) 

The value U (X, 0) of the energy in the absence of a macro¬ 
scopic field •will be denoted by U{X), and the zero of po'tential 
energy •will be chosen so that U(X) vanishes in the groxmd state, 
i.e. •when all displacements vanish. The energy U{X) can be 
considered as composed of (i) the energy of interaction ^(X) 
be-tween the particles of the spherical region, and of (ii) their 
energy of interaction Ug{X) with the external region, 

U{X) = Ui{X)-\-U,{X). (7.4) 

C^(X) ■will depend not only on X, but also on‘the dielectric 
cox^ant €g because the external region is to be treated on a 
macroscopic basis. 

The fact that UJ,X) contains the parameter which may 
depend on -temperature requires some consideration. Our 
present investigation is a special case of a wider group of pro¬ 
blems—^the sta-fcistioal mechanics of systems containing tempera- 
t-ure-dependent parameters—which has been investigated by 
Gross and Halpem [01]. Applied to our case their investigations 
show that €g has to be treated as constant parameter whenever 
difiEerentiations "wifch respect to temperature are required; 
UJJX) is then the energy required to establish at constant 
temperature the equilibrium polarization in the outside region. 
TVom the point of view of the whole system (spherical region 
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-1-outside region) UJi^X) is thus the oorresponding &ee energy 
composed of the energy of interaction of the charges inside the 
sphere with the polarization induced by them in the outside 
region, and of the free energy of this region (so far as it depends 
on X), as shown in the subsequent example (cf. 7.17 and 7.18). 
This is reasonable because U{X) is the energy required to estab¬ 
lish the displacement X at constant temperature. 

Now for a given configuration X, assume the macroscopic 
field E to be applied (which, of course, alters the probability 
of this configuration). Inside the spherical region this leads to 
an additional homogeneous field G (the cavity field) given by 
equation 5.11. It follows from electrostatics that its interaction 
with the charges of the spherical region is given by 


—M(X)G == 




M{X)E ooaO. 


Thus in the presence of a field 


(7.6) 


U(X, E) = U{X) — ■ M(X)E cos B. (7.6) 

2ea-l-l 


This expression for U{X, E) must be inserted into equation 7.3. 
It should be remembered now that E was to be considered as 
sufficiently weak to prevent saturation. Thus if the right-hand 
side of 7.3 is developed into a power series in Ey only the first 
term need be considered, i.e. 


g- u{x.mkT ^ e- A t -cos ^-h...V 

\ 2e,-(-l kT J 


Now 


(7.7) 

(7.8) 


J M(X)coBe = 0 

because this integral is proportional to the average moment in 
the E-direotion in the absence of a macroscopic field. Therefore, 
inserting from 7.7 into 7.3, making use of 7.8, 




3e. EJ 


2ea-t-l kT 


J M^{X)ooa^e e- u(X)Iht dXy (7.9) 


where 1/J = f (7.10) 

Assume now that the field E may have any direction relative to 
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a given ]VI(JSr), and average over all these directions. Then 
oos*d will be replaced by its average value, i.e. by •J. Inserting 
this value into equation 7.9, and mairing use of 1.9, yields 


4»r 3e, 

3F 2€,-}-1 kT* 


(7.11) 


where M* ^ J { dX (7.12) 

is the average value of M\X) in the absence of a macroscopic 
field. Equation 7.11 represents our first general result. It 
sho ws that the dielectric constant can be expressed, in terms of 
-ftf the mean square of the spontaneous polarization of a sphere 
of the dielectric embedded in a large specimen of the same 
material. 

Before developing the general theory further, equation 7.11 
will be shown to be self-consistent. This means that 7.11 must 
be fulfilled identically if the spherical region is treated on 
a macroscopic basis. For this purpose the following theorem 
will be employed (of. reference Tl)^ Suppose the free energy 

of a system depends on a number of macroscopic 
parameters Then the probability of finding the system 

in a range doL-^dot^... is given by 


ern«u«t.„yikT J 

Also it will b© noted that the quantity {jdoty^da^...) in the above 
expression represents the volume of the ^cx-space’ element 
between (a^, ocj...) and a 2 +da 2 ,...)- the present case, 

since the free energy depends on the absolute magnitude of M 
* this element is the shell between the two spheres of radius 
and ]i£*\-‘dM. in the -3f-space, and has the volume dM* 

Since -If may in theory have any value between zero and infinity, 

the probability that the spherical region will have a moment 
between M and M-\-dM is 


00 


M^dMj j » dM. - 

0 

_ J? 

Hence Jlf * == J j[f j J jf 2 


(7.13) 

(7.14) 
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CompaTUig this with equation 7.12, and TemembeTing that 
according to 3.3 the entropy S is given by S = {U — F)IT, it 


follows that 




(7.16) 


represents the number of states of the spherical region, in a range 
between M and M -{-dM. This equation might have been used as 
an alternative starting-point for the macroscopic considerations 
based on 7.12. 

Following the procedure used in the case of the energy U(X), 
the free energy F{M) can be considered as composed of an in¬ 
ternal fipee energy Ji(Jf), and an external one, FJI^M), due to the 
interaction "with the surroundings. The former is the self-energy 
Fg calculated in the appendix (A 2.37), 


Fi{M) = 


^ttM^ €,+ 2 

3F eg-V 


(7.10) 


the latter can be obtained with the help of the reaction-field R 
(cf. equation 6.10). It is composed of the energy —MR of the 
dipole M in the field R, and of the free energy required to polarize 
the surroundings accordingly. This latter should bo propor¬ 
tional to 22®, say jS22®. Then 

Fg{M) = —MR4-i822®. (7.17) 


To find j8, Fg can be considered to depend on the parameters 
Bg., By, Bg [for a fixed value of Jlf], and hence since for a system 
in eqidlibrium the free energy should be a minimum 


dFg 

322* 





or 


2)3R = M. 


Inserting this into 7.17, using 6.10, 

Fg{M) == — MR-f iMR = — iMR = 


47r2lf® €g —1 

"W2i;;+T* 


Thus from 7.16 and 7.18 


F = Fi{M)-^Fg{M) = 


2^Jf® 36, 

V (26,-hl)(€.-l)* 


(7.18) 

(7.19) 
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and henoe 'with.f 7.14, 

'W == pyiliHfsdblKfszill, (7.20) 

27r 

which, shows that equation 7.11 actually is fuUSlled identically. 
The derivation of this result demonstrates the importance of the 
long-range interaction temas—^inoluded in FJ^M )-—in equation 
7.11 . This means that in view of this interaction the fluctuations 
of a sphere embedded in its own medium are different from 
the corresponding fluctuations for a sphere m vacuum. In the 
latter case F — Fi, and 7.14 leads to 


(Ml) 

Comparing 7.20 with 7.21 shows that becomes relatively 
large for substances with hi g h dielectric constant in contrast to 
which does not depend appreciably on if tg ^ 1. 

In returning now to the further development of equation 
7.11 a method will be used which forms a generalization of the 
method employed by Kirkwood {K4) in the case of liquids con¬ 
sisting of rigid dipoles. Assume that the spherical region is 
composed of N molecules or other groups of atoms in such a way 
that each group has the same average polarization in an external 
field. In a pure liquid, for instance, each group contains one 
molecule; in a or 3 ^tal it contains all the particles of a unit cell. 
The spherical region can then be divided into N ‘ units ’ each of 
which makes the same contribution to In view of equations 
7.9-7.12 it follows that Jf® is also composed of N equal terms. 
This finds its mathematical expression in equations 7.31 and 
7.32. A more direct proof that all these terms are equal is given 
below afber equation 7.30. Each such unit contains the same 
niunher, say fc, of elementary charges which in each case can be 
arranged in a similar way relative to each other, and relative to 
their surroundings. Let r^ 2 ,-..r^fc denote the displacements 
of the k elementary changes of the jth unit, and let 


OO 00 

t Using J da? = J J dSa?. 

0 0 
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denote the whole set of displaoements of the jth unit. Then the 
set of displacements X of all the charges izx the whole spherical 
region is composed of the sets of displacements os^, 
of all the N miits. Correspondingly the volume element dX. 
(cf. 7.21) is the product of the volume elements 

dxj = dv^dTj^,.-. dTjjg 

of all the units 


dj^ — dx^dx^*** dx^* d^jc^m 

This means that in an integration over the displacements of all 
the elementary charges of the spherical region one can start 
integrating over the displacements of the charges in the first 
unit, denoted by dx-i, and so on. Alternatively, the volume 
element dX may be written as 

dX = dXjdxj, (7.22) 

where dXj = dx-^dx^...dXj^-ydXf+x‘’'d^N (7.23) 

refers to integration over the whole spherical region except the 
jth unit. 

Now let tix{Xj) be the electric dipole moment of the j’th tmit. 
Then in view of 4.7 the moment M(JC) of the whole region is com¬ 
posed as the vector sum of all the to.{x^), 

M{X) = f tn(xj). (7.24) 

i-i 

Therefore M^{X) = M(X)M(Z) = 2^m(a!^)M(X). (7.26) 

Inserting this into 7.12 leads to 

== 2 f in{xAM{X)e-^i^^^dX. (7.26) 

j-i J 

The integrations in each term of this sum will now be carried 
out in two steps, as indicated by equation 7.22. In the jth. 
term the integration will first be carried out over the whole 
spherical region except the jth unit, and subsequently over this 
unit. For the first step (volume element dXj) the quantity m(a;^) 
must be treated as a constant because it depends only on the 
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displaoements of the jth xinit. The term of 7.26 corresponding 
to this nrii t can be written as 

J J dX 

= J J m(a;^)( J dXi)dXf 

— f tn{x^)rxx*{xj)p{xj)dxp (7.27) 

where 

m*{Xi) = J (7.28) 

and (of. 7.10) 

= J e-c^CsyfcT dxJ f e- ucxvfcr dX (7.29) 

have been introduced. rx%*{xj) represents the average moment 
of the whole sphere for a fixed set x^ of displacements of the jth 
unit leading to a moment m(a:^). p{Xj) is the probability of finding 
the^th imit with this particular set of displacements. Inserting 
firom .7.27 into 7.26 3 delds 

!lf* = 2 f ta{Xj)m*{Xf)p{Xj) dXj. (7.30) 

j-i •' 

Practically all terms of this sum are equal. For m* represents 
the moment of a large spherical region polarized by one (say 
the^th) of its units whose moment is kept at a value m. Accord¬ 
ing to electrostatics (of. Appendix A 2.ii) the same moment xn’*' 
is contained in any other sphere enclosing the jth. unit (not 
neoessanly oo-oentrio), but with no further restriction of its 
position and radius so long as it is sufficiently large to be treated 
macros copically. The actual value of m* is thus determined by 
short-range interaction and is independent of the position of the 
^th unit as long as its distance from the surface is sufficiently 
large to allow its interaction with the outside to be treated on 
a macroscopio basis. The number of units for which this does not 
hold can be made very small compared with the total number N 
if the spherical r^on is sufficiently large. 

The equivalence of all N terms in the sum 7.30 means, of 
course, that all units have the same average polarization. This 
fact might indeed have been used to show the equivalence of all 
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the above terms in a much simpler way. The above oonsidera- 
tions give, however, some insight into the forces determining 
m’*'. It follows from these considerations that a region makes a 
contribution to m* only if the average moment induced in it 
by the ^*th unit cannot be obtained by treating the Jth unit as 
either a point dipole or as a polarized sphere with moment m. 
Thus it is essentially the short-range forces and the deviation of 
the shape of molecules from a sphere which makes m'*' different 
from m. 

It now follows that 

mm* = J m(aji)m* (0:1)33(0:1) dx ^ 

==.., = J m.{Xj)Tsa*{Xf)p{Xf)dXj — ... (7.31) 

and hence with 7.30, 

'W = Nvaxa*. (7.32) 

Inserting this into equation 7.11 leads to our final expression 
for the static dielectric constant, 


, 3€- ^Nn mm* 


(7.33) 


whei*e ^0 == NjV (7.34) 

is the number of units per unit volume. By its definition 7.28 
m* represents the average dipole moment of a spherical region 
embedded in its own medium, if ono of its units is kept in a given 
configuration leading to a dipole moment m. mm* is the 
average value of the luoduot mm* taking into account all 
possible configurations and weighing them according to the 
probability of finding the unit in such a configuration, m* 
differs from m because of the existence of short-range forces 
or because of the non-sphorical shape of molecules. 

£)(piation 7.3.3 is perfectly general. It will now be specialized 
by separating the contribution to duo to elastic displacements 
of electx’ons by treating them on a raaoroseopio basis. This con¬ 
tribution can be meaHuro<l with the help of the optical refractive 
index n because at t>ptioal frequencies all other contributions 
have ceased to exist in view of the higher inertia of the heavier 
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particles. If is the electric moment due to electronic dis¬ 
placement by a field C, on the assmnption that atomic nuclei 
are not displaced, then 


n^—1 = 4^MJVE. (7.36) 

Since according to Maxwell *s law the dielectric constant at 
high frequencies is equal to the square of the refractive index, 
this relation follows directly &om equation 1.9. 

The spherical region now consists of charges which follow 
the laws of statistical mechanics and which are embedded in a 
continuous medium with dielectric constant n^. Thus if M;b> is 
the moment due to aU other displacements, total 

moment of the substance. Hence equation 1.9 becomes 


or using 7.36 = 4^nMj3<iVE, (7.36) 

Instead of calculating the electronic polarization we may 
introduce it through 7.36 as an empirical quantity on the under¬ 
standing that T/t, is the optical re&active index. Nearly all of the 
above developments remain unaltered if M(J!C), m, and m* now 
refer to non-electronic displacements. Alterations have to be 
introduced only at the following two points: 

(i) In equation 7.6 the cavity field G must now be replaced by 
O', the field inside a spherical cavity with dielectric constant 
w* instead of in an empty cavity. Then according to Appendix 


A 2.14, 



(7.37) 


i.e. the denominator (2e,-|-1) in 7.6 and in the following equations 
must now be replaced by (26^-(-«*). 

(ii) In deriving 7.11 firom 7.9 equation 7.36 should be used 
instead of 1.9. Hence 7.11 will be replaced by 


_ 2 _ 47r Jlf B 

3F hT' 

All further developments remain unaltered, leading to 

€ —n* = ^ttATo mm* 

* 2e*-l-TOB 3 hT 


(7.38) 


(7.39) 
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instead of equation 7.33. It should be realized that in consider¬ 
ing the interaction energy U (X) required for the calculation of 
m and m’" (cf. equations 7.27—7.31) the electronic polarization 
has to be included on a macroscopic basis. 

Mnally, it should be noted that the diSerence between 7.39 
and 7.33 is of rather a trivial nature whenever ^ w®, i.e. when 
the electronic contributions are small. In this case 


3e, 


3e 


a 


■ 26^—j—71® 

and 7.39 becomes 


2€^-|— 1 2 


e. >71® 


€g—n‘ 


„ ,T min* 

STrJyj,- » 


€,>71 


2 


(7.40) 

(7.41) 


while 7.33 would have e,— 1 on the left-hand side, which thus 
differs from 7.41 by the optical contributions ti®—1. 

Mixtures 

Some substances contain a number of distinctly different 
groups of charges such as the different types of molecules in a 
mixture, or the positive and negative ions in ionic crystals. In 
these cases a unit cell would often be imcomfortably large or 
have an undesirable shape. Whenever these groups are well 
separated from each other it is possible to derive a formula 
for e, in which the contributions of these different groups are 
separated. Suppose that there are z different types of these 
groups. Then z different kinds of units will exist, each being 
representative for one type. If the substance contains per o.o. 

such units of the 1st, 2nd,...zth kind, then 
7.24 can be written as 


—— = 2 2 m(a;^)-|-.... (7.42) 

' /“I 

For the whole further development each group can bo treated 
separately, so that instead of 7.33 the final result is 


1 = 


3€, 47r 




2€,-H 1 3kT 2^ 


m, m* 


(7.43) 


where m* is the average dipole moment of a spherioal region of 
the substance embedded in its own medium if one unit of the 
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«th kind is kept in a configuration corresponding to a moment 
m^. is the average value of m^m^. As before, the spheii- 

oal region must be large compared 'With a region just big enough 
to have the same dielectric properties as a maorosoopio specimen. 


Summa/ry 

A short summary of the results of § 7 is given below for the 
benefit of the reader who has not studied it in detail. 

(i) The static dielectric constant e, according 'to 7.11 and 7.21 
satisfies exactly the relations 


47r 3e, Jf*_ 4n- €,+2 

3F 2€g+l ~~3~ ~W 


(7.44) 


where is 'the mean square of the spontaneous dipole moment 
of a suf5.oiently large sphere of dielectric material of volume V 
embedded in its ovm medium; i® corresponding quan¬ 
tity for a sphere in vacuum. 


(ii) If the sphere consists of components all of which on an 
average make the same contribu'tion to the polarization, equa¬ 
tion 7.44 can be developed finrther into equation 7.33 which is 
also of a very general nature. 

By treating the electronic contribution on a macroscopic 
basis 'the less general equation 7.39 follows. 


8. Special cases 

The results of the pre'vious section are contained in formula 
7.33 [also 7.11] for the dielectric constant which is valid in a 
very general way, and in the more specialized equations 7.39 
and 7. 43. Al l these formulae represent in terms of quantities 
mm* or Jkf* which refer to properties of the material in the ab¬ 
sence of a field. To calculate these quantities requires a detailed 
knowledge of the struct'ure of the substance in question and of 
'the interaction between the particles of which it is composed. 
In general such a calctdation cannot be carried out 'without the 
use of approximations. The importance of the general form'ulae 
lies in the possibility of deri'ving from it approximate formulae 
of various types, each with a clearly defined range of validity. 
It also shows that even for a restricted class of substances, such 
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as dipolar liquids, one cannot hope to find a simple formula 
representing e, in terms of a few parameters (e.g. the Onsager 
formula), which is valid over the whole stable range of these 
substances. Thus it will be found, for instance, that for dipolar 
liquids the Onsager formula should hold asymptotically at high 
temperatures only, though the deviations from it at lower 
temperatures may be very different for different substances. 

As a first step to get acquainted with the handling of the above- 
mentioned formulae one may assume that all short-range forces 
can be neglected. For spherical molecules this should according 
to § 6 lead to either the Clausius—Mossotti or the Onsager formula, 
depending on whether one considers non-polar or polar mole- ’ 
oules. The calculations proving this are carried out in the 
appendix (A 3). They demonstrate, as the main difference 
between these two cases, that for elastic displacement (non¬ 
polar molecules) the energy of a unit depends on its dipole 
moment m; hence it is shown that m® (which for negligible short- 
range interaction is equal to mm*) increases proportional to 
the absolute temperature T, and the dielectric constant becomes 
independent of the temperature. For rigid dipoles, to take the 
other extreme, the moment of a unit is equal to the dipole 
moment /x, and hence m® = /li® independent of temperature. 
Polar molecules, of course, show always both dipolar and elastic 
contributions to 

Polar liquids; Kirkwood’s formula 

We shall now proceed to consider the general case of polar 
liquids consisting of molecules with an intrinsic dipole moment 
and with a polarizability a. Liquids of this type have been 
already investigated in § 6 on the assumption that short-range 
interaction can be neglected; this led to the Onsager formula 
6.36. This assumption will no longer be made in the present 
calculation. 

It was found in § 6 that a molecule in the liquid state has a 
dipole moment which is different from that in the gaseous state 
(p^). For spherical molecules the difference is due to the polariza¬ 
tion of the molecule by the reaction field of the surroundings. 

4980.IX ^ 
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For no other case than for spherical molecules is it possible to 
calculate the ratio of the moments in any simple precise manner. 
In view of this difElcrQty it seems reasonable to treat the effect 
of the polarization of molecules in a macroscopic way, assuming 
(as in § 6) that the main contributions are due to electron dis¬ 
placement. The liquid will thus be considered as consisting of a 
continuous medium with dielectric constant (ji = optical 
refiractive index) containing dipoles with a moment (i. In this 
model a spherical molecide in vacuum consists of a sphere of the 
continuous mediiun containing a dipole (a at its centre. Accord¬ 
ing to the appendix (A 2.32) the moment of such a molecule is 



spherical molecules. 


( 8 . 1 ) 


For spherical molecules the present model permits us, there¬ 
fore, to express the moment p. of the dipoles by the moment 
of a free molecule. 

On the basis of our model the general foxmula for the dielectric 
constant is given by equation 7.39. A unit contains just one 
dipole p, and its moment is, therefore, 





( 8 . 2 ) 


Since all contributions to m* are due to dipole orientation we 
introduce u.* bv 

m* = p.* (8.3) 


as the average moment due to the dipoles of a spherical region 
if one of its dipoles is kept in a fixed direction. Equation 7.39 

contains the quantity mm*, where the bar indicates averaging 
over ail possible values of m. In the present model the only 
variable is the direction of the dipole p.. Since in a liquid all 
dipolar directions are equivalent, pp.* will have the same value 
for all these directions. Thus 


mm* = pp.* = pp*. (8.4) 

Insertmg 8.4 into 7.39 leads to the Eirkwood formula 


•s 


—n^ = 


Sgg 4Trj&?J>pp* 




( 8 . 6 ) 


In order to make use of this formula we must calculate p*. 
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To simplify this calculation it seems permissible to assume that 
short-range interaction between nearest neighbours only need 
be considered (of. first section of § 5). In this case p.* is the 
vector sum of the moment p. of the central dipole, kept in a 
fixed direction, and of the average of the sum of the moments of 
the nearest neighbours. Therefore if s; is the average number of 
nearest neighbours, 

pp.* = p*(l-|-»oosy), (8*fi) 

where cosy is the average of the cosine of the angle between neigh¬ 
bouring dipoles. This average value has to be calculated by 
appl 3 dng equation 7.28 for m’*' to the present case. Since the 
directions of the dipoles are the only variables, this equation 
reduces to 

cosy = r oosye“*^/*^da)ida >2 j J* (8.7) 

where U is the part of the energy of interaction between neigh¬ 
bouring molecules in the liquid which depends on the angle 
between their dipoles; this energy may depend on the state of 
other molecules of the liquid, and U is then the energy averaged 
over all states of the other molecules (considering their proba¬ 
bilities). dcDi and dco^ the surface elements of the solid 
angles of the directions of the two dipoles. 

The actual value of cos y depends on details of the interaction 
between the two molecules which may be very different &om 
the electrostatic interaction between point dipoles as pointed 
out in § 5; repulsion forces, chemical bond, and other types of 
interaction must be considered. It should be noted, however, 
that a large value of U does not necessarily lead to large values 
of cosy; to do so U must also have the correct symmetry. Thus 
if, for instance, U is proportional to an even power of cosy 
^-uikT jg also an even function of cosy, and cosy therefore 
vanishes, in contrast to the case in which U is an odd function 
of cos y. This means that an interaction tending to direct dipoles 
with equal probability either parallel or antiparallel does not 
play any role in the determination of pp-’". Either of the two 
types of interaction leads, however, to a restriction of firee 
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rotation of a molecule termed by Debye [D3] as bindered rota¬ 
tion. Hindered rotation, as we see, does not necessarily have an 
influence on the dielectric constant. 

The energy TJ need not only depend on the angle between the 
dipoles, but may also be a function of the direction of the vector 
joining the dipoles. Even so, we can always put 

V = !7„„+ CUa, 


where reversal of the direction of one of the dipoles does not 
alter but changes the sign of 

Now let us consider that the temperature is sufficiently high 
to make ^ (8 9 ) 

for all values of Uf^. Then 

Q-UlkT ^ g-Uwen/Jer/l (8.10) 

and hence, using 

J cosyfi“^«^*^da>idc«>2 = 0, 


one finds 


cosy ci: 


hT^ 


( 8 . 11 ) 


where 


Uq= j dco^dcoj { dto^ ( 8 . 12 ) 

since J doDj^dto^ = 0. 

The energy Uq may be positive or negative depending on 
whether the interaction tends to orient neighbouring dipoles 
anti-parallel or parallel. 

Inserting 8.6 into 8.6, we find a specialized form of the Kirk¬ 
wood formula. 


€»—= 


3e 


8 


2eg-)-w® ZkT 


(1+goosy). 


(8.13) 


In the particular case of spherical molecules ^ can be expressed 
by the vacuum moment with the help of 8.1, and hence 

3cs /w*-l-2\® pS 






ZkT 


(l-(-zoosy) 


for spherical molecules. (8.14) 
This equation differs from the Onsager formula 6.36 by the 



n,§8 


SPECIAL OASES 


63 


2 0 osy term. According to 8.11 this term tends towards zero at 
temperatures at which TcT’^ |C^|. Therefore it is foimd that 
the Onsager formula should hold in liquids at temperatures for 
which kT is large compared with the directional part 
interaction energy. This energy may have different orders of 
magnitude for different liquids; and this may lead to very 
different ranges of validity of the Onsager formula. In some 
liquids this range may fall into the whole liquid range at normal 



Fig. 6. Temporature>dependen.oe of the dielectric oonstaiit in. the high- . 
temperatiire region. The full line represents Onsager’s formula. Short-rang^ 

interaction tending to orient dipole s par allel (cos y > 0} leads to larger values 
of e, (-); the opposite case (cos y < 0) leads to smaller values (- • - -). 

pressures, whereas for others zcosy may be appreciable in the 
normal range of the liquid and may have either sign (cf. Fig. 6) 
thus making the Onsager formula inapplicable in these cases. 
In gases, in view of the large distance between molecules, 
dipolar forces only need be considered. For very small densities 
interaction may be completely neglected, leading to equation 
6.16 for €g. For higher densities Onsager’s formula should hold. 
This has been shown more directly by van Vleck [V2]. 

Dipolar solids 

The average potential energy of a polar molecule in a crystal¬ 
line solid in general depends on the direction of its dipole relative 
to the crystal axes. In liquids in contrast, though the dipoles of 
neighbouring molecules have a tendency to orient themselves 
in definite directions relative to each other, the average energy of 
a single dipole is the same in aU directions because a liquid does 
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not possess any preferential direction. The average potential 
energy of a molecule in a solid is often described as being due 
to a ciystaUiue field aotiug on the dipoles.f The crystalline 
field, being due to the interaction between molecules, normally 
depends on temperature. Usually there are several dipole direo- 
ilons for which the average energy of a molecide has a relative 
minimum; they will be called equilibrium directions. The poten¬ 
tial barrier between these equilibrium directions is in general 
very high and prevents hee rotation of the molecrdes even at 
temperatmes near the melting-point. The hypothesis that above 
a critical temperature molecules should be able to rotate tireely 
was introduced by Pauling [PJ] to explain the sudden change 
of the dielectric properties of many solids at a critical tempera¬ 
ture. This idea was a useful working hypothesis, but the present 
evidence seems to indicate that rotation does not occur in most 
solids (e.g. Z/J), Instead, the transition has to be considered as 
leading horn an ordered arrangement of dipoles to disorder. 

This transition horn order into disorder will be used as a guid¬ 
ing principle in discussing the general behaviour of the dielectric 
constant of polar solids. Its detailed properties, of course, 
depend on details of the crystalline field, but the unai’n features 
can be disou^d without such a detailed knowledge.$ To de¬ 
scribe the main features of an order-disorder transition consider, 
a sunple two-dimensional model consisting of dipoles arranged 
m a cubic face centred lattice. The crystalline field is assumed 
to le^ for each molecule to two equihbrium positions with oppo- 

^^^olute zero of temperature tho 
^ Its lowest energy state. In this state the 

ftf ordered way, but there are a number 

deSi on ordered states 

etc.). These ordered 
«t«gements oen be oJaesifled aooerdiog to whether or not they 

*«»*e w l»oa 
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lead to a residtial dipole moment of the crystal. Thus, if aJl 
dipoles are parallel (hig. 6u), the solid has a permanent dipole 
moment, but if, for instanoe, dipoles in the comers have the 
opposite direction from those at the centres (Fig. 6 6) the dipole 
moment of the crystal vanishes. In the former case the energy 



Fia. 6. Equilibrium positions of dipoles for a simple model of a orystalline 
solid. Ordered state: (a) in the ease of permanent polarization; (&) for vanishing 
polarization (two dipoles per unit oell); (c) disordered state, both, directions 

have equal probabilities. 



Fig. 7, Average potential energy of a dipole (due to short-range interaction) 
as a function of its direction for the model of Fig. 6, assuming that all the 
other dipoles stay in the fixed directions indicated in Fig. 6: (oc) for the central 
dipole of Fig, 6 6; (^) for a corner dipole of Fig. 6 & or for any dipole of Fig- 6 a; 

(y) for the disordered case. Fig. 6 c. 


of a dipole is lower when it is directed to the right than for the 
opposite directions (of. Fig. 7 j8). In the latter case the energy of 
the comer-dipoles behaves in a similar way, but for the centre- 
dipoles the two directions are exchanged (Fig. 7 a). The internal 
fields for these two types of dipoles are, therefore, different; 
they transform into each other by a reflection on a plane perpen¬ 
dicular to the dipole direction. In. the former case (a) all lattice- 
points are equivalent cmd a unit oell contains one single dipole. 
In the latter case (6) two different types of sites, comers and 
centres, have to be distinguished, and a unit coll oontams two 
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dipoles. Alternatively, the crystal may be considered as com¬ 
posed of two lattices of type (<x), but with opposite directions of 
the dipoles. 

Clearly the division of ordered states according to whether 
or not they lead to a residual moment holds quite generally 
(i.e. for three-dimensional lattices as well), although there may 
be structures leading to more than two equilibrium directions 
for each dipole. In the following the case of two equilibrium 
directions only will he considered; most results hold qualita¬ 
tively, however, for more complicated structures as well. 

Now, starting from an ordered state at the absolute zero of 
temperature, imagine the temperature to increase gradually. 
This will cause some dipoles to turn into the second equilibriiun 
position in which they have a higher energy, thus creating some 
disorder. This in turn will cause a decrease of the average energy 
difference between the two directions because the energy of 
interaction has its lowest value in the completely ordered state. 
The average energy difference of a molecule in the two equi¬ 
librium directions is thus a function of temperature, say 

V{T) > 0, 

which decreases with increasing temperature. It is found that a 
critical temperature Tq exists at which it vanishes 7 y). 

Let us now define the direction which a dipole at a given 
lattice-point has in the completely ordered state (at ST = 0) 
as the right* direction, and the opposite direction as the ‘wrong’ 
direction. Thus in !Fig. 6 a,->is the right direction; in !Fig. 
6 6, is the right direction for corner-dipoles and is the right 
direction for centre-dipoles. Let the probability of finding 
a dipole in the wrong direction, and (1— w) therefore, the pro¬ 
bability of finding it in the right direction. Since V{T) is the 

energy difference between the two positions, it follows from 
statistical mechanics that 


and hence 


uo 

——— 

1 —w ’ 


( 8 . 16 ) 


W s= 


^-vmiJcT 


1 - w r=s 




l-\-e-V(T)lkT» 


(8.16) 
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A calculation of the temperature dependence of F( T), and hence 
of w is very difficult, but approximations have been devised to 
simpliiy it. A summary of these methods can be found in 
reference We shall not give details of such calculations 



Fig. 8. Temperature-dependence of the average energy diSeronoe V{T) 

between opposite dipole directions, schematically. 



Fig. 0. Temperature-dependence of the probability w{T) of finding the dipole 

in the ‘wrong’ direction, schomatioally. 

here. It is sufficient for us to note (of. Fig. 8) that near T — 0, 
F(T) is nearly constant, but as T aj[3proaolies Tq it decreases 
rapidly towards zer<^, which it reaches i\.t T ~ T^. Thus 

V{T) = 0 if T > Ji¬ 
lt follows (cf. Fig. 9) that for low temperatures 

w ^ 1 if AjT < F(0), (8.17) 

whereas to — if y ^ J^. (8,18) 

Thus above JJ, both dipole directions have equal probabilities, 
i.e. the lattice is disordered (cf. Fig. 6 c). Calculations also show 
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that Tq is of the order of F(0)/X;; the exact value depends on 
details of the stmoture. 

The quantity (1— 2w) is a measure for the degree of order of 
the lattice because 1—=1 at jT = 0, and 1— 2w = 0 if 
T ^ 7^. This t 3 ?pe of order is ofben called long-distance order 
because it defines right and 'wrong directions for any lattice- 
point. In contrast, shoirt-distance order is the order of neigh¬ 
bours relative to each other. It means that, in view of the 
in'teraotion, the direction of a dipole is always infiluenced by the 
directions of its neighbours. £!ach dipole tends to orient itself 
in a certain direction relative to its neighbours. Long-distance 
order vanishes in the disordered sta'te, or in liq'uids. Shoirt- 
distance order persists, however, 'though it decreases wi'th 
increasing temperatiure. In fact the absolute value of cosy 
introduced above for liquids is a measure for the short-distance 
order. 

Let us now discuss the dielectric properties of ordered solids 
near T = 0. If the order is of type (a), the solid is permanently 
polarized. This case 'will not be further considered here. For 'the 
other kind of order the only contribution of the dipoles to the 
dielectric constant is due to elastic displacement of the dipole 
direction by an external field. Thus near T ~ 0 the dielectric 
constant is larger than n®, say where n is the optical 

refractive index due to elastic displacement of electrons (of. §§ 6 
and 7). The difference €„— is usually smaU. 

At higher temperatures, where dipoles are capable of changing 
their directions, increases with 'temperature, as will be shown 
presently. As in the case of liquids, contributions due to elastic 
displacements will be treated in a macroscopic way, and the 
following considerations are, therefore, based on equation 7.39. 
For the type of order which we consider at present a 'unit cell 
contains 'two dipoles as mentioned before, and at 7^ = 0 these 
dipoles have opposite direction. For the two-dimensional case 
discussed above the unit cell contains a comer- and a centre- 
dipole as shown in Fig. 6 b. The moment m of the unit cell 
depends on the directions of the dipoles. Since the potential 
barriers between the equilibrium directions are very high, it 
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seems sixffioient to consider equilibriom directions only. The 
moment m can thus have a discrete number of values mi, 

«rig . If is the probability that the unit cell has the moment 

nif, equation 7.31 applied to the present case becomes 

mm* = 2 pi, (8.19) 

where the stun extends over all values of m^. mf is the average 
moment of a sufficiently large spherical region around the unit 
cell if the moment of the latter is kept at m^. 

In our particular case of two equilibrium positions with 
opposite dipole direction the following four states exist: 


• 

% 

Configuration 


Pi 

energy 

1 


wix =*» 0 

3>i =» (1—w)* 

0 

2 


s=3 — 

=s U>(1—w) 

F(r) 

3 


m, 2 /jl 

p, = toll—Uf) 

v\t) 

4 


m4 — 0 

JJ* = to* 

2V{T) 


Heire the probability p^ was assumed to be the product of the 
probabilities of fi ruling either of the two dipoles in its respective 
direction. Thus for i = 1 both dipoles are in ‘right’ positions, 
each having the probability (1— w). This way of calculating Pf 
is not free from objections because it does not take into account 
the correlation between neighbouring dipoles giving rise to 
short-distance order. It is, however, a sufficiently good approxi¬ 
mation for our present qualitative treatment. 

With the above four states equation 8.19 becomes 


mm* = 8w(l— 


( 8 . 20 ) 


where 2yx* is the moment of a spherical region around the unit 
cell if the moment of the latter is kept at rb2|i.. 

Inserting from 8.20 into 7.39 [it should be remembered that 
this equation holds for poly-crystalline material], we find by 
adding €„—w* (see above) as contribution of elastic dipole- 
displacement 



2e,-l-w* 


ZkT 


4w{l — w). 


( 8 . 21 ) 


Here Nq is the number of dipoles per unit volume (not the number 
of unit cells as in 7.39). 



60 


STATIC DIELECTRIC CONSTANT 


n,§8 


For temperatures below the transition point, w{l — w) is 
probably the determining factor leading to an increase of c, 
with temperature. Above the transition point i.e. 

4w(l— w) — 1. Equation 8.21 then becomes nearly identical 
with the Kirkwood equation 8.6 for liquids. A difference remains 
(i) because 8.21 contains €„ instead of n® at the left-hand side, 
and (ii) because fx* is dejSned differently, f 

A more exact calculation can be carried out, however, above 



Fio. 10. Temperature-depezideace of the dieleotrio constant of a dipolar solid 

showing am order^disordeT transition» scheznaticailly. 

Tq since the unit ceU then contains only one dipole. This calcula¬ 
tion follows the lines given above for liquids and leads exactly 
to the Kirkwood formula 8.6. The difference from 8.21 should 
thus be due to the approximations made in its derivation. 

Thus, the dielectric constant of a solid (of. Fig. 10) which has no 
permanent polarization shoidd start at T = 0 with a value 
which is slightly larger than w®, and should rise very slowly at 
first and more rapidly as T approaches T^, provided Tq is below 
the melting-point. Above T^, Cg should decrease with T and 
there should be no appreciable variation at the melting-point. 
If, on the other hand, cannot be reached below the melting- 
point, then Cg increases up to the melting-point, and then de¬ 
creases. 

The similarity of the dielectric behaviour of disordered solids 

t Kirkwood: (i* ia the moment of a spherical region if one of its dipoles is 
kdpt ia a fixed, direction p/ft. Squatioa 8i21: 2^* is the correspoadiog monicnt 
if two dipoles (of one unit cell) ore kept parallel with a total moment 2|i. 
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and of liquids is not surprising. In faot we have seen before 
(§§ 4 and 6) that a dipole having two equilibrium positions with 
equal energy behaves similarly to a dipole with a continuous 
range of such positions. 

Finally, it shbuld be remarked that as for liquids the quantity 
is a measure for the short-distance or^er. Above the 
transition temperature Tq this quantity can be found experimen¬ 
tally if €g is known over a sufficiently large temperature range. 
Near Tq it may be expected to depend sensitively on the structure. 

Intramolecular rotation 

Consider briefly the case of a liquid consisting of molecules 
in which intramolecular rotation is of importance. Again the 
unit cell contains one single molecule, but its dipole moment 
may depend on the set of coordinates x describing intra¬ 
molecular rotation. Thus we put m(ir) = (i.(a;), and introduce 
correspondingly m*(a:) = [!.*(«) as moment of spherical region, 
In the Kirkwood formula 8.6 the term pp.* must then be replaced 
by the average p(a:)p*(a;). If in particular the molecule has a 
number of semi-stable states with moments pg,..., besides the 

ground state whose moment is po, then similar to 8.19 

H{a:)li(*)» = (8.22) 

i>0 

where represents the probability of flnding the molecule in 
the state i. 
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9. The establishment of equilibrium 

Thu present chapter will be devoted to an investigation of the 

dynamio properties of dielectrics. This will be found to be a 

much more difficult task than the development of the theory 

of the static properties, and in fact it will be seen that only for 

dflute solutions of dipolar molecules has it been possible so far to 

carry out quantitative calculations. The reason for the greater 

difficulties in dealing with the dynamic properties is evident if 

we remember that for the static case it was not necessary to 

investigate the kinetic properties of molecules. In the present 

section a qualitative discussion of the problems involved will be 
given. 

It *will be remembered that in § 2 it was found that in altemat- 
ing electric fields the electric displacement shows a phase-shift 
with respect to the eleotnc field. This leads to the introduction 
of two dielectric constants € 3 ^ and € 2 , both of which depend on the 
frequency <of2n and which according to 2.8 can be considered as 
real and unaginary components of a complex dielectric con¬ 
stant €(«>). The two quantities and €2 are not independent of 
each other, but can both be derived from a function a(i) of time 
with the help of equations 2,14 and 2,16. Making use of 2 . 8 , 
these equations can be written as a single complex equation, 

00 

= €«+ J a(a:)e»*“* dx, ( 9 . 1 ) 

0 

where is a variable of integration. 

The function a(f), according to § 2 , describes the decay of the 
polarization of a dielectric with time if the external field is 
suddenly removed. Alternatively the gradual increase of polari¬ 
zation Tidth time to its equilibrium value if the dielectric is 
brought into a constant field can also be described with the help 
of the decay function a{t). It is thus seen that with the use of 
the macroscopic equation 9.1 the frequency dependence of the 
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complex dieleotrio constant (and hence of dielectric loss, of. 3.16) 
is imiquely connected with the way in which equilibritim is estab¬ 
lished In a dieleotrio brought into a constant field. We shall 
choose this latter process for the qualitative considerations 
which follow because it is very suitable for the purpose of showing 
in a simple way the difficulties involved. 

As in § 4, consider the two characteristic types of displace¬ 
ment of charges, namely (i) elastic displacement, and (ii) dis¬ 
placement to another equilibrium position, and assume that all 
interaction forces between particles can be neglected. From § 4 
we then know the average position of a charge in the presence 
and in the absence of a field, or rather its average displacement 
by the field. It is now ova task to consider the infiuenoe of the 
field in more detail by assuming that the dieleotrio consists of 
such an assembly of (non-interacting) charges which, in the 
absence of a field, are in theimal equilibrixun. 


Case (i). Elastic displacement 

Fach particle of mass m and charge e is bound elastically to its 
equilibrium position and will carry out harmonic oscillations 
with frequency ojg about it. If r is the displacement, then in the 
absence of a field 

(9.2) 


d^r 


= —cugr. 


and hence r = GoOOs(a>oit-|-8o), (9.3) 

where the maximum amplitude Gq and the phase Sq are indepen¬ 
dent of time, and the energy is given by imojQ G%. If the motion 
is not perturbed, the oscillator will' keep this same energy in¬ 
definitely. It should be realized that this is in contradiction to 
the postulates of statistical mechanics, according to which 
within a sufficiently long interval of time the oscillator should 
be found with various energies according to the Boltzmann 
theorem, and its average energy at a temperature T should be 
kT. To achieve this some kind of interaction with other particles 
or with the surrounding medium—permitting an exchange of 
energy—must be assumed. Frequently it is assmned that this 
interaction takes the form of collisions of extremely short dura¬ 
tion. This means that the equation of motion 9.2 holds except 
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during a oollision. Theitefore the solution 9.3 holds between two 
collisions, but at a collision both the amplitude Oq and the phase 
So usually change their value. The assumption of collisions of 
very short duration means that each oscillator actually satisfies 
the equation of motion obtained by neglecting all interaction 
terms, but as a result of a collision it makes a transition into a 
state with diff erent energy (and phase). This, as may be shown 
by statistical mechanics, leads to the correct average energy if 
averaged over a time which is long compared "with the time 
between two collisions. It is of great importance that in equili¬ 
brium this average is independent of the nature of the poUisions. 
Equilibrium properties can, therefore, be derived without con¬ 
sideration of collisions—quite in contrast to the question of 
attaining equilibrium which we intend to investigate presently. 
It should be remarked first that the collisions need not occur 
between the oscillators themselves as one would expect in a gas 
of such oscillators. They might be considered as dissolved in 
another medium (liquid or solid) which is in thermal equilibrium, 
and the collisions to be considered occur between the oscillators 
and the particles of the medium. 

Assume now a constant field f to be applied at the time i = 0. 
Then the displacements satisfy the new equation of motion 4.8 
instead of 9.2. A particle now oscillates about a new position of 
equilibrium displaced by a vector r relative to the former one, 
but the motion is still harmonic (of. 4.9 and 4.10). Both the 
nr>a.ximun> amplitude C and the phase 8 usually change, however, 
at 2 = 0, although r remains continuous, as is showm in Tig. 11. 
This assumes, of course, that no coUisions occur during the 
represented time interval. 

To obtain the polarization P of a dielectric substance con¬ 
sisting of a large number of oscillators the vector sum of all 
displacements r has to be formed according to 4.6. At a given 
instant of time f c 0 in the absence of a field the phase-angle Sq 
can be assumed to have any value between 0 and 27r wdth equal 
probabilities. The vector sum of all displacements, and hence 
the total electric moment of the dielectric, will, therefore, vanish 
at any instant of time. Now let us assume as a first approach 
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that no oollisioiis occur after the field has been applied at i = O. 
Then since all individual displacements vary continuously the 
same must hold of the polarization P. On the other hand, the 



(«) 



Fia. 11. Timo dependonoe of the displacement r of a harmonic oscillator if a 
constant field / is applied at the time i = 0. The displacement and. its timo 
derivative remain continuous at t = 0, but the average displacement is 
altered. Two characteristic cases are shown: (a) both phases before and after 
application of the field vanish; the maximium amplitude Q for t 0 is smaller 
than it was for i < 0; (6) both phases ai»e equal to —w and the maximum 

amplitude is larger for ^ > 0 than for t <z 0. 

time average of the displacement of any single oscillator is equal 
to r so that the average polarization too is different from zero. 
It will therefore oscillate about its average value with the 
frequency coJ27r. This means that after application of the field 
the oscillators are no longer distributed uniformly over all values 
of the phase 8, as can also be shown by direct calculation. 

Collisions, as has been pointed out above, alter the phases. 

4980*11 
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They will tend to suppress the oscillations of the polarization 
about its equilibrium value. One should expect, therefore, that 
after application of a held the equilibrium polarization is reached 
steadily if the avera>ge time between two collisions of an oscillator 



0 b-»- 

( 6 ) 



Fia. 12. Time dependence of tlie polarization of a dielectrio if a £eld is applied 
at the time f = 0. (a) and {b), the material consists of dipolar oscillators^ 
(a) without, (6) with collisions, leading to oscillations of the polarization about 
its average, (c) The material consists of rigid dipoles with several equilibrium 

positions. 


is shorter than its period. Otherwise transient oscillations of 
the polarization will be excited which last for a time of the order 
of the time between two collisions (cf. Kg. 12 a, b, and c). 

Case (ii). Displacement ivio another equilibrium position 
As in § 4 (of. Kg. 4), a charged particle is assumed to possess 
two eqtiilibrium positions A and J5 at distance b from each other 
and separated by a potential barrier. In the absence of an 
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external field the potential energy of the particle is the same in 
A and B. In the presence of a field f according to 4.12 thei-e is a 
difference in the potential energies at A and B. In an assembly 
of a great number of particles of this kind in the absence of a 
field an equal number will oscUlate about each of the two 
equilibrium positions with an average energy hT, if thermal 
eqxiilibrium is assumed to exist, and if the height of the potential 
barrier is large so that the potential energy jfiT of a particle on top 
of it is large compared with kT, 

H^kT. ( 9 . 4 ) 

This means that the fraction of particles with sufficient energy 
to go over the top of the potential barrier is extremely small. 
lt.s order of magnitude is given by the Boltzmann factor 

oxp(— HjkT). 

Assume again that at the time t — 0 an electric field f is 
applied. Then in the absence of collisions the number of particles 
oscillating about A and about B is not altered, because the only 
action of the field is to alter slightly the equilibrium position— 
as pointed out above, case (i). The field is not able, however, to 
lift a particle over the potential barrier (cf. § 4). On the other 
hand, in equilibrium the number of particles near A is larger by 
a fraction of the order ebfjkT (independent of H) than the num¬ 
ber of particles near B, as was shown in § 4. In the absence of 
eolUsions, therefore, equilibrium cannot be established. This 
was also found to be so in case (i), but in the present case lack of 
equilibrium has a much more senous effect. For in adding all 
the displacements in order to calculate the polarization, lack 
of collisions meant in case (i) that the polarization was 
oscillating about its e<iuilibrium value. In the present case, 
however, the required displacements cannot be carried out 
without collisions, which means complete lack of polarization of 
the required type. 

Collisions will tend to establish equilibrium. The time re¬ 
quired for this process depends on details of the model. It is often 
assumed that collisions on either side of the potential barrier 
are very frequent, so that the particles on each .side may alw'ays 



68 


DYNAMIC PROPERTIES 


III, §9 


be considered to be in. equilibrium amongst themselves. By this it 
is meant that a particle having sufficient energy to flow over the 
potential barrier to the other side will suffer so many collisions 
that it is very unlikely to oscillate back to the side from which 
it came originally because on an average over many collisions 
a particle has an energy kT H. 

Starting now from equilibrium in the absence of a field, the 
immediate effect of the application of a field in the JL—B direc¬ 
tion is to lift the potential near A by the amount efb. As a con¬ 
sequence the fraction of particles with sufficient energy to move 
over the potential barrier is approximately given by 

for the A B and the B ~^A directions respectively,'}’ because 
the height of the potential barrier measured from A is now only 
H —ebf. Thus if oi^f^Tr is the frequency of oscillation of a particle, 
the probability per second for the transfer of a particle from B 
to j 4 is given by 

== (9.6) 


and for the transition A B considering 4.18 it is 




'12 


2in 


CO 
^TT 




(9.6) 


Thus if at any instant of time there is a number of particles Ny{t) 
at A and a number N^{t) at jB, a number will flow per 

second from A to JS and a number from J5 to A. Therefore 

the rates of change of N-^ and respectively are given by 


dN^ _ 

dt ~ 


= —NlWi2-\-N^W 


21 


(9.7) 


1 


dt 




dt ‘ 


(9.8) 


t By introducing a difEerent normalization of the potential of the external 

field the two exponents can be written in a more symmetrical way as 77_Je (bf) 

and 77-h^(bf), leaving the results unaltered. 
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It follows that the total number N of particles, 

N = iVi+iSTa, 

is independent of time, as required, because 

dt dt dt 

Subtracting 9.7 from 9.8 with the use of 9.9 3 rields 
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(9.9) 

(9.10) 


(9.11) 


Now according to 9.6 and 4.18, 


^la+^ai = 2uj 


ai 


'SI 


^ ^2 kT) 


2w 


ai> 


ty 


la 


efb 

^ai — 


Using these two equations, 9.11 becomes 

1 d 




1 ef b 


2w2x dt 
If we assume 

•A^i(O) == -Z^a(O) = hN 

at f = 0, 9.14 is solved by 

N f>fb 


2 kT 


N. 


(9.12) 

(9.13) 


(9.14) 


(9.15) 


The induced polarization is proportional to and thus 

approaches its eciuilibrium value exponentially, as shown in 
Fig. 12c. 

In contrast to case (i), equilibrium is reached in the present 
case by the influence of the fleld on the probability of transfer 
of a particle from to .B and B to j 4 respectively, according to 
equations 9.5 and 9.6 or 9.12 and 9.13. This model can be 
generalized by introducing transition probabilities tOjg and lOgi 
measuring the probabilities per second for the transition of a 
particle from A to B and from B to A respectively, without 
specifying them by equations 9.5 and 0.6. This means that 
equation 9.11 holds as before. Instead of calculating 
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Wai from a model, and then deriving 9.12 and 9.13, these equa¬ 
tions can also he obtained from 9.11 by considering solutions 
corresponding to equilibrium which can be obtained from § 4 
without the knowledge of the transition probabilities. Thus 
since in equilibrium — N^)jdt = 0, in the absence of a field 

Ni — N^, and hence from 9,11 = w^. In the presence of a 

field, using 4.13, 4.12, and 4.18, 

= = (9.16) 

This must be a solution of 9.11 if N^—N^ is independent of time, 
which leads to 9.12 and 9.13. 

The model discussed above can be further generalized by 
considering dipoles which can have several equilibrium positions 
with certain transition probabilities between them. In equi¬ 
librium in the absence of a field the number of transitions from 
a given position is just balanced by the number of transitions 
from other positions to it. An electric field alters the transition 
probabilities, and hence, in equilibrium, the distribution of 
dipoles over the various positions. The time required to attain 
equilibrium depends on the transition probabilities. 

The preceding considerations have thus shown that in the case 
of elastic binding the field displaces the charges which then 
oscillate about their equilibrium positions. In the case in which 
charged particles or dipoles possess several equilibrium positions 
the field does not, by immediate action on the charges, transfei* 
them to their new positions, but it alters the transition pro¬ 
babilities between them. This in turn leads to the establishment 
of equilibrium. 

10. The Debye equations 

In the present section it is intended to derive equations for the 
frequency dependence of the complex dielectric constant €(aj) 
which should hold for the case of dilute solutions of dipoles in 
liquids and solids and for a few other cases. These equations 
were first established by Debye [D2] and subsequently have 
been applied to many substances, not always, unfortunately, 
with the necessary discrimination with respect to the intended 
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range of validity. We shall base the considerations of this 
section on the hypothesis that in constant external fields equi- 
librimn is attained exponentially with time, as corresponds to 
case (ii) of § 9. We thus assume for the decay function a(f), 

a{t) oc (10.1)' 


where r is independent of time but may depend on temperature. 
That this assumption leads to the required properties can easily 
be shown from the relationship 2.12 between the electric field 
E{t) and the electric displacement D{t), both of which may de¬ 
pend on time. In 2.12 it was assumed that both D and E vanish 
for times f <; 0. If this is no longer the case, then 2.12 has to be 
replaced by 


JD{t) = €„E(t)-{- j E(u)(x{t — u)du. 


( 10 . 2 ) 


— 00 


This integral equation can easily be transformed into a differen¬ 
tial equation. IFor on differentiating 10.2 with respect to the 
time and making use of 



(10.3) 


which follows from 10.1, we find after multiplication by t 

t 

” eooT^^^-{-Ta(0)i7(i)— r E{u)(x{t — u)du. 

— oo 

(10.4) 

Adding 10.2 and 10.4 yields 

r^(I)-€^E)+{D-€^E) - Ta(0)i?. (10.6) 

at 

To determine the constant a(0) consider the special case of equi¬ 
librium in a constant field. This means that 

^ (£>-€„.&)== 0, D = €,E, 

and hence from 10.6 

Toc(O) = —e*,- (10.6) 
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(10.7) 

is found, as tits differential equation connecting J){t) with JE}(t) 

on the assumption that the decay function ocit") is given bv ^of 
10.1 and 10.6) 


Oi{t) = 


■8 


= co 




( 10 . 8 ) 


We shaJl now use equation 10.7 in the investigation of the 
approach to equilibrium of a condenser. The following two 
oases are to be considered: 

(a) Constant charge on the condenser plates. Then 

dD 


dt 


= 0, D=D., 


and hence, ueing 10.7. 


.dn jq, 


i.e. 2^— €gJS/oce-^', 


where 


= 00 


(10.9) 

( 10 . 10 ) 


«r 

(b) Constant voltage at the condenser plates, i.e. 

dB 


dt 


= 0 , E == B^. 


It follows with 10.7 that 

dD ^ ^ 

^ dt ~ oc c-'/’’, (10.11) 

B^h aus M to exponential approach to equilibrium. 

In ^^c flel^ aeaume £ to be repreeented by equation 

i.8, 1 . 0 . .»oc exp(-,urf). Then, introducing the complex dl- 
clectec constant c (cf. 2.8). we find, ueing 2 10 ^ 


dt 




Introducing this into 10.7 leads to 


( 10 . 12 ) 


e(a>) —€co = 


««—e 


00 


1—W 


(10,13) 
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An alternative way of deriving this equation is to insert lO.S 
into 9.1, 

on 

€M— Coo = (e«— J da:, (10.14) 

0 

which after integration result® in equation 10. Hi. 


Sepaiuting the real and imaginary parts in 10.13 according to 
2.8 we find 


ei(aj) €ao ^ la 2> 

(10.16) 


(10.10) 

and for the loss-angle using 2.6, 



(10.17) 


•Equations 10.15-10.17 (also 10.13) will bo denoted as the 
Debye eqviations and the constant t will b<^ called th<» relaxation 
time. They describe the properti<iH of a dielectric substance in 
alternating fields on the assumption of an exponential decay 
function a(i) (cf. lO.S). S<»me models leading to such a decay 
function will be studiiHl in § 11. Most of them retpiin^ 


<--i: 1. (10.IS) 

a condition which normally is fuliilhul in dilute solutions only. 

In discussing tlu' pnqx'Hh'S of the I)<'hy<’ ecpiatituis it shoidd 
ho noted that the di(i<M'tric constants i, and c .j dtqx'iid on ai- haist 
two parainet<u’s, tin' angtdar fre(|ucn(iy m and llu' t.<un]>crature 
T. The fi'iaiu<'n('y dep«Mid<'n(X’ is {'.\pn\s.s<‘d t'xplicitly, Init l.he 
temperature! a[)pears iin]i]i<Mtly through and r, both of 

which usually <hq>end on 'I’la^v may <hq)(utd <m otlu'r para¬ 
meters as well, variations of which will nol. la* (^onsid<<n'd hen\ 
For the following it will Im assn hum I Ihat hot li i, aiul» „ Jire known 
as functions of 7’. N't wi»re known as well, a new variable 


Z lt»gOIT logoi I log T (1(1.10) 

could he introdiKX'd in t<wins of which 10.1.5 ami )o. 10 heconu^ 
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Kg. 13 shows these functions. It will be noted that 
is a symmetrical function of 2 . 

Actually r cannot be assumed to be known but must be 
determined from measurements of e-i and at various &equenoieB 
and temperatures. Asstuning that the Debye equations are 
fulfilled t(T) can, however, easily be found from the frequency 



Fia. 13. The Debye fiinotionB «i/(€,—(full Une), and (ej — «aci}/(e« — «oo) 

(dotted line), aooordin^ to equations 10.20. 


at which has its maximum. In fact at constant temperature 
the angular frequency to^ of this maximum is determined by 


= 0, if at = <o^ and T = constant. (10.21) 
d<o 

Hence using 10.16, 

( 10 - 23 ) 

Inserting this value into 10.16—10.17 we find for the dielectric 
constants and the loss angle 


= i(«.+0» ^2 == e«,), tan^ = 5*- fs , 

if to == coj^. (10.23) 

Thus £ram the fiequency at which has its maximum we find 
T, and fi:om its height we ought to obtain e,—Coo if ifi® Debye 
equations are fulfilled. 

Instead of considering the maximum of eg it is often preferred 
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to determine the angular frequency at which the loss-angle 
<f> has its maximum. We then require 


dcu 

Oi = 

and T = constant. 

(10.24) 

Using 10.17, this 3 deld 8 

1 

Oi^ = - 

T A 

/-• 

/ ^00 

(10.26) 


In view of 10.18, a>^ and will be nearly equal for most sub¬ 
stances for which the Debye equations can be expected to hold. 
Inserting 10.25 into 10.16—10.17 we find 



= 2 


^8 ^00 




00 



^8 ^00 


^ 00)3 


tan^ = 


^8 ^00 

2V(e8 €„) ’ 



(10.26) 


It is an interesting feature of equations 10.23 and. 10.26 
that the values of and £2 O'f frequency or to^ at which 
either or tan ^ has its maximum are indejiendent of this fre¬ 
quency and of the relaxation time, and are expressible by the 
static dielectric constant and the high-frequency dielectric 
constant 

We shall use tliis occasion to give a more precise definition of 
e,„ than that given previously. According to the Debye equa¬ 
tions decreases from €„ to iu the frequency region in which 
€a has relatively large values. Therefore €00 is the value which is 
asymptotically approached by ei(cu) at frequencies sufficiently 
larger than a>.„t to make €g(cu) relatively small. This value 
of €„ need not coincide with the optical dielectric constant 
because most substances absorb in the infra-red region (see 
also § 13). 

If —€00 is very small, as is required by 10.18, then e,, € 00 , and 
Cl are nearly equal, 

Ca ci 6i(a>), (10.27) 


and very precise measurements would be required to obtain 
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€g—€oo- This quantity can, however, be obtained with the help 
of the phenomenologioal relation 2.18 if 

62 (c(>) = 6 ]^ tan <f> Cg tan <f> 

is known for aU fcequenoies for which it has appreciable values. 
For on inserting firom 10.16, equation 2.18 becomes 


oo 



0 


which is an identity because 


00 

^ \ 2 r 

(e.—fioo)- J > 

® (10.28) 



WT dw _ TT 

1+ojVa "Zr “ 2' 


In the particular case of dilute solutions of dipolar molecules 
in non-polar liquids e,—e® can be inserted from 6.24. In this 
case the dielectric constant of the solvent is nearly equal to 
with 10.27, 

Co — «co ei(co) €g. (10.29) 

Thus from 6.24, 10.16, 10.17, and 10.29 


6 a tan ^ a=L 


2(ea—1)(6,— 

ZkT \ 3 j\ (26a+n2)(6a+2)/ 


for spherical molecules and — 6 oo 1 . 


ayr 

(10.30) 


Here n is the refractive index of a pure liquid of the dissolved 
molecules. The factor containing n is not of great importance 
because { usually differs from unity by a few per cent, 

only. 

In the above discussion of the Debye equations we have been 
mainly concerned with the frequency dependence of the dieleo- 
tno constants using the temperature as a parameter. In practice 

and €2 are often measured as functions of the temperature jT 
with the frequency as a parameter. This does not permit such an 
immediate comparison with the Debye equations because the 
latter do not depend on T explicitly. Assuming that €q — 
is known as a function of T, it is advisable in tliis case to plot the 
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function e 2 /(®«—^oo) against T. Its maximum is determined by 
(cf. 10.16) 

^ _ 0 / €2 \ _ 0 COT _ 0 / COT \ cir . _ , 

0T\€g-e<o/ 1 + CoV® 0T\l-f-CO®T®/dT'’ 

as before 10 . 22 . This means that the maximum for a given 
frequency occurs at a temperature at which 

t{T^) = l/o,. (10.31) 

In the case that and €2 are measured as functions of co for a 
given temperature T, the Debye equations can be checked 
immediately once t(T) has been found from the position of the 
maximum of according to 10.22. For t does not depend on 
«o, and its value can thus be inserted into equations 10.16—10.17 
which then contain no unknown quantities. If, on the other 
hand, and ^2 are measured for a given frequency as functions 
of T, the value found for t{T^ must not be introduced into 
10.16—10.17 because t varies with T. Now if 10.18 is fulfilled 
and thus an exact measurement of is difficult, knowledge 

of the temperature dependence of at a given frequency does 
not permit a check of the Debye equations. If, however, 
and € 1—600 are known as functions of the temperature, cut can 
be derived from 10.16, i.e. 

(a,x)a = iaUfL, (10.32) 

€i — €„ 

and then inserted into 10,16, resulting in 

eg = V{( 6 g—ei)(ei— 6 „)}. (10.33) 

Equations 10.32 and 10.33 are equivalent to the Debye equations 
10.16 and 10.16 in their original form. They have the advantage 
of showing clearly that the relaxation time t is a quantity which 
can be calculated from the measurable quantities according to 
10.32, whereas 10.33 is a relation between measurable quantities 
only. Their disadvantage is that they require measurements of 
€g —€i and —Coo which in the cases where the Debye equations 
are supposed to hold are very small (cf. 10.18 and 10.27). When¬ 
ever this is possible, however, 10.33 affords a very simple test 
of the Debye equations whether temperature, frequency, or both 
are varied. 
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11. Models for the Debye equations 

In the previous section it was shown that an exponential 
decay function leads to the Debye equations. In the present 
section, models for which these equations hold will be studied 
and expressions for the relaxation time will be obtained. A very 
simple model has already been discussed in case (ii) of § 9 and 
§ 4. In this model the dieleotric material contains an assembly 
of charged particles whose interaction can be neglected. Each 
charge has two equilibrium positions separated by a high 
potential barrier. It is then assumed that each particle collides 
with the surrounding medium and that the number of collisions 
per second is such that the average time tq between two collisions 
is small compared witli the average time -r which a particle 
spends near one of its equilibrium positions before jumping to 
the other one. Tliis leads to the linear differential equation 9.1 I 
for the difference of the number of particles occupying the two 
positions. Its solution 9.16 approaches its equilibrium value 
exponentially with a relaxation time (cf. 9.6) 

T = ^ = H>kT, (11.1) 

IWgy a>o 

where H is the height of the potential barrier and the 

frequency of oscillation about either of the equilibrium positions. 
For this model, according to 9.6 and 9.6, the transition pro¬ 
babilities tOjg and of a particle between the two equilibrium 
positions are equal in the absence of a field, but are altered 
slightly when an external field is applied. 

It should be pointed out that from the assumptions made it 
follows that the derived exponential law can be true only on an 
average over time intervals wliich are large compared witli 
Therefore the derivation of the Debye equations from this law 
can hold only if the period of the field is large compared with to- 
It must, therefore, be concluded that this model leads to the 
Debye equations only when the foUowing conditions hold: 

To < l/a>. (11.2) 

^us, m order to ascertain whether or not the Debye equations 
hold in the mam region of absorption (near oit / 1 ). it would 
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be necessary to determine tq, and this in turn would require a 
detailed knowledge of the interaction between the charges and 
their surroundings. It will also be seen from conditions 11.2 
that even if the Debye equations do hold in this main region 
deviations may be expected at higher frequencies. 

It is important to realize that an exponential decay function 
(and hence the Debye equations) can only be obtained when the 
particles can be considered as independent of each other. Inter¬ 
action between particles (which is to be distinguished from that 
between a particle and its surroundings) would mean that the 
transition probabilities are not constant, but depend on the 
positions of the neighbours. Hence equations 9.7 and 9.8 would 
be no longer linear and, therefore, could not be solved by ex¬ 
ponential functions. 

Dipolar solids 

Models more directly applicable to solids, but leading to a 
mathematical treatment similar to the above, can easily be 
devised. Consider, for instance, the model for dipolar solids 
used in § 8. It consists of dipolar molecules, each of which, owing 
to the crystalline field, has a number of equilibrium positions 
with different dipole directions, which are separated by potential 
barriers. In the simplest case only two equilibrium positions 
with opposite dipole directions exist. In such a model at low 
temperatures the dipoles, because of interaction with each 
other, form an ordered arrangement. At a temperature Tq an 
order-disorder transition occurs and for T > Tq long-distance 
order vanishes. Short-distance order, i.e. order relative to 
neighbours persists, however, and can be neglected only at still 
higher temperatures. It will be shown that at these temperatures 
the Debye equations hold. They should be invalid at tempera¬ 
tures near Tq, but become valid again \i T <^Tq. 

Let us consider first the high-temperature region. Then in 
the absence of a field the lowest energy-lev el of a dipole is the 
same for both equilibrium directions. To carry out a transition 
between them a minimum energy, say H, is required to lift 
the molecule over the potential barrier separating the two 
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equilibrium positions. In the calculation of the probability of 
a transition between the two positions equation 9.6 cannot 
now be applied because in general a molecule is also capable of 
internal excitation, a fact which was not considered in§ 9. We 
shall take account of this by introducing a factor A which only 
varies slowly with temperature. Then if 1/2t is the transition 
probability, 11.1 will be replaced by 


2o>a 


(11.3) 


where Trlcog^ is the average time required by an excited molecule 
to turn £rom one equilibrium direction to the other. In the 
presence of a field f the molecule has an additional energy +pi 
and —pf for the two equilibrium directions of its dipole respec¬ 
tively. The minimum energy required by the molecule to enable 
it to carry out a transition in the presence of a field is, therefore. 
Thus, assuming 

fifjkT <1, (11.4) 

the two transition probabilities are given by 


= 

— 


JL JL/i . 

2t ~ 2tV ^kTf’ 
2t 2t\ kT/ 


( 11 . 6 ) 

( 11 . 6 ) 


respectively. Hence if and are the numbers of dipoles in 
the two directions, equations 9.7,9.8, and 9.11 still hold, provided 
these new values of and used. Therefore inserting 

11.5 and 11.6 into 9.11 yields 




For constant fields this equation leads to exponential approach 
to equilibrium and hence according to § 10 to the Debye equa¬ 
tions with a relaxation time r. 

Assuming that a molecule interchanges energy very quickly 
with its surroundings, it can be considered as passing through 
a number of levels above the energy H during the transition 
between the two equilibrium positions. A quasi-thermal 
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equilibrium of the excited molecule with its surroundings is 
therefore established. Under these conditions the probability 
of finding a molecule with an energy above H is approximately 
(2^/I>0)exp(— HjkT), where is the number of energy-levels 
in a range of the order JeT near the groimd level and is the 
corresponding number in the excited state above H. !lTow the 
transition probability 1/2 t must be equal to this quantity 
multiplied by ‘tt/coq, which is the average time required by an 
excited molecule to turn by an angle tt, and divided by 2 because 
only half of the molecules move in the required direction. Hence 
by comparison with equation 11.3 it is found that approximately 

^ = (n.8) 

Expression 11.3 for the relaxation time is sometimes (e.g. 
Frank, F2', Kauzmann, Kl) related to similar expressions for the 
rate of unimolecular chemical reactions. In general a formula of 
this type will be obtained for any process which requires excita¬ 
tion to an energy H, but calculation of the absolute value of a 
rate of reaction (reqmring the factor Afa>^) is difficult. The first 
calculation of this kind is due to Pelzer and Wigner [F^]* Usually 
can be assumed to be of the order 

(Oa 10^®—10^* per second, (11"9) 

BO that measurement of the temperature-dependence of t leads to 
semi-empirical values for the order of magnitude of .d (11.8). 
Sometimes (e.g. Eyring, JE/J, JE2) it is assumed that the trans¬ 
lational or rotational motion of the excited molecule can be 
separated from the other types of motion. In this case if 
< kT, there are about kTjUuif^ energy-levels, in an interval 
kTf connected with the rotation of the molecule as a whole. 
Denoting by the number of energy-levels due to internal 
excitation of the molecule, we then obtain, using 11.8, 



A ^ Do n 

kT’ 


( 11 . 10 ) 


In view of the assumptions that have been made in the deriva¬ 
tion of this expression one should not expect it to give more than 
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an order of magnitude. An d if >> TcT, a result might be 
obtained whioh is not even of the correct order (cf. Pelzer, 
P3). 

Let us now consider our model at lower temperatures. In the 
neighbourhood of the transition temperature Tq interaction 
between dipoles becomes important. The eneigy of a molecule 
then depends on the dipole directions of its neighbours whioh 
invalidates the assumptions whioh led to the Debye equations. 
It will be shown presently, however, that this is no longer the 
case when the temperature is well below the transition point 
and the dipoles are therefore in an ordered state. As pointed 
out in § 8, the lattice then contains two types of sites; and if the 
state is completely ordered, these sites are occupied by molecules 
•vrith opposite dipole directions. They maiy form, for instance, 
the centres and the comers of a body-centred cubic lattice. Slaoh 
molecule has a second equilibrium position with opposite dipole 
direction in whioh its average energy is higher by an amount 
V{T) than in the origmal position. If the temperature is sufB.- 
oieutly low (well below 2),), it is possible to replace V(T) by V (0), 
the. value of V when T = 0. .It is in this region that we can 
assume the energy of a dipole to depend on its own direction 
only, and not on the directions of its neighbours, because the 
latter are nearly always in the state of lowest energy (ground 
state). Hence we expect the Debye equations to hold. 

To prove this it is sufficient to show that equilibrium is 
approached exponentially with time since, as shown in § 10, 
the Debye equations then obtain. It will simplify the calcula¬ 
tions if we assume at a given time the existence of a dipole 
moment of the substance in the absence of a held and calculate 
ii® variation with time. Let us distinguish quantities referring 
to the two types of site by the use of -f- and — indices. Thus, 
for instance, is the probability per second for a transition of 
®' dipole on a -f- site &om the ‘ 1 ’ into the ‘ 2 ’ direction. Similarly 
N}' and N-^_ are the numbers of dipoles in the * 1 ’-direction for the 

two types of site. The total dipole moment is proportional to the 
quantity 

^ = Nt-Nt+Nr~N^. 


(11.11) 
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Also if ^ is the total number of dipoles 

= Nr-\-N^ = (U.12) 

Furthermore, in the absence of a field the ‘1’ (‘2’)-direction 

plays the same role for +Mtes as the ‘2’ (‘ 1 ’)-dixeotion for 
—sites. Thus in equilibrium at T — 0, and 

Nx = ^2 = 0. Also 

= ^ia- (11.13) 

Thus by the same method as that used in deriving equations 9.7 
and 9.8 we now find 


dNj- 

dt 





(11.14) 


from which, with the help of 11.12, we obtain 

= (11.16) 

The same equation holds for the — sites, if we exchange + for 
— indices. Using equation 11.13 this means that 

== —(«?i^+wi)(i^^r—ivr)—i(wi—'w&)-zv’. 

(11.16) 

Therefore by adding 11,16 and 11.16 and making use of 11.11 
we find 


= —(i/;i^+?z>^) AiV, i.e. AN cc 
dt 

(11.17) 

which proves that AN approaches its equilibrium value {AN = 0) 
exponentially. 

Dipolar liquids 

It was shown in § 8 that when a disordered solid melts there 
should be no appreciable change in the static dielectric constant 
because in both the liquid and solid phases the average energy 
of a dipole is the same for all equilibrium directions. The fact 
that—^in contrast to liquids—^there is only a discrete number of 
equilibrium directions in solids has no influence on the static 
dielectric constant. 
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There is, however, an essential difference in the dynamic 
properties of solids and liquids. In solids in view of the interaotion 
of a molecule with its neighbours a dipole has a number of equi¬ 
librium directions. They are separated by potential barriers 
over which the dipole must pass in turning from one such direc¬ 
tion to another. In liquids the average distance between neigh¬ 
bours, and hence the interaotion, is about the same as in solids. 
Therefore, if the positions of all molecules but one could suddenly 
be fixed, this selected molecule would behave similarly to a 

A 

dipolar molecule in a solid, i.e. it would probably possess a 
number of equilibrium directions separated by potential barriers. 
It is an essential property of a liquid, however, that its mole¬ 
cules have no fixed positions. Thus, if we imagine one molecule 
to be turned out of its momentary equilibrium directions, its 
neighbours ■will tend "to rearrange themselves in such a way as to 
make this new direction an equilibrium position. So we may 
have •two conceptions of the way in which a dipole in a liquid 
alters its direction. It may jump into a different direction in a 
similar way as in solids. This requires—^at least for the duration 
of this jump—that the arrangement of its neighbours should 
remain rmaltered. The second possibility is that such jumps may 
occur very rarely and that a dipole may alter its direction only 
in conjimction with a rearrangement of the positions of i'ts 
neighbours. In fact a dipole might be considered as fixed fairly 
rigidly rda'tive -to its neighbours. The rotation of a dipole would 
then affect the motion of molecules at some distance from it. 

The average motion of these neighbouring molecules might be 
described by replacing them by a continuous medium with the 
properties of a macroscopic viscous fluid. This possibility leads 
to the model used by Debye [Dj 8] in which a dipolar molecule is 
considered to be a sphere of radius a moving in a continuous 
■viscous fluid with ■viscosity r) and obeying the macroscopic 
equations of flow. The fluid is considered to adhere to the sur¬ 
face of the molecule. On these assumptions the frictional 
constant $ of the sphere is given by Stokes’s law. 




( 11 . 18 ) 
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This means that when an eleotrio field f, making an angle 6 with 
the direction of the dipole {jl, is applied, 



—.ja/sinfl. 


(11.19) 


if no other forces act on the molecule. In this equation it is also 
assumed that inertia effects of the dipole—^which macros copi- 
calLy would lead to an additional term Id^djdi^ {I = moment of 
inertia) on the left-hand side—can be omitted. This omission is 
justified in so far as this form of the term is concerned since it 
would hardly be expected that such effects could be described 
adequately by this macroscopic expression. 

The macroscopic equation of flow 11.19 can yield correct 
values for 6 only if account is taken of thermal fluctuations due 
to the interchange of energy between the molecule and its 
surroundings. For the above model the energy of a molecule is 
the kinetic energy of rotation of the sphere and has an average 
value of order kT. This rotation—corresponding to Brownian 
motion in translation—changes its direction and magnitude very 
frequently due to collisions with the molecules of the surround¬ 
ing liquid. However, if no external forces act on the dipole its 
average displacement from a given direction vanishes because 
displacements in all directions have equal probability. The 
mean square of these displacements, on the other hand, will 
increase steadily with time. 

[Equation 11.19 may, therefore, be considei’ed to hold for the 
average value of 6 (for a single dipole) when a field is applied. 
The actual value of 6 may, however, only be expected to be near 
its average value as long as the mean square displacement due 
to thermal fluctuations remains a small quantity. For a macro¬ 
scopic dipole this should practically always be so. According 
to 11.19 such a dipole will, therefore, grad-ually approach the 
direction d = 0 parallel to the field. In equilibrium it will remain 
in this direction apart from small fluctuations. 

It is to be noted that the value d used in 11.19, and assumed 
above to be the time average of the 0-values of a single 
dipole, may also be considered as the average of such angular 
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dispiaoemeuts of a large number of (non-interacting) dipoles all 
having the same ^-vaJue within a small range. 

In contrast to a macroscopic dipole the fluctuations of a 
molecular dipole are very large. This follows at once from con¬ 
sideration of the equilibrium distribution of an assembly of 
dipoles which, according to 6.9 (replacing E by/), is proportional 
to exp(/i/cos 6jhT). Since ixf kT for all practicable fields 
(cf. § 4), the field has only a slight influence on the distribution 
function. In other words, a dipole has only a slightly higher 
probability of being near 0 = 0 than in the opposite direction. 

In order to derive now the macroscopic dielectric properties 
of dilute solutions of dipolar molecules in non-polar liquids, 
consider an assembly of noii-interacting dipoles and let 

t)BmddS 

be the number of dipoles per unit volume in a range dd near 6 
at time t. The function iV'(0, /) in general varies with time because 
individual dipoles continuously change their directions. These 
changes are either due to thermal fluctuations or to the action 
of the field. The latter, according to 11.19, causes all dipoles 
to be displaced by an amotmt S0., 

S0 = ^ S< = —“.sin0Si, (11.20) 

at § 

in a short time interval 8f. Therefore within ht seconds a number 
N{d, i)sin 6 Sd will pass through the surface of a cone with angle 0. 
Owing to the action of the field the rate of change of the number 
of dipoles within an interval dO near 0 is, therefore, using 11.20, 
given by 

*)sin0 80} = t^de-^{N(e, «)sin30}. 

( 11 . 21 ) 

This equation follows because a number N{d,t)shid 86 will 
leave the range dd within Si seconds, whereas the number of 

dipoles entering it is given by the same expression replacing 6 
by 6 — dd. 
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£Vom 11.21 it follows that the rate of change of the function 
Ni0, t) itself due to the action of the field is given by 


In the case of large fluctuations with which we are dealing at 
present this equation is more suitable than 1J..1& to express the 
effect of the field on dipoles. 

^Equation 11.22 can be used to find the rate of change of the 
dipole moment Mf in the field direction, of the assembly of dipoles 
under the influence of the field. Assuming the dipoles to be rigid. 


Mf = /X J cos 0 f)sui 5 <20 (11.23) 

o 

per unit volume. Therefore using 11.22 the rate of change of 
Mf due to the action of the field is given by 





o 


sin 0 (20 



or integrating in parts. 


^ J <20 = 

* (11.24) 

where Nq is the total number of dipoles per unit volume and 
wn®0 is the average value of sin®0. As discussed above, the field 
influences the angular distribution of dipoles only very slightly. 
Therefore the value of in equilibrium in the absence of a 

field can be inserted in 11.24, i.e. 



sin*0 sin 0 <20 



sin 0 (20 



2 

3 ^ ■ 



(11.26) 


Hence 


(11.26) 
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The total rate of ohange of Mf contains in addition to 11.26 a 
term due to thermal motion. This term will tend to restore the 
equilibrium distribution in the absence of a field, in which case 
Mf = 0. On tihe assumptions which we have made, namely (i) 
that there is no interaction between dipoles, (ii) that in a short 
time interval 8^ the 0-value of a dipole is altered only very slightly, 
this second term must be a linear function of the deviation from 
equilibrium for / = 0, i.e. it must be proportional to — Mf. 
Therefore introducing 1/t as proportionality factor, the total 
rate of change of Mf is given, using 11.26, by 


dMf_ Mf,2iL^fN^ 
- ~+ 3 —- 


(11.27) 


In view of the linear relationship between Mf and dMfjdt 
equation 11.27 leads, of course, to exponential approach to 
equilibrium and hence, according to § 10, to the Debye equations. 
This result in fact should always be expected when the following 
conditions hold: 

(а) Absence of interaction between dipoles. 

(б) Only one process leading to equilibrium (e.g. either transi¬ 
tion over a potential barrier, or fHctional rotation). 

(c) All dipoles can be considered as in equivalent positions, 
i.e. on an average they all behave in a similar way. 

The value of r can be derived from 11.27 by making use of the 
equHibrium value of in the presence of a field. Thus, using 


6.16, 



(11.28) 


This expression must be equal to the equilibrium value for 
Mf resulting from 11.27 with dMfjdt — 0, i.e. 




(11.29) 


Hence comparing 11,28 and 11.29 and using ii.is, 

_ i _ 

2kT ~ ~TT~' 


(11.30) 


This important formula duo to Debye [D2] is often used 
to discuss the relaxation time of dipolar liquids. It seems 



ni.fll MODELS FOB. THE DEBYE EQUATIONS 89 

appropriate, therefore, to recapitiilate the asBumptions imder 
which it shoiild hold. First of all the above-mentioned assump¬ 
tions (a), ( 6 ), and (c) require 

(a) dilute solutions of dipolar molecules in a non-polar liquid; 

(jS) axially S3anmetric molectdes; 

(y) isotropy of the liquid—even on an atomic scale in the 
time average over an interval small compared with t. 

Here (a) follows from (a), and (fi) and (y) are both connected 
with {b) and (c). All three assumptions are necessary for the 
Debye equations to hold; but they are not sufficient to obtain 
the value 11.30 for the relaxation time. This relation, as has been 
discussed above, is based on the further assumption that the 
dipolar molecule is fixed fairly rigidly relative to its neighbours 
so that large jumps of the dipole direction are unlikely. Now 
the empirical temperatme-dependence of the viscosity, 

170 c (11.31) 

suggests that jumps over a potential barrier of height JSL are 
carried out by the molecules of the liquid in processes connected 
with viscous flow. Therefore if is the height of the potential 
barrier related to jumps of the dipolar molecule, the condition 
that these jumps happen only very rarely suggests that 

(S) H > 

This condition, however, is correct only if the coefficients 
A (cf. 11 . 3 ) have the same order of magnitude for both types of 
transition. In general it should be expected \_F12, /Si5J that low 
viscosity liquids have small values of and they are more likely, 
therefore, to satisfy (S) than are high viscosity liquids. 

Formula 11.30 can thus be used to find the temperature 
dependence of t which is essentially the same as of 17 , i.e. 

Tocc^*^, (11.32) 

because the other terms vary only slowly with T. It follows that 
this temperature-dependence is independent of the nature of 
the dissolved molecules and is a function of the viscosity of the 
solvent only. To obtain the absolute value of t- as well would 
require knowledge of the effective radius a which is the radius of 
a solid sphere having the same frictional constant as the dipolar 
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moleeule. No theoTetical investigations on the correlation of 
this radius to the molecular radius normally used have yet been 
carried out. One should expect, however, that a depends on 
both solvent and solute and, therefore, caimot be assumed to be 
a molecular constant. Thus at the present stage of development 
equation 11.30 cannot be used to obtain absolute values for the 
relaxation rime. 


12. Generalizations 

The derivation of Debye *s relation 11,30—connecting the re¬ 
laxation time of a dipole molecule with the viscosity of the 
liquid in which it is dissolved—^required the assumption that the 
dipolar molecule is bound so strongly to the surrounding mole¬ 
cules that large jumps of the dipole direction are very unlikely. 
This may be true for a number of oases, but others may exist in 
which the opposite is more likely. A dipolar molecule will then 
make many jumps over the potential barrier separating it jErom 
another dipole direction during the rime required for an appreci¬ 
able change in direction by viscous flow. Clearly this holds for 
solids where flow may be considered as entirely absent; but some¬ 
times it may also be expected in liquids, and in particular in 
amorphous substances for which the viscosity is so high that 
flow is practically negligible. In liquids it might also happen that 
the process which prevails is different for different kinds of 
dissolved molecules. A further possibility mentioned by Sohalla- 
mach is the coexistence of both types of transitions. 

As an example of the type of substance exhibiting the second 
typo of behaviour mentioned above we may consider a dilute 
solution of dipolar molecules in a liquid or in an amorphous solid. 
In this case, in contrast to § 11, we assume that the dominant 
process for changing the dipole direction is that involving many 
large jumps. A single dipole will then behave similarly to a 
dipole in a crystalline solid j as in 11.3, the probability for a Jump 
over the potential barrier (height H) wiU be proportional to 
SfJcT)* In contrast to conditions in crystalline solids, 
however, the arrangement of the nearest neighbours is not 
exactly the same for all dipoles, and hence the heights H of the 
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potential barriers will also differ, thus producing different 
values of the transition probabilities. The dipolar mole¬ 
cules in such a substance may then be classified according 
to the height H of their respective potential barriers. If the 
substance has been polarized, and the external field has been 
removed, the contribution of the dipoles in a small range of 
energies near H will decay exponentially—as shown in § 11— 
with a relaxation time t related to H by equation 11.3. There¬ 
fore, instead of H the individual relaxation time r may be used 
to classify the molecule. Let y^r) dr be the contribution to the 
static dielectric constant of the group of dipoles having individual 
relaxation times in a range dr near r. Since interaction between 
dipoles can be considered to be absent (dilute solution), the 
contributions of the various groups superpose linearly. Their 
total contribution to the static dielectric constant is, therefore, 
given by * 

•^8—€« = J y(r)dr. (12.1) 

0 

The function y(r) which describes the distribution of relaxation 
times will be called the distribution function. 

To obtain the i complex dielectric constant we first consider 
the decay function «(<) (of. §§ 2 and 10). The dipoles with relaxa¬ 
tion times in a range dr near r make a contribution to (x{i) which 
is proportional to exp(— tjr) and to y{r)dr/r, which corresponds 
to the coefficient of the exponential term in equation 10.8. 
Therefore the total contribution of all the dipoles is given by 


OO 


«(0 = J* e-^''y{r)^. 


( 12 . 2 ) 


Using relation 9.1, the complex dielectric constant is now 
obtained from 12.2, 

CD OO OD —I 

)—«<*, = J ocix)e^*^^dx j J ^c-^l-'yir)^ 


OO 


= J ^y{r) J dx (12.3) 


00 


0 



92 


DYNAMIC PROPERTIES 


III, § 12 


The last step is due to an interchange of x and t integration. 
The a;-integral is then identical with the integral in 10.14. There- 
fore 

= f (12.4) 

J 1 —^orr 
0 

or separating real and imaginary parts according to 2.8 similarly 
to 10.16 and 10.16, 





y{T)dT 

l+a>M’ 


( 12 . 6 ) 


€ 


2 



y{T)wT dr 
l-ha>®T* ‘ 


( 12 . 6 ) 


!For a substance of the type described in this section these equa¬ 
tions will replace the Debye equations (10.15, 10.16, 10.17), and 
the latter will not be satisfied. A clearer insight into the meaning 
of equations 12.5, 12.6 will be obtained by a discussion of the way 
in which their solutions deviate from those of the Debye formulae. 
These deviations are best considered in connexion with the 
shape of the power loss—^frequency curves € 2 ( 0 )). For a detailed 
discussion, knowledge of the distribution function ^(t) is, of 
course, essential. As a preliminary step, however, we may note 
that y{T) is always positive and that e 2 (a>) therefore consists of a 
superposition of Debye curves cuT/(l+a>M) (of. Fig. 13) with 
diiferent positions of their respective maxima. The resulting 
curve € 2 (a>), therefore—supposing it has a single maximum— 
has a larger half-width than that of a single Debye curve whose 
maximum coincides with that of 62 (u>). As a simple example 
[FT] consider a model in which each molecule has two equilibrium 
positions with opposite dipole directions, and with equal energy 
in the ground level, as in the case of disordered solids; but in 
contrast to the latter the potential barrier between the two 
positions has a different height for each molecule. It will be 
assumed that the heights II of the potential barriers are equally 
distributed over a range between and i.e. 


H - 



(12.7) 
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Thus if Nq is the total number of dipoles per unit volume, 

( 12 . 8 ) 

is the fraction with jy-values in a range dv near Hq-\-v. 

For dilute solutions interaction between dipoles can be 
neglected. As shown in §§ 4 and 6, the contribution of a dipolar 
molecule to €g is then independent of H, and therefore is the same 
for all molecules, i.e. (e,—per molectile. The individual 
relaxation time t, however, depends on H according to 11.3 and 
11.8. Hence, using 12.7 and considering A as constant. 


2a>„ 


(12.9) 


The individual relaxation times t therefore cover the range 

tq ^ T ^ Tj where (12.10) 

To determine the distribution function ^(t) we note that 
^(t) = 0 outside the range 12.10. Now consider r as a function 
of V so that from 12.1, using 12.9, 






( 12 . 11 ) 


This means that 2 /(t)t dvjkT is the contribution of the molecules 
in the range dv to the static dielectric constant. On the other 
hand, we have just seen that this contribution per molecule is 
(e,— €go)INQ independent of v. Therefore, since 12.8 represents 
the number of molecules in dv. 


y(T)T dv _ eg—€„ dv 

- .-^ 


kT 


No 


V, 


( 12 . 12 ) 


or 


kT 1 


y(T) == (eg—O—-» if To ^ T < Ti = 


Vn T 


(12.13) 


2/(t) == 0, 


if T <C, tq and t 


The above distribution function—corresponding to an equal 
distribution of potential barriers over a range Vq- -depends on 
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temperature. In partiotilar its relative width (tj— de¬ 
creases with temperature since 


—Zo *= (12.14) 

'To 

The dielectric constants and can now be obtained by 


1 



I/too 1/10 t 10 100 

Co/COfn 


iFxci. 14* Depeoadenoe of dieleotrio loss €j|(ct>) on angular frequency cu acoord- 
ing to equation 12.19 for the three values 1, 5, and 10 of the parameter 
They correspond to a range of heights of potential barriers of width 
» ^logrx/T0. ci(u>)/€s(a>f„,) is plotted against co/c<i^ on a logaiithmio scale; 
xepresenta the v^ue of io where reach^ its maximum. 


inserting y(r) from 12.13 into 12.6 and 12 . 6 . All u 
be carried out in an elementary manner leading 


ei(a}) 


'CO 


— ®ot>)f !• ““ 




2v. 


log 


l-fa>*7^ / 


(12.16) 


Tel' 

ea(tu)= («,—€«) —(tan-i(toToe»W*' 2 ’)—tan-Wo). (12.16) 

^0 

These equations replace the X>ehye formulae 10.16 and 10.16 
for omr present model. In this case ei(o) and e 2 (a>), considered 
as functions of fiequency, depend on two parameters—^a relaxa¬ 
tion time T 0 , and the factor v^fTcT determining the width of the 
range of relaxation tunes by 12.14. The Debye formulae are 
obtained for v^jhT = 0 ; they contain, therefore, one para¬ 
meter only. The shape of cg can be seen from Kg. 14, where 
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is plotted against oijco^, being the jBreqnency at 
which ^2 hs.s its maximum. 

It can be seen that •with increasing values of the second para¬ 
meter vJkT the curves become more and more flattened out. 
The position of the maximum is obtained from 12.16, using 


day 



if a> = a>^. 


This leads to 


CO 


m 


= — _ 


0 




(12.17) 



Fia. 15. Dependenoe of the Tnayi-miim value of the loss, ei(€t>f^), on the 
width Vq of the range of heights of potential barriers. 


which for == 0, of course, becomes identical with the corre¬ 
sponding value 10.22 for the Debye formulae. Inserting 12.17 
into 12.16 yields the maximum value of cg, 


h'V 

[tan-^e»««/*^—tan-ie-4««/*2’]. 


(12.18) 


Again for Vq = 0 this quantity becomes equal to 10.23, but it 
decreases with increasing v^kT, as shown in Pig. 16. 

With the help of 12.14, 12.17, and 12.18 equation 12.16 can 
be written in the form 


g2(“>) 

€2(a>wi) 



(12.19) 


which shows clearly the dependenoe of €^{oj) on the parameter 
tq/ti < 1. For to/ti — 1, equation 12.19 is identical with the 
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corre^nding equation of § 10, and the smaller is, that is 

the larger w,, the larger is the deviation from a Debye curve. 

^ above equations are based on the assumption of ecjual 
distribution of the height of potential barriers over a ran^e 
leading to a range of relaxation times given by equation 12.14 
and to the distribution function 12.13. This may seem to be 
inther a specialized assumption. It can be shown quite readily, 
owever, that the behaviour of in the main absorption ranee 

IS altered but little, if the distribution function 12.13 is replaced 
by any other smooth function which is large between r„ and t, 
and small outside this range. Quite generaUy, therefore, in the 

^ maximum of ea(a>) this function can he 

consider^ to be determined by two parameters, namely the 

a-t which the maximum occurs, and the 

times r,-T„. It should bo 

1 ^ _ . . * ® * solute value of e,(eu) is then determined 

^In ***® phenomenologiool relation 2. 1 8. 

re^u in the main absorption 

o^tS31 *“«ti°'> Vi-) for r-valnos 

relaiatiou timef. Bhr tholS;^^'* 

only small contributions numerous, make 

Sr“°‘r “ 

small value which changes verv sl3 ’ remains at a very 

'^is may be explained^a svmienZlir ^ frequency [02], 
times is used. Garten ron haa ^ range of relaxation 

absorption may be due to the fA ^’^S^e^ed that this residual 
eqnflibrium portions of tht 

w^ofdepthmmarC'^L of suoh a potential 

porbonal to i • -Big. 16) is assumed to be pro- 




pro 
( 12 . 20 ) 
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This is a plausible assumption, though no proof for it has been 
given. W^e shall proceed, to show that such an assumption does, 
in fact, lead to a value of « 2 (to) which is substantially independent 

of frequency for frequencies well outside of the main absorption 
range. 

Liet H be the height of the potential barrier measured from a 



Eio. 16. Potential well of depth to due to thermal fluctuations. 

normal (permanent) position and consider a molecule near which 
a temporary position of depth w has been formed. According to 
11.3 the time it spends in the temporary or in the permanent 
position is given by 

T = and tq = (12.21) 

respectively on the further assumption that the temporary 
position exists for a time which is long compared with t+to. 
The coefficients B and B^ are considered to be independent of w. 
Assuming t tq, the relative probability of finding this mole¬ 
cule in the temporary position is given by cri t/tq. The 

average number of molecules in such positions, using 12.20 and 
12.21, is thus proportional to 

—e-^l^^dw 2 = — — oc —. (12.22) 

Tq Tq T fltr Vq T T 

£!ach molecule makes a contribution to proportional to 

<*>T/(l-l-ajV*). Also since w > 0 the smallest value that t can 
take is T = J5. For large values near t == tq the developments 
become invalid, and tq will be used as the upper limit for t; this 

4980.11 u 
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should not be an essential limitation, however, beoause near 
T = To the main contribution to ea(a)) is due to the transitions 
between permanent positions which are much more numerous. 
Thus 

T« 

60 ( 0 )) oc f ■ L — = tan~^a)To—tan“^jBai (12.23) 

s 

for frequencies higher than the main absorption range, cuto ^ 1. 
If, in particular, 

1/to 1 / B , (12.24) 

then ei(a}) is practically independent of < 0 . Here IJS will 
probably be a very high frequency. 

The preceding discussions were all based on the assumption 
that interaction between dipoles can be neglected. It seems 
likely that the influence of interaction will always tend to 
broaden the Debye absorption curve, because the energy of a 
moleoide in an equilibrium position may have a whole range of 
values depending on the positions of its neighbours. One would 
thus expect the heights of the potential barriers to vary accord¬ 
ingly. In the present case, however, the time which a molecule 
spends in a given equilibrium position is of the same order as 
the time during which the height of the potential barrier remains 
constant. If an attempt is made to include interaction under 
these conditions, great mathematical di£Q.culties are encountered 
and no solution for this case has yet been obtained. 

13. Resonance absorption 

In § 9 it was pointed out that two types of power loss should 
be expected to occur in dielectrics: (i) The loss due to displace¬ 
ment of charges bound elastically to an equilibrium position. 
Such charges have a proper frequency of oscillation, say a>o/2ir, 
and the power loss (and hence e^) will be expected to have 
a maximum, near this frequency (the resonance frequency), 
(ii) The loss due to transitions of charges or dipoles between 
equilibrium positions separated by a potential barrier. Such 
transitions axe described by a relaxation time t, and in this case 
the power loss will have a maximum near a frequency 1 as 
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discussed in §§ 10—12, This frequency is usually strongly depen¬ 
dent on temperature, in contrast to the resonance frequency in 
case (i). The cliief aim of the present section is to derive expres¬ 
sions for the complex dielectric constant when there is resonance 
absorption (case (i) above). 

It is our intention to consider only the simplest case, and it 
will be seen presently that the required formula can be derived 
from very general considerations without choosing a particular 
model. We shall employ the same method that has been used in 
§ 10 to derive the Debye equations. The polarization of the sub¬ 
stances considered there was assumed to approach equilibrium 
exponentially with time. This cannot be expected to hold in the 
present case of elastic displacement. One would rather expect 
damped oscillations with frequency cuJ'Itt about the equilibrium 
polarization, as has ali’eady been discussed in § 9. Therefore 
instead of 10.1 we shall now assume a decay function 

a(0 = (13.1) 

and then, making use of equation 9.1, derive an expression for 
the complex dielectric constant e. Equation 13.1 contains two 
constants y and ifj which will be determined below. 

The reader is reminded that a{f>) is a macroscopic quantity 
referring to the behaviour of the polarization of the whole sub¬ 
stance. It would be wrong to describe the amplitude of the 
oscillations of a single molecule (in the absence of a field) by a 
function similar to 13.1, because in equilibrium this would lead 
to vanishing amplitude and thus would not take account of the 
thermal motion. 

By inserting 13.1 into 9.1 we now find for the complex 
dielectric constant e. 


CO 


e — € 


oo 


= y J c“®/'^coa(a»o^+0)<’*‘'''d:t’ 


^ yzl _ 

2 \l—i 






I ] 

^ 1-f-l(coo— co)t/ 


VT 

= -^cos 


J .y (13.2) 

-;(co„-f-a>)r^ 1+/(a*o-^)T/ 
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As the next step, the two constants y and ^ will be determined 
by considering the ILtniting cases of very low and very high 
frequencies. Low frequencies, cu oig, are consideired to be 
outside the main range of absorption, i.e. in a region in which 
^2 is very small and, therefore, the dielectric constant is prac¬ 
tically real. In this range e will approach a constant value, say 

€ = eooH-Ac, oi a»o, (13.3) 

where Ae is a real quantity. This value of the dielectric constant 
will hold down to frequencies at which another absorption range 
starts. In the absence of further absorption regions 13.3 repre¬ 
sents the static dielectric constant. The high-frequency dielectric 
constant Co,, on the other hand, is the value of e at the high- 
frequency side of the absorption peak, again in a region in which 
^2 is very small. 

Comparing 13.3 with equation 13.2 and neglecting eu compared 
with ojg, we find 



cos ^— wq t sin ift 
l-l-<ogT® 


= yrGOBifs 


1—euQrtan^ 

l-t-a»§T* 


(13.4) 


as a first condition for the determination of y and iji. 

The second condition is of a more subtle nature and concerns 

the behaviour of e at angular frequencies oj ^ wq. At these 

frequencies coq can be neglected compared with a>, so that 13.2 

becomes , 

cos W o r V 

e—Coo = y'T- -r^, O) > OJO- (13.6) 

1 - tOiT 


On the other hand, for a> ^ coq the influence of the restoring force 
should be negligible during one period 1/co. Therefore a be¬ 
haviour according to the Debye theory should be expected wdth 
a static dielectric constant given by 13.3, Thus by replacing 
€,—Cflo by Ac in 10.13, 

€“■“€00 . , CO ^n. (13.6) 

1- ^COT 

Combining 13.5 and 13.6 we find 

Ae = yrcost/i (13.7) 

as the second condition for the determination of y and ijj. 
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From 13.7 and 13.4 we now conclude that 

coqt = —tan 0 . (13.8) 

The insertion of this expression and 13.7 into 13.2 finally leads 
tof (cf. references V5 and F9) 

— i tOpT 1 + 

- 0 -|-a>)T *^ 14 - 4 ( 01 , 

or separating reaJ and imaginary parts according to 2 . 8 , 


'00 


lA i 1 — 
= ^Ael-- j— 

\1— 




(13.9) 


ei(<«) 


— €00 = 


l~h<»*o(<^H~u>o)T^ 1—a>o(<o— 


+ 


1 ~l~(<U-|-<Oo)®T® l-j-(6lJ 


and 




2—2 ~ 


0>T 




< 0 —a>o)Tf\ 

)■ 


(13.10) 

(13.11) 


These equations represent the dielectric constants in the 
simplest case of resonance absorption. In many practical cases 
Ae is a small quantity so that is approximately equal to 600 . 
In these oases the loss angle according to 2.6 is given by 


tan^ , ifA€<l. (13.12) 

®00 

If Ac is small it will be difficult to determine this quantity 
experimentally with reasonable accuracy by direct measure¬ 
ment. Following 2.18, however, the relation 

OO 

Ae = — (* C 2 («*>) = — r 62 d(]oga>) (13.13) 

77 J to 7T J 

o 

can be used instead. This equation holds for the present case, 
as can be chocked by inserting Cg from 13.11. 

We shall now proceed to discuss power loss (oc € 3 ) in more 
detail. At constant temperature as a function of a» has a maxi¬ 
mum when CO = co,„, 

= (13.14) 

T 

as can be found from the condition dejdcj = 0 . The maximum 
value of €3 is then given by (inserting 13.14 into 13.11) 

€2(to,n) = ^Ac.y/(l-f-eogT*) = iAeco,„T. 

I 

t Another derivation of 13.0 is given in the Appendix A 4. 


(13.16) 
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Consider now that the time of relaxation t depends on tem¬ 
perature and may he expected to decrease with increasing 
temperature. 

Then 

1/t and ea(<o^) Cii if coqt 1 , (13.16) 

i.e. at very high temperatures. But 

co^ C:i oJq and e 2 (^m) — ‘iAectiQT if <OqT ^ 1, (13.17) 

I 

i.e. at lower temperatures. Comparison with 10.23 shows that 
case 13.16 leads to behaviour similar to that for Debye absorp¬ 
tion. Without detailed consideration of a special model it is 
impoE»ible, however, to say whether this region can be reached 
at temperatures corresponding to the solid state. The low- 
temperature limiting case 13.17 leads to typical resonance 
absorption. In contrast to Debye absorption, the angular 
frequency cjq of maximum absorption is independent of tem¬ 
perature, but the resonance peak becomes narrower and higher 
the lower the temperature, i.e. the larger coq t. Figures 17 and 18 
give a comparison of Debye absorption and resonance absorp¬ 
tion for various temperatures. 

A number of models have been considered in greater detail, 
and equations 13.9—13.11 have been found to hold. These deriva¬ 
tions are much more complicated than the one given above, 
but they may be considered to be more rigid (of. Appendix A 4). 
Transitions between rotational levels of gases have been con¬ 
sidered by van Vleck and Weisskopf [F5] and by van Vleck [F4J. 
In gases the relaxation process is due to collisions between mole¬ 
cules. Therefore t will be expected to decrease with increasing 
pressure and with increasing temperature. At temperatures or 
pressures at which 1 the absorption peak according to 

13.16 is hear the frequency = 1/t. In this region, therefore, 
<o ^ < 0 g, so that cuq can be neglected. 

In the main absorption region, therefore, the Debye equations 
are fulfilled whenever cuqt 1 . At small pressures, however, 
<i>gT ^ 1 and the main absorption peak according to 13.17 will 
be foimd near the resonance frequency ojq. 

Biesonance absorption may also be of importance in solids 
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oontaining dipolar molecules with large moments of inertia 
\F9, H3, The resonance in this case is due to rotational 



Fig. 17- Frequency depezidenoe of dieleotrio loss for three different tempera¬ 
tures (sohematioally)- (a) Debye case ; (b) resonance case. 



its maximum on the relaxation time. JD refers to the Debye cause* J2 to the 

resonance case. 

oscillations of the moleovdes about their equilibrium positions. 
The behaviour of the relaxation time in this case, however, has 
not yet been investigated. 




CHAPTER IV 
APPLICATIONS 


The present chapter is designed to show by extunples how the 
general theory developed in previous chapters can be applied 
in the discussion of the properties of dielectric materials. It is 
the intention to cover some important types of materials (e.g. 
non-polar solids, ionic crystals, dipolar liquids, etc.); this should, 
however, not be considered as an attempt to discuss systemati¬ 
cally the properties of dielectric materials, but rather as a 
selection of t 3 qpioal examples. 

14. Structure and dielectric properties 

AUma 

In the present section we shall commence by discussing the 
structure and the dielectric properties of atoms and molecules 
leading up to a classification of dielectric materials according 
to their structure. An atom, as the reader will know, consists 
of a positively charged nucleus surrounded by electrons which 
just compensate its charge and whose combined mass is very 
small compared with that of the nucleus. By adding or removing 
one or several electrons one obtains negative or positive ions. 
In a stationary state the electric dipole moment of an atom (or 
of an ion relative to its nucleus) vanishes. 

A very general quantum mechanical theorem says that under 
the influence of external (static or alternating) electric fields the 
electrons of an atom behave like an assembly of classical har¬ 
monic oscillators whose firequenoies and other properties can be 
specified, but will not be required by us. Usually the frequencies 
are equal to or higher than those of visible light, and, therefore, 
are much higher than the frequencies of electric fields contem- 
plated in the present book. Under the influence of an external 
electric field f an electric dipole will, therefore, be induced 
which has all the characteristics of a dipole produced by elastic 
displacement of electric charges (electrons) as discussed in case 
(i) of §§ 4 and 9. In the ranges of field-strengths and frequencies 
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to be considered here, the strength of the induced dipole is 
proportional to the field-strength /, but independent of the 
frequency, and it shows no phase-shift relative to the field. These 
properties of an atom can be described by a constant polariza¬ 
bility Ota, ^ (14.1) 

The resulting polarization will be called optical polarization, 
to indicate that the characteristic frequencies are usually in the 
optical region. Often it is also described as electronic polariza¬ 
tion because it is due to the displacement of the electrons; thia 
description may, however, be misleading as will be seen presently. 


MolecvZes 

Consider now a diatomic molecule A B composed of two atoms 
A and B. In view of the interaction between the atoms—Pleading 
to the chemical bond—^the distribution of electrons in the mole¬ 
cule is different from a mere superposition of the electron distri¬ 
butions of the free atoms A and B. The distribution should, 
however, be axially symmetrical around the line joining the 
nuclei of the two atoms. It must, therefore, be expected that a 
diatomic molecule will possess a dipole moment in the direction 
A—B. An exception occurs in the case of two equal atoms, 
A = B, where the dipole moment must vanish for reasons of 
symmetry. Thus, for instance, the molecules HCl or CO have a 
dipole moment in contrast to Hj, O 2 , or Clg. The magnitude of 
the dipole moment gives valuable information about the distri¬ 
bution of electrons. Thus the largo <lipolc moment of HCl 
(fi 1 X 10~i®e.K.u. — 1 Debye unit) compared with that of 
CO (fi. ~ ()• 1 Debye units) gives some justification for considering 
the HCl molecule (but not the CO molecule) as composed of a 
positive and a negative ion, H+-t-Cl~. 

Under the influence of an external field f an average dipole 
moment m is induced in a molecule which, according to §§ 4 and 
9, can be considered as composed of contributions of two types; 
(i) those due to elastic displacement of charges, and (ii) those due 
to a change in the average orientation of the permanent dipole 
of the molecule. The former again are composed of the contri¬ 
butions of the various normal modes. In view of the great 
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differenoe in the mass of electrons and nuclei there is one group 
of normal modes ■which is practically entirely due to displace¬ 
ment of electrons relative to the nuclei. Their proper firequencies 
are usually in 'the op'tical region, visible or ■ultra-violet. The 
proper feequenoies of the normal modes connected with 
oscillations of 'the nuclei, on the other hand, are usually in the 
in&a-red region. They must con'tain displacements of electrons 
relative ■to ■their respective nuclei as well, because of the variation 
of the in'teraction be'tween electrons and nuclei with the distance 
between the latter. Thus 

where the jBrst two terms are due to elastic displacements with 
proper frequencies in the op'tical and in the infra-red region, and 
is the moment due ■to dipole orientation. The quanti'ties 
and nifr are often denoted as electronic and atomic polarization, 
respectively, which may be considered as somewhat misleading 
beoaiase the lat'ter also contains terms due to electronic dis¬ 
placements. Since the frequencies contemplated in this book 
are much smaller ■than all proper frequencies, the relations 

= OLi^f ( 14 :. 3 ) 

should hold with constant values of the polarizabilities and 
cq,., independent of the frequency of the field or of temperature. 

Polyatomic molecules behave in a similar way 'to diatomic 
molecules in that ■the average dipole moment induced by an 
external field can be considered as a superposition of three terms 
according to 14.2: 'two ■terms representing the elastic displace¬ 
ments ■with proper frequencies in the optical and infra-red 
regions respectively, and a dipolar contribution. The former 
can be obtained from the optical and the in&a-red polarizabilities 
do and (Xf,. according 'to 14.3; the value of the latter depends on 
the dipole moment of the molecule (e.g. § 6). The determination 
of the dipole moment is of great in'terest in investigations on the 
structure of molecules. Por instance, the fact that CO 2 does not 
possess a permanent dipole moment leads to the conclusion that 
the three atoms must be arranged in a straight line ■with the car¬ 
bon atom half-way be'tween the two oxygen atoms. In contrast 
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to this the H 2 O molecule has a permanent dipole moment, which 
means that HgO is a triangular molecule. 

A great number of interesting applications to steireoohemistry 
have been made, but a detailed discussion of this subject is 
beyond the scope of the present book. The reader will find 
examples and reference to the literature in the fairly recent book 
by liO Fevre \L2'\ or in the article by Sutton The earlier 

developments are to be found in the hooks by Debye [D2] and 
by Smyth [/Sf5]. 

Two qualitative rules concerning large molecules should be 
mentioned. Such molecules often contain a number of dipolar 
molecular groups such as the hydroxyl group O—, or the ketone 
group C=0. To a fair approximation these groups contribute 
the same dipole moment in whatever molecule they occur. The 
total moment of such a molecule is then the vector sum of the 
moments of all groups present in the molecule. Clearly this rule 
cannot hold exactly for several reasons; one of these is the in¬ 
fluence of the interaction between the various groups on their 
individual moments. 

The second rule refers to the optical polarizability of large 
molecules. As mentioned above, optical polarization is due to 
displacement of electrons on the assumption that all nuclei are 
kept in fixed positions. Clearly the largest contributions will 
arise from the electrons with the smallest binding energies, i.e. 
from valency electrons. The behaviour of these electrons is 
different in a compound than it is in the atomic state of the atoms 
contained in the compound. Thus the optical polarizability of 
CO 2 is different from the sum of the polarizabilities of a carbon 
and two oxygen atoms. On the other hand, the distribution of 
electrons in certain molecular groups is to a fair approximation 
independent of the type of molecule in which they occur, as in 
the above case of the dipole moment of such groups. In applying 
this rule it should be remembered that frequently the polariza¬ 
bility of such a group depends on direction; for the ketone group 
(!=0, for instance, the polarizability in the direction of the 
line connecting the two atoms would bo different from that in a 
direction perpendicular to it. 
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Classification of dielectrics 

The division of polarization into three types which was intm- 
duoed above is of a very general nature. It leads at once to the 
following division of dielectric materials into three classes: 

(i) Non-polar substances showing optical polarization only. 

(ii) Polar substances having optical as well as infra-red 
polarization. 

(iii) Dipolar substances which in addition show also polariza¬ 
tion due to dipolar orientation. , 

In the first class of materials, non-polar ones, an electrical 
field produces elastic displacement of electrons only. This class 
contains all dieleotrios consisting of a single type of atom, 
whether they form gases, liquids, or solids. Examples are 
diamond, oxygen (solid, liquid, or as vapour), the inert gases, 
and many others. The dielectrics of the second class, polar 
materials, are capable of infra-red polarization as well as optical 
polarization. Substances of this t 3 rpe may contain dipolar 
groups of atoms, but these groups must show only elastic dis¬ 
placement. If, on the other hand, there are several equilibrium 
positions for the dipole, then the substance belongs to the third 
class of materials. 

The second class contains, first of all, substances consisting 
of molecules whose total dipole moment vanishes, though they 
contain dipolar groups of atoms. Examples are CO 2 , paraffins 
CHg—(CHg)^—CH 3 (cf.§ 15), benzene CgHg, carbon tetrachloride 
CCI 4 , and many others in the solid, liquid, or vapour phase. In 
most of these substances the infra-red polarizability is only a 
small fraction of the optical polarizability. From a practical 
point of view, therefore, their behaviour is very «imi1a.r to that of 
non-polar substances. At sufficiently low temperatures many 
solids consisting of dipolar molecules fall into the same category 
because the dipoles in these solids freeze in, i.e. the thermal 
energy is msufficient to turn them in a recisonable time into 
other equihbrium positions. The most representative substances 
of the polar type are, however, ionic crystals which may show 
very large infra-red polarizability. Examples are rock-salt 
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NaCl, other alkali halide crystals, Ti02 crystals, a>n.d most other 
crystals of salts. In contrast to molecular lattices, which have a 
whole molecule at each lattice-point, ionic lattices contain one 
ion at each lattice-point. Thus rock-salt, for instance, forms a 
simple cubic lattice of Na+ and Cl- ions (cf. § 18). In contrast 
to other dielectrics, most salts on melting become (ionic) con¬ 
ductors. 

All materials consisting of dipolar molecules belong to the 
third class, except in the very low temperature range men¬ 
tioned above where the dipoles are frozen in. Often in solids 
this process starts just below a critical temperature at which the 
substance undergoes a phase transition. In nearly all of these 
materials the turning of a dipole into another equilibrium posi¬ 
tion is connected with a turning of the whole molecule (e.g. for 
ketones, § 17). Exceptions are ice and some other crystals where 
a turning of the direction of a dipole may be obtained by a 
transfer of a H+ ion from one equilibrium position to another. 

15. Non-polar substances 

The simplest type of dielectric substances show elastic dis¬ 
placement of electrons only, and in the classification of § 14 are 
designated as non-polar. For such substances the lowest fre¬ 
quency Vq at which appreciable absorption occurs is usually in 
the visible or int he ultra-violet region. For all frequencies v 
which are less than by a sufficient amount, the dielectric con¬ 
stant should be independent of frequency. Thus for v vq the 
dielectric constant e should be equal to the static dielectric 
constant and should satisfy the Maxw'ell relation 6 = n*. That is, 

(16.1) 

should hold between the static dielectric constant and the refrac¬ 
tive index n at fioquencies v vq. Whitehead and Hackett [ W6] 
have recently checked this relation on diamond and found that 
it holds within the range of accuracy of the measurements. They 
measured the dielectric constant at frequencies between 600 
and 3,000 cycles per second and obtained the value 6*68±0*03. 
The refractive index n was obtained by extrapolating measure¬ 
ments at various wave-lengths in the optical region to long 
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waves where w is independent of wave-length. The accuracy of 
tliis procedure can be seen from Fig. 19 and leads to the value 
71 * =s= 6*66. This agreement means that absorption of electric 
waves at fr^quenraes between 3,000 cycles per second and optical 

frequencies must be so weak that (cf. 2.18) ( 2 / 77 -) J e^dvfv <. 0*03 

if integrated over this range. Actually diamond absorbs at 
infra-red frequencies, but the reason for this absorption has not 



Fia. 19. Dependence of the refractive index n of diamond on tlie wave num¬ 
ber according to measurements collected by Whitehead and Hackett . 

yet been found. The above condition, by the way, need not 
imply that is small provided the absorption bands are suffi¬ 
ciently narrow. 

Gases of non-polar molecules fulfil fairly closely the conditions 
required for the validity of the Clausius—Mossotti formula, 
(of. §§ 6 and 8 and Appendix A 3). These conditions require 
(i) elastic displacement only, (ii) absence of non-dipolar (short- 
range) interaction between molecules, (iii) isotropy of the 
polarizability of a molecule, (iv) isotropy or cubic symmetry in 
the arrangement of the molecules. Condition (i) is fulfilled for all 
non-polar molecules; (ii) holds so long as the distance between 
molecules is sufficiently large; (iii) holds for spherical molecules 
only, and (iv) holds always in gases. Hence deviations may be 
expected for non-spherical molecules (condition (iii) violated) 
and at very high pressure (condition (ii) violated). All conditions 
are thus fulfilled for the rare gases. The main check of the 
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Clausms—Mossotti formula consists in measurement of the 
dielectric constant of a gas as a function of its density. Por, 
since the number of molecules per c.c. is obtained with the 
help of the Avogadro number (6’02x 10 ®*), 

•^0 = i 6*02 X 10®» (16.2) 

(d = density, W = molecular weight)- 
It follows with 6.34 and 6.32 that 




2*62 X 10®* 



(16.3) 


Often the quantity 

p = 2*62xl0®*a (16.4) 

is defined as molecular polarizability. Using 15.3 and the Max¬ 
well relation (of. 16.1), 


^ We—1 
^ d €+2 


(16.6) 


should be a constant independent of density, temperature, or 
frequency. For the rare gases the following values for p have 
been found: 

He Ne Ar Er X 

p = 0'6 1*0 4'2 0*3 10 0*0* 

The increase of polarizability with atomic weight is mainly due 
to the increase in the number of electrons per atom. Also, 
relation 16.1 has been found to hold very accurately. 

Agreement of equation 16.6 with experiment holds for many 
other non-polar substances over a very wide range of densities. 
Thus we notice with van Vleck [F3] that for the right-hand 
side of 15.6 has the magnitude 3*860 for the gas and an almost 
identical value 3*878 for the liquid despite the fact that the den¬ 
sities differ by a factor of over a thousand. Similar agreement is 
found in many other cases, and since the Maxwell relation holds, 
e in 16.5 can be replaced by w®. In nitrogen, for instance, it was 
found from measurements of the refractive index that the right- 
hand side of 15.5 varied by less than 1 per cent, in a range of 
pressures between 1 and 2,000 atmospheres. 
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In judging the signifioanoe of such a striking agreement it 
should he noted that in spite of the large variation in the magni¬ 
tude of € —1, the largest values of this quantity are usually 
smaller than unity. In this case, even with the assumption of 
complete absence of interaction between molecules, there is 
fair agreement with experiment. This assumption (of. 6.14, 
using 16.4 and 16.2) leads to 


e—1 



(16.6) 


In order to compare 16.6 and 16.6 we note that if we introduce 
the molecular volume v given by 


V 


W 
d ’ 


(16.7) 


equation 16.6 is equivalent to 


Thus 


e—1 = 


Sp/v 

1—PIv' 


(16.8) 




€ 2 


1. (16.9) 


The first term in this series is identical with 16.6. The higher 
terms which are characteristic of the Clausius-Mossotti formula 


are usually small. Thus in the above-mentioned caee of nitrogen 
p/v is only slightly larger than 0-1 at 2,000 atmospheres. Com¬ 
pared with equation 16.6 the aaxisius-Mossotti formula intro¬ 
duces a correction of only about 10 per cent. 

Kirkwood [KSj has investigated the influence of the anisotropy 

of the polarizability of molecules and foimd that 16.9 has to be 
replaced by 


' * ~ + +—). (16.10) 

moMnrmg the enieotropy of the polarize- 

y. The effect ofthis and other corrections on the magnitude 

of € are usuaUy small, since they affect the smaU second-order 
term only. 
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There are a great number of materials which according to the 
classification of § 14 are polar, but which, for all practical pur¬ 
poses behave like non-polar substances. They comprise mole¬ 
cules which contain a number of polar groups in such a way that 
their resultant dipole moment vanishes. The inira-red polariza¬ 
bility of these molecules is very ofben so small that its contri¬ 
bution to the dieleotiic constant is negligible. Important 
examples are the para£B.ns, and in the following we shall show 
why the dipole moment of a paraffin molecule vanishes. 

A molecule of a normal paraffin CHg—(CHa)„—CHa consists 
of a chain of CHg groups with a CH 3 group at each end. The bond 
angle C—C—O is very nearly equal to the tetrahedral angle 
d= 2oos“'(l/V3) 109®. In solids the chain forms a plane 

zigzag (e.g. § 16, Fig. 28), but in liquids and gases this is not the 
case in general. Here rotation around the C—C bonds as axes 
may occur, but it will bo assumed that this does not alter the 
bond angle. Fach CHg and CH 3 group has a dipole moment 
whose direction will be discussed on the assumption that all 
four bonds of a C-atom make angles of 2 oos“^( 1 /V3) with each 
other.f Thus, if a C-atom is considered to occupy the centre of a 
regular tetrahedron, the bonds are directed towards the vertices. 
Now, according to § 14, the C—C bonds should have no dipole 
moment in contrast to C—H bonds which will have a moment. 
The resultant moment of a CHg group is then always in the 
C—C—C plane (of. Fig. 20 ) bisecting the C—C—C angle, where¬ 
as the CHa moment has the direction of the adjoining C—C bond. 
To show this we consider a cube with a C-atom at its centre 
(Fig. 21 ). Then the fpur bond directions point towards four of 
its corners which are not neighbours. These corner points form 
the vertices of a regular tetrahedron. Thus, if /4 is the C—H 
moment, the CHg moment is 2 /Ltoo 8^/2 — 2/ti,/V3. This may be 
considered as composed of two dipoles of strength /x in the direc¬ 
tions of the C—C bonds (dotted arrows in Fig. 20 ). Furthermore, 
the moment of a CH 3 end-group is also equivalent to a dipole /u, 
in the direction of its adjoining C—C bond. For the resultant of 

t It should be noted that in Fig, 28 the diatanoe refers to CHg groups os 
found from X-ray inoosuromonts and not to tho C—C diatanoe. 

4980.11 T 
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any pair of H-atoms of the OHs group is a dipole of strength 
2/i/V3 in the direction perpendicular to a surface of the cube, 
as has just been shown (this assumes that C—H dipoles in a 
CHg and in a OH 3 group are equal). The resultant thus forms an 
angle oos“^l/V3 with the diagonal in which the residual dipole 
lies, and thus contributes (1/V3) (2^t/V3) = 2fi/3 to the OH 3 
dipole. Now three CHg groups can be selected from OHg, but 
then each OH dipole is counted twice. Thus the moment of the 



1^0. 20. jHe 21. 

represents the dipole of a OH, group seen to lie in the 
plaae of a parafiQn molecule. It can be decomposed, into the two 
dotted arrows in the C—O directions. 

C-atom of a parafiQn molecule with its four neighbourSs two O- 
^ "atoms. The former are indicated by rings, the latter by dots; the 

^pole of the group. The lines joining the central atoms 
j. indicate the directions of the latter rather than their 

ce, actually the C—distance is smaller than the C—C distance. 


Hg group is f (2/i/3) = It follows that the paraffin chain 
can be considered as composed of C—C rods with a dipole fj, 
on each end. These end dipoles point in opposite directions, so 

at t e dipole moment of each rod, and hence of the whole 
chaan, vanishes exactly. 

M an example we mention pentane CgH^ which at 30° C. 
u^er a jm^sure of 1 atmosphere has a dielectric constant 

^ refractive index is n = 1-36, i.e. w® = 1-86 which 

shows that the Maxwell relation 16.1 is fulfiUed—probably 
wi experimen^ error. Hence the contribution of the 
a-red polarization must be very small. The following table 
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of measurements by Danforth shows that the Clausius— 

Mossotti formula also holds to a good approximation. 


Pressure 

€ 

(€-l)/(e+2) 

l/« 

p = «(e—l)/{e-(-2) 

1 atm. 

1-82 

0-216 

0-613 

0*356 

12000 atm. 

2-33 

0-308 

0-907 

0*339 


Again it can be seen that (of. 16.9) pjv 1 so that the zero 
approximation of 15.9 accounts for more than two-thirds of 
6— 1. More accurate comparison with equation 16.10 shows that 
a value o- c=l —0-06 is required to account for the experimental 
data. This shows that the next approximation—i.e. the first one 
which is significant for the Clausius—Mossotti formula—deviates 
by only 6 per cent, from 16.9. 

For an example exhibiting a larger contribution of infra-red 
polarization we turn to the measurements of Michels and 
Hamers [Jf2] and Michels and Kleerekoper [MS] on COg. This 
molecule is rectilinear and has, therefore, no dipole moment. 
Dielectric measurements give the value 7-6 for t?(e—l)/(e-l-2); 
from refractive index measurements, on the other hand, the 
value 6-7 is found for v{n^ —1)/(%®-|-2). Hence about 10 per cent, 
of the static polarization is due to infra-red contributions in this 
case. 

Relatively large infra-red contributions are obtained in dipolar 
solids at low temperatures at which dipoles are frozen in, so that 
these substances behave in a similar manner to non-polar 
materials (temperature-independent dielectric constant in a 
certain range). As an example, the case of ketones will be dis¬ 
cussed in more detail in § 17. 

16. Dipolar substances 

Oases and dilute solutions 

Provided the pressure is not too high, the static dielectric 
constant e^ofa dipolar gas satisfies the relation Cg—1 1. There¬ 
fore €g should satisfy equation 6.16, which, with the help of 16.2 
can also be written as 





(16.1) 
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or, ixisertmg numerical values and expressing in Debye units 
(cf. § 6) and T in degree Kelvin, 


= 1.83xiO*A^. (16.2) 

This equation, according to § 6, has been derived on the assump¬ 
tion that the interaction between dipoles is negligible. This 



Fia. 22. Temperature-dependence of the dielectric constant of 

various gases according to Sanger [SI'], 


involves an error of the order of — 1 )*, as can be seen by com¬ 
parison with formulate where the interaction has not been 
neglected (e.g. 6.36). 

Since no assumptions about the interaction have to be made for 
the derivation of 16.1 or 16.2, measurement of the temperature- 
dependence of the static dielectric constant of gases should make 
possible a reliable determination of the dipole moment of free 
molecules. As an example Fig. 22 shows the results of measure¬ 
ments of the dielectric constants of CH 4 , CHgCl, CHaCla. OHCla, 
and CCI4 gas plotted against the inverse temperature l/T (from 
Sanger, SI). Since the highest value of —1 is only about 10 “®, 
the above equations should hold very well. As required, the 
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experimental points lie on a straight line whose slope allows us 
to calculate the dipole moment with the help of 16.1 or 16.2 if 
the density is also measured. 

The results are also of interest in the held of stereochemistry. 
Thus the fact that the dipole moments of methane (CH 4 ) and of 
carbon tetrachloride (CCI4) both vanish (because is indepen¬ 
dent of T) is evidence in favour of a molecular structure in which 
the four hydrogen (or chlorine) atoms form the corners of a 
regular tetrahedron with the carbon atom at the centre. 

Figure 22 also permits a determination of Co, which is obtained 
by extrapolating the straight lines to 1 /T = 0 . The value 
obtained should be slightly larger than the square of the optical 
refractive index n. But, as already mentioned in § 15, the differ¬ 
ence 60 a—is very small for most molecules; in gases it is prob¬ 
ably too small to be measurable using present experimental 
techniques. 

There should be no appreciable frequency dependence of the 
dielectric constant, and hence no loss, up to the absorption 
frequencies of the molecules, which often lie in the infra-red 
region. For a number of molecules absorption begins at the still 
longer wave-lengths of the ultra-short electric wave region. 
Such gases will then show resonance absorption in this region 
following the laws derived in § 13. Relevant experiments will be 
discussed in detail in the book by Jackson and Saxton [«/^]. 
Here we merely wish to discuss briefly the shape of the absorption 
spectrum of ammonia which, following Cleeton and Williams 
\_C2], has been investigated in the centimetre region in great 
detail by Bleaney and Penrose [B3]. Ammonia, NH 3 , forms a 
pyramidal molecule and the absorption is connected with the 
swinging of the nitrogen atom through the plane of the three 
hydrogen atoms. The spectrum shows considerable structure, 
which is resolved at pressures below 5 cm. Hg, as shown in Fig. 23. 
The shape of a single absorption line should be given by equation 
13.11, provided the concept of a single relaxation time r can be 
applied. Assuming 1/t to be pi-oportional to the gas pressure, 
the shape of the absorption sijeotrum at a piessure of 10 cm. of 
mercury has boon calculated by Bleaney and Penrose with the 
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help of 13.11 using data from measurements at 0*5 mm. Figure 
24 shows excellent agreement with the measured absorption. 
At a pressure of 60 cm., however, agreement is no longer so 
good. 

We shall now turn to dilute solutions of dipolar molecules in 
non-polar substances. The concentration will be assumed to be 



Fia., 23. The ammonia absorption speotruin at 1^ 2, 6, and 10 cm. pressure 

according to Bleaney and Penrose 


suHiciently small so that interaction between dipoles may be 
neglected. It then follows from equation 6.19 that the static 
dielectric constant satisfies an equation ainnilfl.r to those of 
gases, but with an effective dipole moment which depends on the 
nature of the solvent and on the structure of the dipolar mole¬ 
cule. This effective dipole moment can be calculated in a simple 
way only when the molecule can be approximated by a point 
dipole at the centre of a sphere with refractive index n. In this 
case equation 6.24 should hold. In general, we may put (cq = 
dielectric constant of solvent; JN'o = number of dipole molecules 
per unit volume) 



4rr/4-^o 

ZkT 



(16.3) 
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where y is unknown except in the case of spherical molecules for 
which, according to 6.24, 


^(^0— l)(^o— 


for spherical molecules. 


(16.4) 


It follows that \y\ is usually much smaller than unity, but, in 
view of the factor —u®, may be either positive or negative. 



Fia. 24. Calculated and measured absorption spectrum of ammonia at 
10 cm. pressure according to Bleoney and Penrose [Bd]. 


The main conclusion to be drawn is that by using experimental 
values for the temperature-dependence of it is not possible to 
calculate the dipole moment fjL„ of the free molecule but only the 
quantity /u.^(l— y) which differs from though usually only by 
a small amount. 

Besides the solvent effect, one of the main differences between 
gases and dilute solutions is their behaviour in alternating fields. 
The interaction between the dipolar molecules and the mole¬ 
cules of the solvent leads to an energy loss whose maximum at 
most temperatures falls within the frequency region of electric 
waves. This affords a new way of obtaining the important 
quantity —€«, without measuring the static dielectric constant 
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*«» for if th© dieleotrio loss, and hence the ima^ginary part 
of the complex dielectric constant e(a>) (cf. § 2) is known as a 

^notion of fiwquenoy, is obtained from equation 2.18, 

which may also be written as 


J €8(a))d(loga>). 


( 16 . 6 ) 

and which holds independently of the special shape of the func¬ 
tion ^^(eu). The integration has to be carried out over th© whole 
range of frequencies in which «a(a>) has appreciable values. The 
present method can be expected to yield more accurate values for 
than measurement of th© static dielectric constant 
bemuse Ca, is nearly equal to th© dielectric constant of the 
so vent and e, is small compared with in the case of dilute 
so utmns. Very high accuracy in the measurement of €g would 

ere ore be required to give reasonably accurate values of 
e. — 


The investigation of dilute' solutions permits a study of the 
a^^our of single molecules independent of dipolar interaction; 
an e results obteined will be of interest in attempts to under- 
s t ® properties of less dilute solutions, or of pur© dipolar 
^ stances. To obtain full information the dieleotrio loss should 
meas over a large range of firequenoies at different 
^mpera urea. Absorption can usually be described by th© 
Debye formulae (§ 10) or by a generalization (§ 12). IVequenoy 

the Debye equatioxS hold^ 

otherL ^ ^ for the relaxation time, but 

AvoU^rmit an estimate of the width of a band of 

thA + ^ measuremwits will then yield 

^e tem^raWde^^dence of the relaxation time-or of the 

^ of relaxation times. No complete 

system of measurements is available at present ^ 

experiments by Jack- 

/I e* X ^ measured at one temperature 

(19 C.) the frequency dependence of dieleotrio loss of dilute 

molecules in benaene and in para^ In 
E^. 26 the expenmental valnee of tan ^ = lo£ migle) for 
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solutions (1 gm./lOO o.c. at 19® C.) are plotted 
against frequency. For the solution in benzene these values 
lie on a curve which represents the Debye formula. According 
to 10.17, 10.22, and 10.23, making use of 10.29, this formula can 
be written in terms of the angrdar firequency at which the 
ma^cimum loss-angle occurs: 

tan<^ = tan^„,- ( 16 . 6 ) 

wh® 37 e tan ^„, == i €«, < 1 - (1®.7) 



10 ’ 10 ® 10 ® 10 ^° 10 ^’ 
Frec^uency in cycles per second 


!Fia« 25* Loss angle of dilute solutions of benzophenone in benzene (maximum 
necLx* 10^<* cycles) and parailln; experimental values according to Jaekson an<l 
Po-wles [t7t3]. The full line represents a Debye curve (equation 16.6, or 16.8 
with p ^ 1) I the dotted line represents equation 16.8 with )3 => 5. 

From 10.22 the relaxation time is then given by t = l/a>^. 

In contrast to benzene the solution of benzophenone in 
panaffin leads to a broader absorption curve than could be satis¬ 
fied by a Debye curve. According to § 12 this can be described 
by introducing a whole range of relaxation times. If, in par¬ 
ticular, we introduce a model where the relaxation times are 
distributed between the two values t 0 and according to the 
distribution fimotion 12.13, then tan^ can he obtained after 
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division by Ci(a>) Cg (of. 10.29) from 12.16 or from 12.18 and 
12.19. Thus with the use of 12.17, 


where 

and 


tan^ = ta>n^^ 

k 


tan~^(<o)3/a>,^)—tan~^(a>/j3aj^) 
tan“^j8—tan~^(l/j8) ’ 


tan 


1 €g —Cqo tan~^/3—tan~^( 1 /j3) 

2 €g logj5 

^ > 1 . 


(16.8) 

(16.9) 

(16.10) 


For p — 1, equations 16.8 and 16.9 are identical with the Debye 
equations 16.6 and 16.7. For j8 > 1, however, the curves are 
broader and flatter. Figure 25 shows that the curve with 6 
satisfies the experimental points fairly well; appreciable devia¬ 
tions occur only at jhequenoies very far £rom eo^ where tan^ 
is relatively small. This is rather to be expected, as follows £rom 
the discussion in § 12 following equation 12.19, because the 
distribution function 12.13 can be considered to represent 
position and width of the distribution of relaxation times but 
not finer details. 

It wiU be of interest to see how the time 1 /o)^ which in the 
Debye case represents the relaxation time is now connected 
with the two limits tq and between which the relaxation 
times are distributed. According to .12.17, 

= V('^o^l) = jSTo = rjp. (16.11) 

Thus, in the present case l/a>„j = 6 tq = tj/S. 

Since no experiments at other temperatures are available at 
present, it is not yet possible to make use of the temperature- 
dependence of the relaxation times to discuss the relaxation 
mechanism. On the assumption of a model in which the relaxa¬ 
tion time is determined by the frequency of jumps of a dipole 
over potential barriers (of. § 11), the factor jS allows us to calcu¬ 
late the maximum difference Vq in height of these potential 
barriers on the assumption that they are equally distributed 
between the values Hq and (of. 12.7). 

Thus from 12.17, using 16.10, 


vJkT = log(Ti/To) = 21og)S. 
In our case, therefore, vJkT = log 26 c::i. 3. 


(16.12) 
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It would be of interest to investigate the variation with tem¬ 
perature of the ratio as well as that of (or t®) itself. 

Measurements at various temperatures were carried out by 
Whiflen and Thompson [TF^], though unfortunately, up to the 
time of writing, the range of frequencies has'been rather re¬ 
stricted. The dielectric loss in solutions of chloroform in hep¬ 
tane seems to follow a Debye curve as indicated in Hg. 26, but it 



Fio. 26. Loss angle of a dilute solution of chloroform in heptane acoording to 
Whiffen and Thompson IW2']; the curve represents the l>ebye formula 16.6. 

must be realized that in the available frequency range there is 
only a variation of a factor of three in tan0. The temperature 
was varied in a range between —70° C. and 80° C. It was found 
that the logarithm of the relaxation times is proportional to 1/T. 
A similar result was obtained for solutes, as shown in h’ig. 27. 
To test the validity of the Debye model (§ 11), according to which 
the temperature variation of r should be the same as that of the 
viscosity rj of the solvent, logrj is also plotted in the same figure. 

Putting (16.13) 


Debye’s model according to 11.32 requires Hr == Theexperi- 


mental values are as follows: 

Hnf. in k,tiaL 

a-Droiiio ^naphthalene 

. 1-8 

Methyl benzoate 

. 1*8 

Camphor 

. 1-7 

Chloroform 

. 1-6 

On the other hand, for heptane H^ = 2*0, which thus is larger 
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jj £qj. gjj th© above molecules. Tliis indicates that in the 
above ^s the I>©bye model probably does not apply exactly. 



2-5 3-0 3-5 4-0 4-5 S-0 5*5 

I/T X 10^ 


Fig. 27. Temperature-dependence of the relajmtion time r of various solutions 
in heptane^ and of the viscosity ^ of heptane, according to WhifFen and 

Thompson [W2'[. 

In all experiments on dilute solutions the value of can 

be obtained from dielectric loss measurements by carrying out 
the integration indicated in equation 16.5. It should be remem¬ 
bered that for dilute solutions both quantities ei(a>), are nearly 
equal to cq, the dielectric constant of the solvent, so that (of. 
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2 . 6 ) € 2 (^) replaced by c^tan^. From, e,—€„,> then, the 

dipole moment oan be obtained according to 16.3 if y can be 
neglected. This method was first used by Sillars [/Sd] in the case 
of dilute solutions of dipolar long-chain molecules in solid 
paraffin-wax. Although equation 16.4 for y is certa.inly not 
applicable in this case, it nevertheless indicates that y should be 
particularly small if the square n® of the refractive index of the 
pure solute is nearly equal to the dielectric constant of the sol¬ 
vent. For long-chain molecules this should be very well fulfilled. 



Fig. 28. Paraffin chain. The indicated distances refer to X*'ray scattering 
centres of CHg groups according to Muller [Jkfd]; they are slightly Afferent 

from C—C distances. 

The investigation of the solid solutions mentioned above 
are of great interest for the elucidation of the mechanism of 
dielectric loss and will, therefoic, be discussed in greater detail. 
We shall start with a discussion of the structure of substances 
built up of long-chain molecules. Most of these structures can 
be derived from that of paraffins, which was investigated by A. 
Muller [Md]. A paraffin molecule, in a crystal, consists of a 
plane zigzag with a CH^ group at each comer and a CH 3 group at 
each end (of. Fig. 28). The distance between neighbouring 
CHg groups is about 2 A (1 A = 10 "'®cm.), while the distance of 
their projections to the chain axis is about 1 *25 A. In a paraffin 
crystal the chains are arranged in layers whose thickness is 
approximately equal to the chain length. Within such a layer 
the molecules form rectangular cells with side length a, 6 , c, 
where a 6 A, 6 crL 7*6 A, and c is slightly larger than the chain 
axis and is parallel to it. Figure 29 shows how the chain planes 
cross the a —6 plane which is perpendicular to the chain axes. It 
is of importance to notice that in the subsequent layer the whole 
arrangement is shifted by about 1 A in the direction of the 6 -axi 8 . 
This is indicated in Fig. 29 by the dotted lines. 

A paraffin molecule carries no dipole moment, as has already 
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been mentioned in § 16. Now suppose that some dipolar It 
chain molecules, e.g. esters or ketones, are dissolved in a pare 
or 3 ?Btal. This means that some of the paraffin molecules 
replaced by dipolar long-chain molecules. These latter can 
usually derived from paraffin molecules. A ketone molecule, 
instance, is obtained from a paraffin molecule by replacing t 
CHg group by a CO group. This latter carries a dipole momt 
whose direction is perpendicular to the chain axis and proba 



Pio. 29. Flu. 30. 

Pio. 29. Seotion of paiafOn crystal perpendicular to chain axis indicating th 
chain planes, according to Hnller [Jif fi]. The dotted lines indicate th© sectioi 

of c hain planes in a subsequent layer. 

Fio. 30. The two equilibrium poaitiona of a ketone molecule^—jndioating^ 
dipole directiouB—leplaoiag a parafiSn molecule of longer chain length. Th* 
suse of the xeotangular cell is the same as in 29. 

lies in the chain plane. Also assume that the length of the ketone 
molecule is shorter than that of a paraffin molecule so that it 
can easily replace a paraffin molecule. There should then be two 
possible directions for the dipole which differ by 180^, as indi¬ 
cated in Kg. 30. For, if the ketone molecule is shorter than th© 
paraffin molecule, the former should still fit into th© structure if 
it is turned through 180“ and shifted hy the length of on© link 
along the chain axis (of. Pig. 31). In a dilnte solid solution such 
a dipolar molecule has thus two equilibrium positions with 
opposite dipole directions. Prom the similarity of the two 
positions it can be assumed that the energy of the molecule is 
the same in both. 

In general, a molecule will oscillate about either of the two 
equilibrium positions. Occasionally, due to a thermal fluctua¬ 
tion, it will acquire sufficient energy to turn into the other 
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equilibrium position. Such asubstance thus may well be described 
t»y the high-temperature model of dipolar solids discussed in 
§ 11, where a dipole was supposed to have two equilibrium posi- 
-bions separated by a potential barrier; the TniTiiTmiTn energy 
required to lift it over the barrier was denoted by H. The energy 
loss of such a substance is described by the Debye formulae 
(e.g. 16.6 and § 10) and the relaxation time is given by equation 
11.3 (cf. also 11.8). 

Loss measurements on dilute solutions of the type described 
above have been carried out by W. Jackson \J1, 2], SUlars [S6^, 



lE'io. 31. Result of turning a ketone moleoule by 180° and ^biftiiig it by one 
chain length. The aorow indicates the dipolar CO group. 

I 

and Pehnore [P2]. These authors find that the dielectric loss 
oturves satisfy the Debye equations fairly well, though not exactly. 
This slight deviation from the Debye equations may be due to 
bhe fact that a dipolar molecule which is shorter than a para£G.n 
molecule, say by z links, has z possible locations (apart from the 
two dipole directions) corresponding to different positions on 
the c-axis. Its energy in these positions may not be exactly the 
same, and the same may be true of the energy H required to lift 
a. molecule over the potential barrier. This would lead to a distri¬ 
bution of relaxation times. It should be realized that in calcu¬ 
lating the relaxation time t according to the Debye formula 
10.22 some error is introduced (cf. 16.11), though it may be negli¬ 
gible if the variation of logr with temperature, and not the 
absolute value of t, is considered. 

The most interesting result of these experiments is the 
measurement of the dependence of the relaxation time on the 
chain length of the dipolar moleoule. If the moleoule is con¬ 
sidered to be rigid, and if therefore each link has to be lifted over 
the barrier at the same time, the total energy required to lift 
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a whole molecule of nhq.iTi length m over t/li© barrier shotild 
increase in proportion to m. Thus, sinoe aooording to 11.3 

logT = oonstant+jOJ^^iT, (10,14) 

the experimental values of log t plotted against vn should lie on 
a straight line. This is not the case. Actually molecules are not 
rigid but have a certain flexibility, and, if this is taken into 
account (cf. reference satisflaiotory agreomcnt with experi** 
ment is obtained. The influence of the flexibility on the value of 
can easily be discussed qualitatively. Suppose for a moment 
that the chain is rigid. Then the whole nbAin bias to be lifted over 
the j^tential bamer at the same time. Thus if H-y is the onei^gy 
squired to lift a smgle link over the barrier, the total enoi^gy 

be equal to where m is the jnumber of carbon 
- -tM. Actually, however, the diain has a certain 

allows the chaiu to be gradually lifted over the 

n» ^ energy less than wiSTj is required. 

^ f'^wistiug of the molecule requires a certain 
fiir «. ofv ^ energy decreases with increasing chain length 
for *^^0 ends of the chain. Thus 

ecMMtant val • j'** for long chains will tend to a 

®*‘*®«**>**on length. A simple 

^==^i»*otanh~, ( 16 . 16 ) 

t^parison with cxiSimeTife, u ^ ^-PProoiable torsion). 

= 28, both «wnrt,Ie fi ~ 

cowrtwt is required acoortW to eoua+f^ obtain r a further 

expenmentally determined oonstAT,^ 10.14. Thus three 

formula fr in ‘seconds) included in the final 

log2m- = — 50 - 4 -l 1 _^^ 00 . , m 
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only if w*" is smaDer than the chain length of the paraffin in which 
the dipolar chains are dissolved. Figure 32 shows that the results 
obtaiii©<3- from equation 16.16 agree well with the experimental 
results. The expression tanh(m/26) is plotted as a frmotion 
of wi/26 and is compared with the experimental values of 
(log —{-60*4) r/13800 for four different chain-lengths (m — 20, 
22, 24, 32). Two of the constants in equation 16.16 have 

been deisermined from the experimental values for m = 20 anti 



Fio. 32. Oomparison of the theoretioal depeodenoe of the relaxatioD. timA on 
olnain length with experuneintal values^ aooording to [^5]. 

m = 22, while the third was determined from the temperature- 
dependence of m = 20. For m =- 24 and 32 no more adjustable 
oonstarLijs are available, which shows that formula 16.16 leads 
to a oorx:eot dependence of the relaxation time on chain lengths. 

Frona the numerical value —50*4 for the constant in equation 
16.14 (oomparing with 16.16) it might appear that the order of 
magnitixde of the constant A of equation 11.3 (cf. also 11.8) could 
be estimated. For asstuning for the frequency of oscillation 
(aJ^TT 10'* one finds log .<4 —26. Not much significance 

should lae attributed to this value, however, in view of the un- 

« 

certain'bies connected with the experimental determination of 
the absolute value of t; for, as mentioned above, the deviation of 
the loss curve from a true Debye curve brings with it some un- 
certain.'bies in determining t which cannot be adjusted by apply¬ 
ing equation 10.22 as in the Debye oase. This would probably 

4B80.1X „ 
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influenoe the value of ^ to a much greater extent than the less 
sensitive 

In aJl examples discussed above the concentration of dipokr 
molecules was assumed to be so small that dipolar interaction 
wotdd be neglected. In this range, therefore, the quantity 
€g—Coo as "well as dielectric loss should increase proportionally 
with the concentration. For farther increases of concentration 
the effects of dipolar interaction should start to appear, and an 
investigation of the deviations of —e* froni equation 16.3 
would be of great theoretical interest. In particular, it would 
be of interest to see to what extent and in which range of con¬ 
centration the Onsager formula 6.38 can account for such devia¬ 
tions. This would be of particular interest for a solute which in 
the pure state does not satisfy the Onsager formula for the 
temperature range contemplated.t It will be remembered (§ 8) 
that the Onsager formula tahes into accoimt the long-range 
dipolar forces, but neglects short-range forces. Then with 
increasing concentration we should expect three stages: (i) all 
interaction can be neglected, leading to equation 16.3 for (e,— e > 
(ii) long-range dipolar forces have to be considered, but short- 
range forces can still be neglected, Ifiading to the Onsager formula 
6.38; (iii) short-range forces have to be considered as well. This 
leads to the Kirkwood formula 8.14. 

Of very great importance would be the investigation of the 
dependence on concentration of the shape of the dielectric loss- 
frequency curve. The deviations from a Debye curve would 
probably increase with increasing concentration owing to the 
interaction between dipoles. As yet no theoretical investiga¬ 
tions on this subject have been carried out. 

17. Dipolar solids and liquids 

Survey 

From the developments of §§ 7 and 8 it follows that no formula 
exists which ropresents in a simple manner the dielectric 

f Measuremonta on the dependence of cg on concentration have been 
carried out, of courae. They do not extend^ however, to the region of very low 
concentration where eg—-e® is best obtained from loss measurements with the 
help of 10.5. 




181 


I 17 DIPOLAB SOIilDS AND DIQtTIDS 

bebavioiir of dipolar solids and liquids. Qualitatively, however, 
most of these substances behave in a very similar way. At very 
low temperatures in the solid state all dipoles are frozen in, and 
if the solid has not become permanently polarized the dipoles will 
not ooiitribute to the dielectric constant. In this region (T 6), 
therefore, the dielectric constant e,(0) should be nearly inde¬ 
pendent of temperature, and its magnitude should be approxi¬ 
mately equal to that of the dielectric constant e^iiT) at a higher 
temperatxue but at a frequency which is so high that the dipoles 
do not take part in the polarization. This equivalence holds 
only, of course, if the voltime of the substance is the same at both 
temperatures, and it includes the assumption that upon turning 
a molecule from one equilibrium position into another its con¬ 
tribution to the high-frequency dielectric constant is not 
altered. Then 

€g{0) c:i 6oo(T) = 9i*- 1-A€, constant volume. (17.1) 

In t-bia formiila it has been indicated that €„ is composed of 
the square of the optical refractive index n (cf. § 15), and of a 
term Ae due to elastic displacement of nuclei or of dipoles. The 
corresponding polarization was called infra-red polarization in 
§ 14. A a pointed out in § 15, this typo of polarization is frequently 
negligibly small in non-dipolar substances. In dipolar substances 
Ae usually has small but noticeable values, mainly due to a dis¬ 
placement of the equilibrium directions of the dipoles by the 
field. The frequency of the rotationary oscillations which the 
dipoles carry out about their equilibrium positions is normally 
in the far infra-red region, but there are indications that in large 
molecules it may approach the centimetre region (cf. reference 
F9). Such substances would thus show resonance absorption 
(§ 13) at ultra-high frequencies. Substances which are expected 
to show such absorption have been discussed by Szigeti 
and a further theoretical investigation on the expected fre¬ 
quency-dependence of the absorption has been carried out by 
Huby leading to a generalization of equation 13.11. In 

this connexion it should be mentioned that Girard and Abadie 
[G3] found that normal liquid long-chain alcohols absorb in the 
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centimetre region besides showing Deby© absorption at longer 
wave-lengths. It was suggested by Magat [-ifI] that this is an 
example of rresonanoe absorption. Further experiments would 
be required to establish the validity of this interpretation. 

Let us now turn to the regions of higher temperatures and 
lower hrequencies where dipole oiientatioii plays an important 
role. Within creasing temperature an increasing number of 



!Fia. 33. Temperature -dependence of the dieleotrio oonstant of solid hydrogen 
sulphide aooordmg to Smyth and Hitohoook [891 sliowing an order-disorder 
transition near 103^ abs. A further transition and melting are of little induenoe* 


dipoles oftTi turn into other eqriilibrium directions, as pointed 
out in § 8. The static dieleotrio constant will thus increase 
with temperature, at first according to the relation 

oc (17.2) 

as follows from 8.17 (F(0) is a constant energy discussed in § 8). 
At higher temperatures the rate of rise of gets increasingly 
larger until at the temperature the solid carries out an order- 
disorder transition beyond which (T >■ 2^) the dielectric oon¬ 
stant falls with increasing temj>erature. The order-disorder 
transition may occur in the solid state; a number of examples are 
shown in Figs. 33 and 34. In this case there is usually no appre¬ 
ciable change ofat melting. In other substances melting occurs 
before the order-disorder transition is completed. In these 
oases the decrease of starts at the melting-point. 
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Ab shown in the discussion of entropy in § 3, a positive value 
of deJdT indicates that the degree of order is d.eoreased. when 
the field is applied, while a negative value of deJdT indicates the 
opposite, namely, that the degree of order is increased by the 
application of the field. It follows that in a substance in which 
the dipole directionB are completely disordered BeJdT must be 
negative, whereas if the dipole directions are completely ordered 



Fia. 34. Temperature'dependence of the dielectric constant of solid ethylene 
cyanide according to White and Morgan [W4] showing an order ^disorder 
transition near —4CP C. without a discontinuity at melting. 

deJdT must be positive. ^riiuB from purely thermodynamical 
reasoning it follows that the transition from positive to negative 
values of deJdT must be connected with a change in the order 
of dipole directions. A large number of examples of such transi¬ 
tions has been reviewed in a recent article by Smyth [/Sfy]. 

Dielectric loss due to dipole orientation is expected to occur 
at T >• as well as for T <iTf. (Figs. 36 a and 6 give examples 
for the latter). Fxcept very close to 2^ the latter will, of course, 
be much smaller than the former because of the much smaller 
value of e,—Thus (of. Fig. 36) at any temperature we expect 
resonance loss in the far infra-red or ultra-high frequency region, 
and loss of the Debye type (but requiring a distribution of 
relaxation times) at longer wave-lengths. For the latter, no 
satisfactory theoretical treatment of the shape of dielectric 





-150 -100 - 60 0 +20®C 



J°“ * '• di--»p«pyi 

-^at4-4Mo 1 * frequency of 1-2 Mo./u., 

'»~*~-a<^20Mo./8. Melting occurs near -73° C. tuid tlioro w 
no previous transition. 
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loss-frequenoy curves exists at present. Compared with dilute 
solutions the difGLculty consists in accounting for the influence 



Pia. 36* Frequenoy-dopendenco of dieleotrio loss,-abovo and- 

'below the transition temperature (schematioally). wol2n represents the 
lowest resonance firequenoy which is usually in the far infra-red or in the ultra- 

high electric frequency region. 



Fig. 37. Frequenoy-dependonoe of loss angle according to Hartshorn, Megson, 

and Kushton [^i'J for various materials. 

(1) Phenolic resins. (2) Benzyl alcohol resins. (3) Rubber-sulphur compounds 
(Scott, McPherson, end Curtis, S5)- (4) Cetyl palmitato in paraffin wax 

(SiUars, S7), (6) Chlorinated diphenyl (Jackson, J2). (6) Represents a Debye 
curve; comparison with Fig. 14 shows that curves (l)-(5) can be represented 
by equation 12.19 with suitable values of the parameter ti/tq. 

of iuteractiou. This, one should expect, would lead to a broaden¬ 
ing of the Debye curves similar to that obtained through the 
existence of a whole range of relaxation times (of. examples in 
Fig. 37). The width of the latter ought to be larger in amorphous 
substances than in crystalline solids. It would be of interest, 
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theiefoTe, to investigate substances which at the same texnpem** 
ture exist in various states of order; e.g. amorphous substances 
which have been quenched from different temperatures, or 
substances existing at the temperature as super-cooled 

liquids and oa crystalline solids. 

Anc^aia of eg~^T marve^ 

Without any knowledge of the structmre of a dieleotrio, in¬ 
teresting oohotujnons about the behaviour of its dipoles can be 
drawn from the experimental Cg-f-iT curves with the help of our 
general theory. For from equation 7.39 it follows that we can 
obtain the quantity 

= (17.3) 

O fC 3 €^ 

proTided tlxatthe optical irefipaotive index known in th.B correct 

region (of. § 16), If this is not the case, we can make use of the 
fact that the infca-red contribution Ac = is usually small 

(except in ionic crystals) compared with and will replace 
by Cod and hence, according to 17.1, by €,(0). Thus, if 


we find approximately 

S(r) ~ fe(T)-..(o)}?s^±^> y. 

^ fhis v^alue 1.7.6 foT Hi'FFcvpa 4VrMfv\ 


(17.4) 


(17.6) 


/ITA,' ^*(y) + 2e,(0)-A€ 

SCa 


(17.6) 


which is of the order iTAc* Using approximation 17.6, it is thus 
possible to obtain S{T) from, measurement of the static dielec¬ 
tric constant at various temperatures. Theoretically 
according to 17.4, is proportional to xnm^, and the value given 
by 17.6 is thus constant (subject to the small correction 1*7.6) 
independent of temperature, if Onsager’s formula holds, i.e. if the 
dipole moment of a molecule can be considered as a constant and. 
if there are no short-range fcrces leading to mutual orientation 
of dipoles. Since for T ->oo all interaction can be neglected 
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and the dipole orientation is completely random, a function 
JB{T) which increases with decreasing temperature points to an 
increasing tendency towards parallel orientation of dipoles. 
Similarly if B{T) decreases with decreasing temperature the 
dipoles tend towards anti-parallel orientation. It should he 
pointed out again (of. §§ 7 and 8) that this orientation is due to 
non-dipolar short-range forces. Long-range dipolar forces too 
lead to a certain tendency for dipole orientation. But averaging 
over all angles between dipole direction and radius vector, this 
orientation cancels in an isotropic substance. Since m* is the 
average moment of a sphere if one of its rmits (molecules) has 
the moment m, dipolar interaction alone leads to in’* = m, 
i.e. to a value of mm** which is independent of temperature if 
m is a constant. 

If the assumption is no longer made that the dipole moment 
m of a single unit is constant, then similar conclusions to those 
above can be drawn &om the temperature-dependence of B{T), 
A decrease of B{T) as the temperature is lowered can now mean 
either an increasing tendency towards anti-parallel orientation 
between the dipoles of neighbouring units (molecules or unit 
cells) or a decrease of the dipole moment of a single unit^ or 
vice versa if B{T) increases. 

In Figs. 38 (i) and (ii) a few examples of this analysis are 
given. Wo shall now show how in the case of a known structure 
the general theory can be applied. 

Water [if4, 02, K6] 

Water no doubt is one of the most important dielectrics, and 
the correct calculation of its dielectric constant must be con¬ 
sidered as a great success of Kirkwood’s [if4] formula 8.5, 8.14 
which applies to dipolar liquids in general. A water molecule can 
be considered as a negative O ion with the two positive H”*" ions 
attached in such a way that the lines connecting the latter with 
the centre of the oxygen ion form an angle of 106°, as can be con¬ 
cluded from the infra-red absorption bands of H 2 O vapour. The 
dipole moment of the H 2 O molecule is a vector directed along 
the lino dividing the H—O—bond angle into two equal angles 
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of 52^°. In the liquid state the position of an individual molecule 
is strongly correlated with those of its neighbours, and X-ray 



38, Ebcamples of €g—T azLalysiB. (ii) Temperature-dependezioe of the 
function JB(T) obtained &oxn. CgC?) (of, i) with the help of equation 17.5. 
(o) Diohloro-propane* (OHa)xCClt» solid, using €^(0) = 2-23. (6) Tertiary- 
butyhohloiide, (CH,)3CC1, liqxiid, using €,(0) = 2'46, (o) Ponta methyl- 
ohlor-benzene, (CHs^^Cl, solid, using » 2-8. Measurements: (a), 
(6) Tucgewioh and Smyth \T2'\i (c) White, Biggs, and Morgan [WS}* 
With decreasing temperature, (a) indioates increasing tendency to 
anti^paraUel—(&) to parallel—orientation of dipoles., (o) points to 
random orientation, and to validity of Onsager’s formtda. 


investigations show that the average number of nearest neigh- 
hours is close to four. Also, X-ray investigations of liquids show 
that in general the nearest neighbours of a molecule are arranged 
in a fairly ordered way corresponding to a certain type of crystal 
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struotuTS. In. contrast to solids, however, this order grardually 
disappears for more distant molecules. In water, acoordmg to 
the model of Semal and IFowler the four nearest neigh- 

hoiurs of a given molecule form a nearly regular tetrahedron with 
that molecule at its centre. The bond 
molecules, called the hydrogen bond by 
chemists, is supposed to be directed 
along the O—bond of one moleoiile 
towards the oxygen ion of the other. 

The arrangement of a molecule and its 
four neighbours may be represented 
schematically as in the diagram, (Fig. 

39) in which the hydrogen bonds are 
indicated by dotted hues. 

To simplify the following considera¬ 
tions we shall assume the bond angle 
to be 2oos“^ 1/V3 Cii 109° instead of 
106°. This should make little differ¬ 
ence to the final result. It enables us to 
consider the tetrahedron formed by the 
four neighbours of a given molecule as regular, with that mole¬ 
cule at the centre. A simple way to picture the arrangement of 
molecules and their dipoles is to start with a cube with four 
of its vertices forming a regular tetrahedron, and the centre 
occupied by HgO molecules (Fig. 40). The four bonds emanating 
from the centre, namely two —O—H* • -O^ bonds and two 
>0 • • • H—O— bonds, follow the lines joining the centre with 
the vertices. A dipole can have any one of the six directions 
perpendicular to the faces of the cube. Once the direction of 
the central dipole is fixed, however, each of the neighbouring 
dipoles can have one of only three possible directions: for the 
two molecules bound by —O—H • • • bond the dipole direc¬ 
tions are away from our cube, and for the other two towards 
it, as can be seen from the figure. The other three directions 
are excluded because they would require —O—• -H—O— 
bonds. It will now be further assumed that correlation with 
more distant neighbours need not be considered. Then the three 


between neighbouring 



Pia. 89. A H^O moleonle 
with its four neighbours, 
BohematioapUy. Full lines 
represent ordinary chemical 
bond, dotted lines indicate 
hydrogen bond. 
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possible dipole directions of the neighbours wiU be equally 
probable. Hence the average value cosy of cosy, where y is 
the angle between neighbouring dipoles, is 

cosy = (l*^**^) 

because in two positions the dipoles are perpendicular to each 



Fig. 40. An HjO moleoule with its four neighbours; dotted lines represent 
hydrogen bond. The arrows at the comers show the possible dipole direotiona 
of the neighbours for a fixed dipole direction of the central molecule. 


other (cosy = 0) and in one they are parallel (cosy =1)- Tii© 
same result can be obtained if the O—bond outside the centre 
is permitted to have any direction which does not destroy the 
O—^H* • *0 bond; in other words, if free rotation around the 
O—H* ■ *0 bond is assumed (of. Kg. 41). For this leads to an 
average value pL = in the direction of the O • • • H—O 

bond, if is the dipole moment of a moleoule, and hence to a 
value ft/VS = /x.^/3 in the direction of the fixed dipole. 

Upon inserting the value of cosy from 17.7 into 8.14, using 
2 = 4 for the niunber of neighbours, we find since 


l+zoosy = 1+4/3 = 7/3, 

^ _ 3€, / n»+2 \g ^ ^ 7 

“ 2e,+TOa\ 3 / 3 3* 


( 17 . 8 ) 


or since the refractive index n = 1*33, i.e. 


3e,/(2e.+w*) c- 3/2, 

* ~ 3 \ 3 I kT ' 


( 17 . 9 ) 
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Now substituting the following numeTica.1 values: Avogadro’s 
number =6x10**, density = l, molecular weight = 18, 
^ == l*4x lO”'-*, = 1*9 X 10"'^®, n® = 1'77, we find, using 16,2, 



19000 

T 


i.e. 


Cg E± 63 


for T = 900® abs.. 


(17.10) 

which compares favourably with the experimental value of 78. 



H 


0 


Pio. 41. Direotion of tlie dipole relative to an O—bond. 

At 80® C., i.e. at about 363® abs., the expeidmental value of 
€g is 60 and the theoretical one is 63. This shows that the ex¬ 
perimentally determined temperature-dependence is also in &ir 
agreement with the l/I* law given by 17,10. 

In judging the approximations made above it must be remem¬ 
bered lhat the rigid structure assumed for a molecule and its 
neighbours does give only an approximate picture of the correla¬ 
tion between neighbours. X-ray investigations show that 
the average number of neighbours is slightly larger than four; 
this may be the reason for the larger experimental value of e,. 
Also, with increasing temperature the restriction placed on the 
relative orientations of neighbouiring dipole by the hydrogen 
bond may occasionally be overcome. This would decrease the 
dielectric constant and account for a stronger temperature- 
dependence than that given by equation 17.10. 

It thus seems that on the whole the Kirkwood formula 
accounts in a satisfactory way for the static dielectric constant 
of water. 

No satisfactory theoretical treatment of the firequency- 
dependence of the dielectric properties of water is available at 
present. This would require a generalization of the Kirkwood 
formula applicable to time-dependent fields. The experimental 
results by Collie, Hasted, and Kitson \GS\ and by Saxton and 
Lane point to very simple properties. It seems that the 
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Debye formulae 10.16-10.17 are satisfied, and moreover that 
the relaxation time t is conneoted with the viscosity by Debye’s 
relation 11.30. This means that t and ij/T have the same tem¬ 
perature-dependence; hence t must be a linear function of tj/T. 
iPigure 42 shows that this holds for water as well as for heavy 
water. 


S 

T 

2 

1 


Fio. 42. Dependenoe of relaxation tizoe t on ij/T (17 •>« visooBity) for water 
fkiid heavy water, according to Collie, Hasted, and Ritson £ 03 }. 

Ketoves [J’S] 

As an example of long-chain substances we shall consider solid, 
ketones, and we shall begin with a discussion of their structure 
in the lowest energy states which is realized at sufficiently low 
temperatures. A ketone molecule is obtained from a parafifixi 
molecule by replacing one (or several) CH 2 group(s) by one (or 
several) CO group(s). In a solid a ketone molecule forms a plane 
zigzag, as shown in Fig. 28, § 16. The dipole lies in the chain 
plane and is directed perpendicular to the chain axis. In the 
simplest case the crystal structure is similar to that of paraffins 
which has been investigated by Muller [Afd]. The molecular 
chains are arranged in layers whose thickness is approximately 
equal to the chain length. Within such a layer the molecules 
form rectangular cells with side length 

a ci 6*10“® cm., b 7*6 x 10“*cm., 
and c slightly larger than the chain length. The dipoles of such. 
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a layer are all in one plane, the dipolar plane, and they are 
arranged as shown in Fig. 43. These dipolar planes are thus 
polarized in the 6-direotion, and it is of great importance in 
studying the behaviour of the whole crystal to know the relative 
directions of the polarization of successive dipolar planes. First 
then consider the position of successive layers in paraffins. 
According to Muller, distinction must be made between two 
cases, depending upon whether the number of C-atoms in a 
oh ain is even or odd. Fig. 44 shows the positions of two successive 



Fio. 43. Dipole direotiops in tbs gcomid state of 

a dipolar plane. 

layers for the two oases. Now let us attach a dipole to each 
molecule at distances and Cg from the two ends so that c^-l-Cs is 
equal to the chain length. For both even and odd chaias there 
now exist two possibilities. Except for the gap between the 
chains the distance between successive dipolar planes is (a) 
C 1 +C 2 , or (6) it is alternately 2ci and Figure 46 shows that 
for odd chains both cases lead to opposite directions of the 
polarization of successive dipolar planes. For even chains, 
however, this is only so in case (6), whereas in case {a) successive 
layers have the same direction of polarization. This latter case 
will thus lead to a strong polarization of the crystal with aU its 
dipoles nearly parallel, while in the other cases the polarization 
of successive layers cancel. 

For paraffins both cases, (a) and (6), are identical. Their 
energies, therefore, differ only by the contributions of the 
dixJoles. 

The interaction between dipoles can be considered as composed 
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of the interaction between the dipoles of each single layer and 
of the interaction between dipolar layers. Since the distance 
between neighbouring dipolar layers is large compared with the 
distance between neighbouring dipoles within a layer, it is possible 
to approximate a dipolar plane by a continuously polarized plane. 
The interaction of a given dipolar plane with all the others is 



odd even odd even 

Fig. 46. The two possible positions of 
Fia. 44. Position of pacafiSn chains dipoles in the groixnd state of odd and 
in suooeesive.layers (schenxatioally)* even ketones (sohematioally)* 

then equal to its interaotion with the surface charge which the 
other planes produce at the surface of the specimen. iFor a 
sufficiently large specimen this interaction energy should be 
nearly equal to the self-energy of a continuously polarized 
specimen of the same size and with the same total moment. 
According to the appendix (A 2.iii) this energy is positive, but 
its value depends on the shape. It thus follows that the energy 
of the spontaneously polarized structure is higher than that of 
the other structures, so that we cannot expect the spontaneous 
formation of permanently polarized crystals, corresponding to 
feiro-magnetio substances in the magnetic case. On the other 
hand, if Such a structure oould be obtained in an even hetone, 
it would probably persist for a very long time, because to reach 
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an unpolaiized state the chains of some dipolar planes would be 
required to reverse the direction of their chain axis; this should 
hardly be feasible in the solid state. It is possible that the 
required structure can be obtained by solidifying an even ketone 
in the presence of a strong electric held. The strength of the held 
required depends on the shape of the specimen; the most favour¬ 
able case would be that of a needle-shaped specimen with its 
pa'»*all el to the held, because such a shape leads to the 
smaUest self-energy. 

Now let us consider the dielectric properties of ketones which 
are not permanently polarized. At low temperature the only 
effect of the held will be to displace slightly the equilibrium 
direction of the dipoles, apart from the displacement of electrons. 
By comparing, at these temperatures, the static dielectric con¬ 
stant of a ketone with that of a paraffin of equal chain lengtii, 
Muller [ATT] found that the contribution of dipole displacement 

is of the order Ae 01. (17.11) 


This value can be used to estimate the frequency of the 

rotational oscillations of the dipoles. Suppose that the held E 
an angle 6 with the equilibrium direction of a dipole ft. 
It then exerts a couple y^E sin 6 on the dipole. The restoring 
couple is Itxil 0 if J is the moment of inertia and 4 the average 
an gnlii.T displacement of the dipole direction towards the held 


direction; hence 


j fLEsiixd 

* -Tar" 


(17.12) 


For weak helds, ^ < 1, so that fi^sinO is the projection of the 
induced dipole in the held direction. Averaging over aU dnw- 
tion between p. and E, i.e. replacing sinafl by f, leads to an m- 

duced moment per unit volume 

2 y^ENg (17.13) 

— 3 Ja>* 

if NJ, is the number of dipoles per unit volume. Hence using 1.9 
the contribution to the dielectric constant is 

(17.14) 

3 Ioj% 

For the ketone used in MuUer’s experiments W. =! 2 X 10» per 


Ac = 


4980.11 


Ij 
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0 . 0 ., jx 2*6 X10““ e.B.u. To obtain I we first assuiae that the 
C-atom remains in its initial position and the O-atom can be 
displaced. Then I — where A = 16x2x 10“®* gramme is 
the atomio mass of oxygen, and r 10“® cm. Hence using 
17.11 and 17.14, the resonance wave-length Aq is of the order 



2vc /SJ A€ \\ ^ 

/X \8nNol - 


(17.16) 


On the other hand, if we were to assume that the chain oscillates 
rigidly, then 1 would be larger by a factor of the order of the 
chain length {r^ 20), and hence Ao would be larger by about •V20. 
Thus we expect a resonance wave-length between 1/10 and 
1/100 cm. This may well lead to a measurable absorption in the 
centimetre region. 

It was shown in § 16 (of. Kg. 30) that for dilute solutions of 
ketones in a paraffin crystal with longer chain length there are 
two equilibrium positions for each molecule with opposite 
dipole direction and about equal energy. In a pure ketone crystal, 
too, it seems likely that an equilibrium position can be obtained 
by turning a chain plane by 180° around the chain axis. This new 
position must, however, have a higher energy than the former 
one, partly due to the alteration of the interaction, between the 
dipoles and partly due to non-dipolar interaction. The inter¬ 
action between dipoles is mainly restricted to those of the same 
dipolar plane. Hence, in view of the strong anisotropy of ketone 
crystals, calculation of this interaction is reduced to a two- 
dimensional problem. In this case the interaction with nearesb 
neighbours is a good approximation to the total interaction. 
Since the non-dipolar interaction is also restricted to nearest 
neighbours it is not necessary to take into account any long- 
range interaction if we consider non-polaiized states of the 
crystal. It follows from these considerations that with increasing 
temperature an increasing number of dipoles will turn in the 
opposite direction, leading to an order-disorder transition in the 
way discussed in § 8 and in the first section of the present para¬ 
graph. The dielectric constant should, therefore, rise with tem¬ 
perature, as indicated in Fig. 10, imtil the transition temperature 
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has been reached. Actually Muller’s experiments show that the 

crystal melts before the transition temperaturehasbeenreached-t 
This possibiUty was envisaged in our previous discussion. 

A transition of the type considered here was found, however, 
by Ubbelohde [UI] in a paraffin by measurement of the specific 


O 

J_600 

“to 

o 


o300 



-60 -40 -20 ^ . 40 "C 

heat c at constant pressure. The turning of molecules into 
^Lnt equilibrium positions mak^ a contribution Ac, to 

. -L -■■-! nf nrder-diaorder transitions. It leads to 

which IS charactenstio or oraei 

a inorease of o„ noar 

by an abrupt fall beyond it. Rgure 46 *or™ that Ac, oan be 
aeparated &om the normal apemflo heat ftn aoonr^y. 

The tranaition U oonneoted with a Aange dS m entropy which 
ia given by the well-known thermodynamio relation 

(17.16) 


AS 


= J 


Ac 


T 


^dT 


integrated over all temperature. If dS ^ the entropy merea^ 
per Lleeule, then in the diaorde^ a m^eoule he 

^tions for each one in the ordered atate. We should thus 

exp60t ^5=fclog2. (17.17) 

t Tto immltion h»i aev tee O'- 

paJaffin-kefcoae solution. Added in pro f. 
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Actually, however, the experimental value of AS is log 20- 

This means that above the transition temperature a chain must 
have many more than two positions. It is likely that this is due 
to twisting of ohains [jFd]. This may be expected in view of 
the flexibility of chains which has been demonstrated previously 
in the discussion of the relaxation times of dilute solutions of 
ketones (§ 16). 

Additional evidence of the twisting of chains near the transi¬ 
tion temperature has been given by Muller He measured 



I _ i _ i _ I_L--1 

10 25 40 55 70 

Pia. 47. Temperature-depeudenoe of the dielectric constant of the two di¬ 
ketones OioB[i80a(l) and CiiHsoOs(2) according to Muller [JlfS]- mole¬ 

cular chainB in the untwisted state are shown as zigzags, the chtowb indicating 

the directions of the dipolar ketone groups. 

the dielectric constant of ketones containing two ketone groups, 
i.e. two equal dipoles, for two different cases: (i) the number of 
links between the ketone groups is odd, and (ii) it is even. In 
case (i) the dipoles of the ketone groups are anti-parallel in the 
untwisted chain, so that the total dipole moment vanishes where¬ 
as in case (ii) they are parallel. Nevertheless, in this case the 
dielectrio constant showed a similar increase below the melting- 
point to that of case (i) (of. Kg. 47). This is possible only if chain 
twisting occurs in this temperature region, for otherwise the 
ketone in case (i) would behave as a non-polar molecule and 
have a dielectric constant independent of temperature. 
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18. Ionic crystals 

Ionic oxystals are materials in which each lattice point is 
occupied by an ion. Examples are the alkali halides such as 
rock-salt, iNTaOl, which forms a simple cubic lattice consisting of 
two interlocking face-centred cubic lattices of Na+ and of Q- 
ions. Normally the negative ions are much larger than the posi¬ 
tive ones and often the negative chaoge of the former overlaps 
into the region of the neighbouring positive ions so that it is 
not always possible to refer to the charge of an ion as an integral 
multiple of ±e. In the case of alkali hedicids the success of Bom’s 
[B4'\ lattice theory suggests, however, that the charges of posi¬ 
tive and negative ions are well separated. For in this theory the 
binding energy and other properties of the crystals are calculated 
on the assumption that the interaction between ions is composed 
of (i) the attraction between the ions of charge ie, (ii) the repul¬ 
sion between nearest neighbours, and (iii) some corrective terms 
due to van der Waals forces. 

The polarization of ionic crystals is entirely due to elastic 
displacements, in both the optical and the infra-red regions. 
Optical polarization, as indicated in § 14, is due to displacement 
of electrons relative to the nuclei and the corresponding reso¬ 
nance frequencies are in the visible or ultraviolet region. Infra¬ 
red polarization is connected with the displacement of nuclei, 
and usually is accompanied by a displacement of electrons 
relative to the nuclei. In contrast to its relatively minor part in 
other substances, infta-red polarization plays an essential role 
in the case of ionic crystals. We shall, therefore, have to consider 
it more closely than has been done hitherto. The most general 
(not necessarily homogeneous) polarization can be considered as 
a superposition of plane polarization waves whose wave-lengths 
range from the order of the lattice distance to infinity, the latter 
corresponding to homogeneous polarization. The frequencies 
of polarization waves fall into two main regions corresponding 
to the frequencies of ultra-violet and of infra-red electromagnetic 
radiation, but the wave-lengths in both regions may range from 
values of the order of the lattice distance to infinity. Polariza¬ 
tion waves whose wave-lengths am smaller than the size of the 
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specimen can be divided into longitudinal and transverse 
waves, arftRnTning the crystal to be electrically isotropic. For a 
plane wave moving in the k direction with wave-length 2'jr/Jfc 
and ficequency a>/27r, 

P C 3 C E C3C D = E-f-47rP oc (18.1) 

where D and E are the electric displacement and the electric 
field-strength due to the polarization wave. If a magnetic field 
is also assumed to be connected with the polarization wave, then 
with the use of 18.1 the Maxwell equations lead to the well- 
known formula k = rua/c {n — reflective index). Thus if we 
exclude polarization waves with this wave nrunber, i.e. if we 
exclude the presence of electromagnetic radiation, then the 
magnetic field and the magnetic displacement must vanish. 
Assuming also the absence of free, electric charges and of con¬ 
duction currents, we have p = 0,j = 0, .B = 0, and hence flx>m 
the Maxwell equations conclude that (cf. A 1.1 and A 1.3) for 
polarization waves 

div D = 0, curl E = 0, (18.2) 

or using D = E+^ttP and 18.1, 

div E —47r div P = —47rtkP. (18.3) 

It follows that the three vectors E, P, and D are parallel. Now 
for longitudinal waves k is parallel to P or D, so that 

div D = ikD = 0 

requires that D = 0. For transverse waves, liowever, k is 
perpendicular to D or P, so that 18.3 yields div E = 0, which 
together with the second condition 18.2 leads to E = 0. We 
have thus shown that 

D = 0, i.e. E = —47rP for longitudinal waves; (18.4) 

E = 0, i.e. D = 477 P for transverse waves. (18.5) 

Consider now a spherical specimen whose radius is large com¬ 
pared with, the lattice distance, but small compared with 
c/v€% where v is the frequency of the applied electric field. Prac¬ 
tically we are thus interested in homogeneous—or nearly homo¬ 
geneous—polarization of the sphere. In this case there is no 
difference between longitudinal and transverse waves. Moreover, 
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it will be shown presently that the frequencies of these long 
polarization waves in a sphere are different from those of a 
specimen which is large compared with the wave-length. For 
this case of homogeneous polarization, if I* and Ej, 
instant, the optical and infra-red polarizations respectively 
(of. § 14), the total polarization is given by 

P = Po+P*.. (18.6) 

As shown in § 14, the optical polarization is entirely due to 
the displacement of electrons relative to the atomic nuclei 
(electronic polarization). The infra-red polarization, however, is 
due in part to the displacement of whole ions (considered as 
rigid; atomic polarization) and in part to a displacement of 
electrons which results directly from this change in position of 
the ions. This electronic displacement ‘induced’ by the motion 
of the ions vanishes only in exceptional oases. There are thus 
two displacements of the electrons superimposed upon each 
other. However, since the second one is considered produced by 
the change in position of the ions it is essentially an infra-red 
frequency phenomenon and should be lumped with the ‘atomic’ 
polarization to form P^^.. The two components of the polariza¬ 
tion P^ and Pif. can then be considered as independent of each 
other, i.e. they superpose linearly without perturbing each 
other. If the total polarization were divided into components 
in a different manner—for instance ‘atomic’ and ‘electronic’ 
components—these could not be considered independently since 
the atomic displacement is usually accompanied by a certain 
amount of electronic displacement. 

Assume now that the sphere is brought into a homogeneous 
field Eq due to sources outside the sphere. Then in equilibrium, 
according to A 2.16, the field inside the sphere is 


if is static. 


E = 

Hence, using 1.9, 


static case. 


(18.7) 



3 

4w 



static case. 


(18.8) 


t Vibrationy aro rallod normal vibrations if they superpose linearly. 
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If the external field is not statio, but has a frequency which 
is large compared with the in&a-red resonance frequency but 
aYna.UftT than the optical resonance fiequency, then only optical 
polarization will be excited, and the dielectric constant is c = n^. 
Thus corresponding to 18.8 we find 



3 I p 

47r 2 ® 


(18.9) 


This equation holds for all frequencies smaller than those 
within the optical resonance band, including frequencies smaller 
than the infira-red resonance frequency, although in the latter 
case the relation £ = is no longer valid. 

Kow let a, oLq and (x^,. be the total, the optical, and the infra¬ 
red polarizability of thesphere, respectively. Then by definition 
(§ 14), if F is the volume of the sphere. 


P = ^Eo, P„ = ^Eo, P*. = ^Eo, static case, 

(18.10) 

holds for the equilibrium polarizations in a static field. Hence, 
using 18.9, 18.8, and 18.6, if is the radius of the sphere, i.e. 
F = 4,70^3, 

«ir __ g«t— ^ 1 _■ 3(€,—W^) nKll^ 

< £,+ 2 na-l-2 (e^+2)(w“-f 2)* ' ' ^ 

We shall now consider vibrations of homogeneous polariza¬ 
tion (i.e. polarization waves of infinite wave-length) of a sphere. 
We shall assiune the materud to have one infra-red frequency 
only for this type of polarization and denote by ciij2rr its value 
in the absence of an external field. Assuming the vibrations to be 
harmonic, must thus satisfy the equation = 0 in 

the absence of an external field. Apart from a constant factor, 
the first term in this equation represents the rate of change of 
momentum of the ions, and the second term is the restoring 
force except for its sign. In the presence of an external field E^y 
which may depend on time, the right-hand side of the above equa¬ 
tion must be equal to the force exerted by this field (apart from 
the constant factor), i.e. it must be proportional to Eq where the 



IV, §18 


IONIC CRYSTALS 


163 

This 


proportionality factor is independent of the fiequenoy- - 

factor can, therefore, he obtained by considering the shabio case 

where 0 and where the last equation of 18.10 musb I?® ful¬ 
filled. Thus 

» 2^E.. (18.12) 




2 


47r< 


shall now return to the study of polarization in a 

specimen which is large compared with the wave-length-- To 
calculate the infra-red frequency in this case we shall ooxisider, 
within this specimen, a spherical region whose radius small 

compared with the wave-length hut large compared with, the 
lattice distance* The polarization within such a spheiro is thus 
nearly homogeneous, so that the equations derived a/bov© can 
be applied. In particular, equation 18.12 will hold, if -©6 
represents the field inside the sphere produced by the i>ol»xiza- 
tion of the surroundings. This means that must be jpexiodio 
with the same period as the polarization wave, and it rttay be 
expected to be proportional to Pir (of. below, equatioxx 18.16). 
Now witbJn the sphere such a field gives rise to an. induced 
optical polarization JJ, (proportional to E^) which can Tbe calcu¬ 
lated &om 18.9 as weU as the in£re>-red polarization jg^ iv'cn by 
18.12. Thus, assuming that 18.12 gives the same valnc for 
as was previously considered to exist iu the small isolabed sphere 
with no applied field, there will be superimposed upon bliis an 
optical polarization due to the field E^ produced by the ncmninder 
of the specimen when the small sphere is merely a region in that 
specimen. Hence in a large (compared with the wav'c-length) 
specimen an i]x&a.-red polarization wave is connected with a 
polarization Pir+lu contrast to the case of a small (compared 
with the wave-length) sphere where a polarization wctwe is con¬ 
nected with a polarization Pir only (this was the dejOnibion of 
This means that the polarization due to an infira- 

red polarization wave in a large specimen is composed, of atomic 
and electronic polarization in a different way than ix-i a. small 
sphere.^ We shall show presently that this composibion is 

t v<f. which was a norm^ vibratdon of the ^here, is no long^ex* ct normal 
vibration of a wav© in a large Gfpeoixnen. 
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difforent for longitudinal and. for transverse waves, and that the 
frequencies of these two types of waves are different as well. 

The external held Eq may be expected to be proportional to 
the total polarization P, say 


E. = = 3^(Po+P<,), 


(18.13) 


wher© iili© constant q ■will b© determined laiter. Inserting 18.13 
into 18.9, and solving for yields 

P = g(^^—1) p 

° n®+2— q{n^ —1) 

Introducing this again into 18.13 leads to 

q{n^ — 


(18.14) 


^ E 


= gpjl+-i. _\ = «R 

“ ^n^+2—q{n^—\)) “ * 


n®+ 


w®+2— q{n^ —1)* 

(18.16) 

Hence with the use of 18.11 and 18.16, equation 18.12 becomes 


i.p -i_p _ 


Zq 


CO 




This 


ea+2 w.®+2— q{n^ —1) 


Icr* 


(18.16) 


CO 


== 0, 


(18.17) 


leading to vibrations with a frequency co/ 27 r given by 

to® - e.—71® Zq 

^2 ^ _LO ».2 I o -,/„a Tv* (18.18) 


■8 


To determine q w© notice that the macroscopic electric field 
E inside the sphere can be considered as composed of the contri¬ 
bution Eo of the region outside the sphpre, and of the self-field 
Eg of the sphere. According to A 2.21 the latter is equal to 


—47rP/3, SO that 


E = Eo-:^P, 


(18.19) 


or inserting E from 18.4 and 18.6, 

477 


E, = E+fp = * 


for transverse waves 


(18.20) 


P for longitudinal waves. 
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(18.21) 


Comparison with 18.13 shows that 

1 for transverse waves 
^ — 2 for longitudinal waves. 

Denoting the angular frequencies of longitudinal and transverse 
waves by and a»^ respectively, we find by insertixg 18.21 into 

a>?_»®+2 (18.22) 

■a 


I 


a »? _ 6 , ro »+2 
* n* c„+2 • 


CO 


2 


«a+2 ’ 

The ratio of longitudinal to transverse frequency is thus given 
by (of. references Fll, L4, L5, K2) 


oit \n>) 


(18.23) 


We thus see that the frequency of longitudinal polarization 
waves is greater than that of transverse waves with an inter¬ 
mediate value for the frequency of a sphere whose radius is 
pTWfl.n compared with the wave-length. We also find, inserting 
18.21 into 18.14, that the induced optical polarization has 
opposite Hign for the two types of waves. Finally it should be 
emphasized again that these conclusions hold only for wave¬ 
lengths which are large compared with the lattice distance. 

We shall now proceed to calculate the static dielectric constant 
in terms of the infra-red frequencies. Wq shall again consider a 
sphere homogeneously polarized by a constant external field 
Eq leading to a static polarization P = P^-l-Ptr* The change 

p _J?o in free energy due to the infra-red type of polarization 

can be considered as composed of a self-energy proportional 
to P|r due to the elastic displacement, and the interaction 
energy —P^,.EoF. Thus following a procedure similar to that 
in Appendix A 2.iii we obtain 

F-Po = -PirEo V-\-^G^VFl, (18.24) 

where the constant y in that appendix has now been denoted as 
In equilibrium F must be a minimum if P^, is treated as 
a vector parameter, because and are indei)endent of each 
other in the case of a sphere. Hence, as in the development in 

the appendix, _ 2 ^^ 
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Inserting this expression into the last equation of 18.10 and 
making use of 18.11 we find 


eg —_ 4frr 1 


(18.26) 


Now for elastic displacement a generalized displacement co¬ 
ordinate Q referring to a single unit cell can always be found so 
that the self-energy of the sphero is given by 

(18.27) 


where il^ed is reduced mass of the ions and Nq is the number of 
unit cells per unit volume. This expression must be equal to 

so that 2 ^ ji/ . 2 Qi 

C-a = ( 18 . 28 ) 


Furthermore, the polarization must be proportional to the 
displacement Q and to the number of cells per unit volume. We 
thus can introduce an efiteotive ionic charge e* byf 

P<,.==e*iViO. (18.29) 

Insertii^ the value of (P from. 18.28 into 18.26 a,ji^ making use 
of 18.29 we thiis fiLndf (of. Szigeti, S13) 


eg~7p _ 477 w»-t-2 

. 3 3 




0 


.a* 


(18.30) 


3 3 

If instead of a>g we introduce the transverse angular frequency 
o>| from 18.22 the above equation becomes 



e*^No 


(18.31) 


t It should bo remeixibered here that was defined as polarization due to 
the infiai-red normal vibrations of a ■small ^here. The efiective charge e* is 
therefore a o^uantity which speoifioally refers to a sphere. 

{ An alternative way of deriving eiquation 18.30 would bo to make use of 
the general theory of § 7 according to which (7.21, 7.44) 

u —I 

e,+2 “ SF 3kT * 

where is the average scjuare of the spontaneous fluctuation of the moment 

of a dielectric sphere of volume V in vacuum. is in opr of harmonic 

vibrations the sum of optical and of infra-red terms and is projiortional to ibT> 
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In deriving this formula we proceeded on macroscopic lines, 
except when introducing means of equations 

18.27 and 18.29. Of these two constants the reduced mass can 


always be obtained from the mass of ions and %om a knowledge 
of the structure with the help of Bom*s [-64] lattice theory. 
Calculation of the effective charge e*, however, would require a 
detailed knowledge of short-range interactioxxs and of the charge 
distribution in the lattice. The reasbn for this different behaviour 


of mass and charge is due to the fact that the mass is concentrated 
in the nuclei and contributions by the electrons are negligible, 
whereas the charge cannot be considered as concentrated at 
points except for long-distance interaction. 

As an example consider a crystal of the NaCl type. Let 
and M~ be the masses of the positive and negative ions and let 
r+ and r~ be their displacements. In an oscillation of the sphere 
corresponding to a polarization wave with infinite wave-length 
the polarization is homogeneous at any instant. This means 
that the displacements and r~ have the same value in each 
unit cell, and that r+ is opposite in direction to r“. The restoring 
force must then be proportional to lr+—r~| and act in opposite 
direction on the positive and the negative ions. This suggests 


the relation 


O = r+—r-. 


(18.32) 


because then if we put the restoring force as equal to i O 

we obtain 18.27 as the self-energy of the sphere. This assump¬ 
tion for the restoring force requires 


Af+f+-filf„aa>|0 == 0, ilf-f--Af;ed"20 = o. 

(18.33) 


Dividing by Jfef'*' and respectively and subtracting the two 
equations we obtain, using 18.32, 




(18.34) 


which leads to oscillations of the required angular frequency 


«*>* if 


1 _ 1 I 1 


(18.35) 
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We liave thus caloulated the reduced maes and now require 
effective charge which has to be obtained by calculating tb^ 
polarization of the sphere for a given displacement Q aooorditi^ 
to 1S.29, On the assumption that no non-dipolar interaotiof** 
exists, and that the charges of neighbouring ions do not overla]p^ 
it is easy to show that e* is equal to c, where ie is the charge 
an ion. The assumption. of the absence of non-dipolar forces 
entirely fictitious, however. In fact the restoring force is mainl3^ 
due to short-range repulsion between neighbours* A satisfaotort^ 
calculation of e* has not yet been carried out. 

W^e have thus shown that at the present stage of developmeiih 
of the theory all quantities in equation 18.31 except e* can 
obtained from experiment. Use of this equation, therefore* 
allows a semi-empirioal determination of e*. The following tabl^ 
(Szigeti, 81S) shows that for alkali halides e*/e is smaller than 
unity. This need not be taken as an indication of pronounced 
overlapping of charges, however. It may also indicate that on- 
setting up a homogeneous polarization by a displacement of 
the nuclei of ions, an electronic polarization in the opposite 
direction is induced by the short-range forces. 


Tabus 


(\ = 27rc/a>^) 



u 

n» 

\i X 10* cm. 

e*/e 

laF 

93 

1*02 

32-6 

0*83 

NaP 

6-0 

1*74 

40-6 

0*94 

NaCI 

6*6 

2-25 

61*1 

0*76 

NoBr 

6*0 

2.62 

74*7 

0*85 

Nal 

6-6 

2-91 

86-5 

0-71 

KCl 

4*7 

2-13 

70*7 

0-80 

KBr 

4*8 

2-38 

88-3 

0-76 

KI 

4-9 

2-69 

10-2 

0-69 

Rbd 

6*0 

2-19 

84*8 

0-86 

RbBr 

6*0 

2-33 

114 

0*88 

Rbl 

6*0 

2-63 

129*6 

0-78 

CaCl 

7-2 

2-60 

102 

0-88 

CsBr 

6-5 

2-78 

134 

0-81 

Tica 

32 

6-10 

117 

1-11 

CuGSl 

CoBr 

10 

8 

3- 67 

4- 08 

63 

67 

1-10 

1*0 

MgO 

10 

2-96 

17-3 

2 X 0-88 

CaO 

12 

3-28 

27*4 

2x0-76 

SrO 

13 

3*31 

47 

2 X 0-60 
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For the oxides a factor 2 is separated from the c*/e value to 
indicate that an ideal oxygen ion is doubly charged. We see that 
ie*/e in this case is somewhat smaller but of the same order as 
e*/e for alkali halides. It is also noticeable that the high dielectric 
constant of TlOl does not lead to an excessive value of e*le, but is 
largely due to the high refractive index. In fact, even in the case 
of TiOji, whose dielectric constant is larger than 100 , a value 
|e*/e 0’7 is obtained (referring to the oxygen ions), as has 

been shown by Szigeti The high refractive index, together 
with the high charges of the ions, are mainly rei^onsible for the 
high dielectric constant. 


In the above discussion we have made use of equation 18.31 
and not of 18.30 because tOf and not < 0 ^ can be obtained from 
experiment. Equation 18.31 is very similar to an equation 
derived by Bom \B4\ by an approximate method. This equation 


does not contain the factor 



and it replaces e* by the 


actual ionic charge. 


It should be remembered that, in view of the relation 18,22 


between and equations 18.31 and 18.30 are equivalent. 
This raises the question of the possibility of permanent polariza¬ 
tion of ionic crystals. For on solving 18.30 for it is found that 



4Tr 

3 3 



(18.36) 


On the other hand, equation 18.22, for ->oo, requires 0 
(if we exclude a» 3 ->-ao), which with 18.31 would also lead to 
€3 00. Investigations on these lines should be of importance 

in view of the properties of crystals like barium titanate. They 
have not been developed far enough, however, to be included in 
the present book. 



APPENDIX 

A 1. EI^CTRO-MAGNETIC THEORY 

■ 

(1) Conservatioii of energy 

Maoroaoopic eleotromagnetic theory is based on the Maxwell 
tions» 

curl E = —1 ^ (A l.l), div D = 4«rp (A 1.3), 

curl H = i ^+~J (A 1.2), div B = 0 (A 1,4), 

where H and B are the xnagnetio field strength and induction respec¬ 
tively, p is the density of the true charge, j the density of the conduction 
current; for E and D of. § 1. To allow a unique calculation of the field 
vectors E, D, H, B &om p and j, the Maxwell equations must be supple¬ 
mented by two fiirther relations of the type 

E = E(D) and H = H(B). (A 1.5) 

These relations are not included in the firndsunentcd equations for the 
eleotromagnetio field, but are ohoraoteiriBtic of the type of material 
employed. For substances of interest in. this book 

H = B (A 1.6) 

can be assumed to hold to a good approximation. 

To introduce the energy law, multiply equation A 1.2 by E, and sub¬ 
tract fipom it equation A 1.1 multiplied by H. Use 

H curl E—E curl H = div[ExH3. 

Then with equation A 1.6, after multiplication by o/4ir, 

Books on eleotromagnetio theory discuss in detsiil the fact that the 
last term of this equation represents, per unit volume, the rate of con¬ 
version of eleotromagnetio energy into other types of energy (heat, 
kinetic energy of particles, etc.) in so far as this is connected with con¬ 
duction currents. o[E xH3/4i7 is the Poynting vector representing the 
rate of flow of energy; the second term, therefore, gives per unit volume 
the rate of efflux of eleotromagnetic energy. Then from oonseivation of 
Gneigy it follows that the first term must represent the rate of change 
of energy content per unit volume, provided that there is no flow of energy 
of another type (e.g. heat currents). 

Now -Lh—= -~ 

47r at at Bn* 

so that H*IBn is readily recognused as the density of the TnA.gnftt.f» energy. 
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It is, then, often suggested that 

*(D)dD 

should represent the density of the electric energy. The integral is con¬ 
sidered to extend from E = 0 to the actual value. In fact, if E == D/e,, 
on© finds, for a field-strength 

(A 1.8) 

0 

i f e, can be considered as constant. Such a procedure is, however, 
only if E is a univalent function of D, for otherwise the integral no 
uxdque meaning. Thus only in this case can the electric energy be defined 
without further investigation. 

In general, the conclusion to be drawn is that 

^EdD (A 1.9) 

represents the change of energy density (not necessarily electric enCTgy) 
connected with a variation of £> by dl>, if there is no fiow of ^ergy of any 
other type. It has been shown in § 3 how the total energy can be calcu¬ 
lated from this expression. 


(il) Conduction current and energy loss in periodic fields 
As in § 2, assume a homogeneous field with JB = EoCostui, and 

, wf . 3JS 

D = €iBQOOBo}t+€^B^em.€ut = 

Then using and assuming the conduction current y to 

vanish, the right-hand side of equation A 1.2 becomes 

1 dD Cl dE , cyw B? 1.10) 


e dt c 0* 


If, on the other baud, j ^ 0 , and if Ohm’s law holds, i.e. if (a is the con¬ 
ductivity) y ^ 

fLTiri if furthermore Cg = 0, then the right-hsmd side of equation A 1.3 

!i ^+tE^E. (A 1.11) 

c at ^ c 

Comparing equations A 1.10 and A 1.11 it follows ttot “ 
fields the introduction of two dielectric constants €i(ai) and «*(a») (tf. ^ 
equation 2.8) is equivalent to using a single e^Cc,). put^ e, ^ 

introducing instead a frequency dependent conductivity <t(c«i). 
representations are coimected by 

, . a)eg(ci) fA 1.13) 


4980.11 
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JoTjle’s law for the rate L of loss of energy i>er unit volume then becomes 

i = ctE* == io--Eo» 1*13) 

where the bar has been used to represent tbe average oyer one period. 
Tnaart.mg a from A 1.12, this expression is found to be equivalent to oqua- 

tion 3.15. 


(Ill) Relation between €i(co) and €^(€ 0 ) 

To derive equations 2.16 and 2.17 apply the theorems of Fourier 
tronsformatioiis to equations 2.14 and 2.15. 

Then from 2.14, 




—e^}oos^ dft. 


(A 1.14:) 


«0 


and from 2.15, 


a(») = J / C2(/x)siiifU]? dfjb. 


<A 1.15) 


Introducing A 1.16 into 2.14 yields 

00 yr OO 

ooscurc 


€l(co)—€oo = “ J r dfjbj 

0^0 ' 

00 / JB % 

= ? lim r J coscorrsin/xa; dxj 


Q 

eo 


=:^lim r d^. 

^rjfr-Mo J ^ u + oj fx —a> / 

0 

(A 1-16) 

Now the integrals containing the oos-terms vanish if jR—oo. The re¬ 
maining terms lead immediately to equation 2.16. 

Finally, introducing A 1.14 into 2.16 gives 


00 / 00 . 

6 t(ai) = ~ J clajfsincoa? J {€i(/x)—€^}oos/Lta: (A 1 . 17 ) 

A ' rt / 


00 . oo 

smcoo? 

0 

which 3rields equation 2.17 by a aimileu: method. 


(iv) Relations between dielectric constants and optical constants 

From the Maxwell equations A 1.1—A 1.4 the wave equation can. be 
derived. For this purpose apply the operator curl to A 1.1 and 1/c &/3t 
to A 1.2. Then, using A 1.6, H can be eliminated. Assuming periodio 
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solutions, we have D = e(ei>)E, and hence in the absence of free charges 
(i.e. p = 0) and of conduction currents {j => 6) we fbad 


V*E—- = 0 

c* 


(A 1.18) 


if use is made of the relation (since from A 1.3 and p = 0, divE = 0) 

curl curl E = graddivE—V*E == —V*E. 

Now assume a solution representing a wave penetrating the dielectric 
m the a^-direotion, ^ 1,19) 

wheire A is a constant vector perpendicular to the :r-direction. Sere by 
the usual definition of the optical constants n is the re&active index and 
K is the absorption coefficient. They can be expressed in terms of the 
complex dielectric constant e (cf. 2.8), for after inserting A 1.10 into 

A 1*18 we find (n+^ic)* = € = ^i+iea- 


Hence 

= n*—/c* 

(A 1.20) 


= 2nic. 

(A 1.21) 

Alternatively, using A l.l2. 

27r 

W#C sss "■"**" (T* 

CJ 

(A 1.22) 


A 2. DIPOLE MOMENTS AND OTHER ELECTROSTATIC 

PROBLEMS 

(1) The basic problem 

"W^e consider an infinite homogeneous dielectric with static di^ectnct 
constant which contains a spherical region of radius a and dielectrio 
constant W^e wish to calculate the electric field due to any of the 
following sources: 

(a) Outside sources leading to a constant field E^^ (s^^y hi the z-direction) 

at a large distance from the sphere. 

(&) A point dipole (jl (say in the 2 -direction) at the centre of the sphere, 
(c) An extended dipole M (say in the 2 -direotion) inside the ^here 
arising fipom a homogeneously polarized sphere of radius a, 

M = ^a»Pe. (A 2.1) 

Pp being the polarization. 

Let O be the electrostatic potential so that 

E = —grad ® 

is the field-strength. ® satisfies the Laplace equation 

V* 9 = 0 

t Iix oontrajst to the denotation in Uie sections, and c, now represent 

static dielectrio constants. 

4080.11 3f2 
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subject to the following conditions: if r is the distance from the centre of 
the sphere and 8 the angle between r and the z^axia, then in 


case (a) 

B == Le. ® = —jEro^roos^ if r ^ a. 

(A 2.4} 

case (b) 

«a 

(A 2.6) 

case (c) 

D = €2E+47rP^ if r < u, 

(A 2.6) 


D == ciE if r > a. 

(A 2.7) 


where D is the electric displacement. In aU three cases the normal com¬ 
ponent Df. of D and the tangenticd component E 0 of E must be continuous 
at r = a. 

The ccdculations are simplified by treating all three cases together 
because they all lead to the same angular dependence of In general 
theory <D is developed into a series of spherical harmonics. Owing to our 
three conditions, and to the boTxndary conditions, only terms proportional 
to COS0 appear. For this angular dependence the general solution of 
A 2.3 is given by 

® 

containing two arbitrary constants A, B. They will have different values 
outside (Aj, B^) and inside (Aj, B 2 ) the sphere and have to be determined 
from the boundary conditions. Using A 2.4 it follows that = Bqo, and 
hence . . * 

® 

Also using A 2.5, we have —Aj. — i-o* 

a> = 

According to the boimdary conditions 

must be continuous at r = a, so that 

*^1 I tji _ 


T > a. 


(A 2.8) 


a. 


(A 2.9) 


a' 


3 --^+^ 2 . 


(A 2.10) 


Similarly, using A 2.6 and A 2.7, continuity of Doosd at t = ct 

yi^ds 

+^i-®oo *= 


a* 

From these two equations one fi-ndg 

® 2 €i +€2 a* 


(A 2.11) 


i3 


'a 


3«1 r, , 2 ei—«, n 


-»«+ 


4irJ^ 


2«i+ea * cj 2€i+e, 0 “ 2€i+«,' 


(A 2.12) 


(A 2.13} 
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Inserting A 2.12 and A 2.13 into A 2.8 and A 2.9 gives the solution of the 
general problem. 

For the discussion we shall now separate the three oases. 

Case (a): No internal sources, /a = 0, fj, = 0. In this case the field 
inside the sphere is equal to and will be denoted as cavity Reid 0\ 
Thus from A 2.13 with A 2.2 


is the field inside the sphere. If in particular = l» i.e, for an empty 
sphere ^ 

GVa = 1) - G « (A 2.15) 

If, on the other hand, the sphere is in vacuum, then = 1, and the field 
inside becomes 

(A 2.16) 


The field outside the sphere according to A 2.8 is composed of the field 
Eqo at infinity and of a dipolar field with potential (using A 2.12) 


— 

2€i+ea 




cos^ 


(A 2.17) 


OcLse (6): Point dipole, no external field, E^q = 0, jFJ, = 0, In this case, 
according to A 2.9 and A 2.2, B 2 represents the deviation of the field 
inside the sphere from a purely dipolar field, i.e, B 2 ib tlie reaction-field 
R' acting on fjL- Hence with A 2.13 


R' 


(A 2.18) 


^ /i. 

2 € i " f “€‘2 

if in particular ^2 = 1, then the reaction field will be denoted by R, and 

€1—1 2 


R = where 17-2e, + Ia»* 

Outside the sphere wo have a dipolar field with potential 

^ _ _ Ai COB0 _ 3 fjLCoad 


(A 2.19) 


Case (c): Homogeneously polarized sphoro, no external field, Egc — 0, 
/X = 0. Consider first the sphere in vacuum, i.o. Ci = €2 Then B 2 

represents the self-field inside the sphere which using A 2.13 and A 2.1 is 
thus given by 

A<Tr rVI 

{A 2.21) 


V - *Et> _ ^ 


3 


a* 


Now surround this sphere by a medium with dielectric constant €|. Then 
the increase of the field inside the sphere, —is the reaction-field R. 
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Thus using A 2.13 with = 1, 

2€i+ 1'^3 2ei-|-l 3 2«i + l o® 

(A 2.22) 


in foimal agreement 'with A 2.19. 

Thus, if the sphere contains a point dipole p. surrounded by a homo* 
geneously polarized sphere 'with moment M, then its reaction field is given 

^ R = |,|^(M-|-p) = er(M-f-|*). (A 2.23) 

where g is defined in A 2.19 euid M+ji is the total dipole moment of the 
sphere. 

For regions outside of the sphere also oaise (&) turns into case (o) if (A 
is replaced by M. This holds even if e, ^ 1. Thus as far as the field in the 
region external to the sphere is oonoemed, cases ( 6 ) and (o) lead to the 
same results. 


(11) Dipole mom^ts 

The potential of a rigid (non-polarizable) dipole with moment in an 
infinit e medium of dielectrio constant is given by 

fl> = dfs (A 2.24) 

ei r* 

Now let us consider a model of a molecule consisting of a rigid dipole |a 
at the centre of a sphere with dielectric constant € 2 - (The same model will 
of course apply to an extended dipole with the same moment correspond¬ 
ing to case (c).) If this molecule is embedded in a medium of dielectric 
con s ta n t then according to A 2.20 (emd to the remarks at the end of the 
discussion of case (c)) the potential outside of the sphere can be repre¬ 
sented by equation A 2.24 if 

Therefore is defined as the external moment of the molecule in the 
medium 

It should be realized that the moment of the molecule in vacuum, (Af;, 
is difierent from {x. For is defined as the external moment in a medium 
for which € 3 ^ = 1 . Thus the vacuum moment is given by 

Pt, 

It follows that the external moment p^ in a medium can be expressed 
in terms of p^. Using A 2.26 and A 2.26, 

^ In contrast to the external moment p^—^which is the moment a point 
dipole must have to produce the same field as the molecule—the internal 
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moment {JU( is the aotual moment of the moleonle if embedded in the 
medium ex* In vacuum (ex = 1) the external and internal moment are, 
of course, equal. In another mechum, the internal moment differs firom 
the vacuum moment by the polarization of the molecule through the 
reaction field. Hence equation 6.18 for jiq is obtained and this, in turn, 
with the help of A 2.19, leeuls to equation 6.20. Equation 6.20 can also be 
derived by using directly the properties of the external moment. As we 
have seen above, the field produced outside the sphere does not depend 
on whether the moment is produced by a point dipole or by uniform 
polarization of the sphere if this leads to the same moment. Therefore 
the molecule (as represented by a sphere with dielectric constant and 
a> dipole fjL at its centre) must have the same moment as a sphere with 
dielectric constant ^x which is in a medium with the same dielectric con¬ 
stant €x and has a moment at its centre, because by definition of 
both produce the same field outside the sphere. Thus 

J P diTf 2.28) 

spliflre 

showing that is composed of the moment of the rigid dipole and the 
moment contained in a sphere around it. Now if points in the z-direc- 
tion, the integral must be a vector with the same direction for reasons of 
symmetry. .Hcuce using 

47ri^==(€x—1)^^ = —(cx (A 2-29) 


and 


OOS0 


_ i _ _ 

dz\r}^ 


we find with the help of A 2.24 that 




(A 2.30) 


smce 


and 


/I 

dx^\r 


j - -4ir. 




Insorting A 2.30 into A 2.28 yields 


€l-l\ 


(A 2.31) 


Tjir'hicll. IS idonticfitl 6*20« A • J 

The above calculation also sho’jro that the di]^le “ 

in a sphere surrounding a dipole is independent o * ® the 

Hence the moment contained between two sphenca 
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dielectric vanishes* This holds even if the two spheres are not con- 
centrio.f 

Thus, if we consider a large sphere inside an infinite dielectric of dielec¬ 
tric constant Ci, and this sphere contains a molecule with internal monaent 
p-i, then the moment of the lai^e sphere is also 

This is, however, no longer the case if the large inhere is not embedded 
in its own medium because then the field outside the molecule is different 
from the field of a point dipole* By analogy with A 2.26 one should expect 
the moment in this case to be 

I*. = 

as is actually confirmed by calculation if the radius of the sphere is large 
compared with that of the molecule. 

We have thus introduced five different dipole moments, p, /Ai, /Xg, 
ft,. To prevent confusion we give a brief summary. 

/t has a meaning only in terms of the special model used. 

ft, = moment of the molecule in vacuum. 

ft| = moment of the molecule in a medium 

= moment of a sphere (containing the molecule) within an ixiflnito 
medium 

p, = mo me n t of a rigid point dipole producing in a medium the same 
dipolar field as the molecule. 

ft, = moment of a dielectric sphere in vacuum containing the molooule. 
To summarize the formulae: 


=* 


ea+2 2 €i+1 

_gi+2 2€^+gj 


Pi = 




«s+2 2ei + l‘**~6,+ 2 


P. 


(A 2.33) 


^ direction, say -with radius vector s, is proportional to 

Tali^ace between the two spheres with 

rj. Thismtegral can be transformed into the difference of two 


rj J cos ^ dn —rj f cos ip dQ. 

m T 


the mirfaGAA -noin'i-tsi js^ • is the angle between s and the normal to 

Into qilMric.1 *' «•»'>• developed 
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(iii) Self-Energy 

The free energy of a dielectrio polarized by a field can be considered 
as composed of three terms: (1) the free energy of the field and of the 
dielectrio before the latter is brought into the field; (2) the energy of 
interaction between the x>olarized dielectrio and the field; and (3) the free 
energy required to polarize the dielectric. The latter will be called the 
self-energy Fg. If the dielectric is polarized homogeneously, then it will 
be assumed that Fg = yJlf where M is the moment of the dielectric. 
The energy of interaction is equal to — (MEq) if Eq is the field (supposed 
homogeneous) before the dielectric is inserted. Thus 

JP-JPo == -{MEo)+yJWr*- (A 2.34) 


Tri equilibrium F must be a minimiim where M = (ikfa;, M^) is treated 
as a parameter. Hence from 


dF _ dF _ BF 


(A 2.35) 


we find 


2y]M = Eq 


or 



(A 2.36) 


Hence our assumption Fg = yM^ implies that M is parallel to Eq. The 
value of y depends on the shape of the dielectric. Thus a sphere of radius 
a brought into a homogeneous field Eq has a moment 

M = a®Eo(€,^l)/(€,+2). 


Hence 


y = o 


1 + 2 1 


2 €g — 1 a 


3» 


leg+2M^ 
*"^ 2 €a— 1 a3 


(A. 2,37) 


A slab of volume V with surface perpendicular to Eq, on the other hand, 
has a moment M == T^EQ(cg— l)/47r€a. In this case 

277 €g — 27 tM^ 




e,-l V 


(A 2.38) 


In both cases, of course, Fg is proportional to the volume if the polariza¬ 
tion M/V is constant. 


A 3. THE CLAUSIUS-MOSSOTTI FORMULA 

In § 6 an equation for the static dielectric constant €g was derived 
(5.13) which holds exactly for the model used in that section. This is 
usually called the Clausius—Mossotti formula. Later in § 6 it was found 
that a similar formula holds approximately for liquids of non-polar 
molecules, but that the approximation required negligible short-range 
interaction. In § 8 it was pointed out that by neglecting short-iange 
interaction the general theory of § 7 will lead to either the Clausius— 
Mossotti or to the Onsager formula, depending on whether one considers 
non-polar or polar spherical molecules. The proof of this will be given at 

the end of this section. 

Discussions as to whether or not the Clausius-Mossotti formula holds 
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exactly have formed part of the literature on dielectrics over a period of 
many years. Controversial conclusions have been drawn, mainly due to 
a misunderstanding of the significance of this formula. In fact one should 
distinguish between a znacroscopic cuid a molecular formula. The same 
mathematical symbols are normally used in the two oases and, as a 
result, they are often ooxzfused with each other. The macroscopic formula 
is valid exactly, but the molecular formula holds only subject to the 
conditions mentioned in § 15. 

To derive the macroscopic formula, consider a sphere of a continuous 
isotropic dielectric to be brought into a constant electric field f. Accord¬ 
ing to equation A 2.16, the field inside the sphere is then given by 



(A 8.1) 


if is the dielectric constant. Let Mjg be the dipole moment induced in 
the sphere. Then the polarizability of the sphere will be defined by 
{the suffix m indicates ^mctcroscopic’) 



(A 3.2) 


On the other hand, in view of 1.9, 

= (A 3.3) 

»7r o 


if F is the volume and is the radius of the sphere. Equating A 3.3 and 
A 3.2 and introducing / from A 3.1 leads to the Clausius-Mossotti formula 


^ 8+2 


(A 3.4) 


An alternative way of deriving this formula is to consider a spherical 
region within a homogeneous dielectric. If a constant field E is produced 
in the dielectric, the part of the field inside the spherical region which is 
due to sources outside this region (the izmer field) is given by expression 
6.9 and is therefore identical with our field f (A 3.1). The remaming 
development is then identical with the one given above. 

Equation A 3.4 is alwa 3 ^ correct if the spherical region is sufficiently 
large, so that the material contained in it znay be oozisidered from a 
znacroscopic point of view. Prom a molecular point of view equation 
A 3.4 has, however, no significance. To give it such a significance the 
polarizability oe must be expr^sed in terzns of other quantities whose 
values will not be found experimentally by measurement of the static 
dielectric constant but by experiments of a different nature. 

Equation A 3.4 finds its main application in the case of dielectrics in 
which polarization is connected with elastic displacement of charges 
(case (i) of § 4). The most general displacement of the dielectric material 
in a spheire may then be developed in terms of its normal vibrations 
(cf* van Vleck, F3), and it is then found that the polarizability a is a con¬ 
stant independent of temperature if the density of the substance is kept 
constant. This procedure is of importance because it shows in a very 
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^neral way that the assmnption of Mastic displacement leads to a 
iemperature-independeixt dieleotrio constant.. 

We shall now disoiiBS the ^molecular’ dausina—Mossotti formula* In 
}his case it is assumed that the maorosoopic polarizability can be 
expressed in terms of a molecular polarizability <x —which is a constant 
independent of temperature or density—by 

ot^ 5f=5 Noi (A 3.6) 


if ^ is the number of molecules in the macroscopic sphere of radius 
Prirthermore, if 


4ir 

T 


a® « — 
® N 


47T 

3 JV 


(A 3.6) 


is the volume occupied per moleotile, insertion of equation A 3.6 and 
A 3.6 into A 3.4 yields 




a 


. 8 * 


(A 3.7) 


€,+2 a * 

This molecular foimula has the same mathematical structure as formula 
A 3.4, but its significance is different. In equation A 3.7 the polariza¬ 
bility CK is a property of a sin^e molecule independent of inaorosoopio 
parameteirs. The quantity 3/4irex® is equal to the nimiber Nq of molecules 
per imit volume which is determined hy the molecular weight TV and the 
density d, 

^ No ^ , (A 3.8) 


8 


47ra® 


where A is Avogadro ^s number. Hence cc/cfi is proportional to the density 
of the dieleotrio, which may be varied by altering the external pressure. 
No such concliision can be drawn firom equation A 3.4, where the macro¬ 
scopic polarizability (Xn^ may depend on the density in an unspecified way • 
The crucial step leading from equation A 3.4 to A 3.7 is contained in 
the hypothesis A 3.6. One should expect A 3.6 to be correct only in the 
absence of short-range interaction between molecules, because such an 
interaction would influence the recwstion of a given molecule to the 
field. Actually in the two cas^ in which we have derived the molecular 
Clausius—Mossotti formula (6.13 and 6.34) sho3rfc-range interaction was 


absent. 

Finally let us derive A 3.7 and the Onsager formula hy applying the 
general theory of § 7 to a model in which short-range forces ax© entirely 
absent. Also it will be assumed that the unit cell is spherical and that the 
average polarization outside the unit cell is equivalent to that of an 
isotropic continuum polarized by the moment m of the unit cell. In 
this case the average moment m* of a sphere (embedded in a larger speci¬ 
men) oontaming the unit ceU is equal to the moment m of the unit cell, 
as was shown in § 7. Hence the sphenoal region introduced in § 7 may he 
taken aa equivalent to the unit cell^ and equations 7.11 and 7,33 become 
identical if 7.34 with N = 1 is taken into account. Thus, using A 3.8, 


3€, ^ 


(A 3.9) 
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where m* is the aveorage square of the gpontaneous dipole moment of the 
spherical cell embedded in its own medium. Assuming first that the unit 
cell contains a rigid dipole with moment m = and w* = ft®, the 
Onsager formula follows ficom A 3.9. Consider, on the other hand, a 
spherical polarizable molecule without a permanent moment. Such a 
molecule can be represented by an elastically bound charge e. Then if r 
is the displacement, the internal energy must be quadratic in r, i.e. 


Ki=|r». (A 3.10) 

where c is a molecular constant. To express this constant in terms of the 
polarizability, assume the molecule to be brought into an external field f. 
Then its potential energy is given by 


|r«-cfr. (A 3.11) 

A 

The equilibrium value r of r is obtained by making this expression a 
minimum. Hence ^ 

e¥ = -t (A 3.12) 

c 


is the average moment of the molecule. Since this must be equal to otf, 
we obtain 



e* 

(A 3.13) 

and hence with A 3.10 

08 

= IS 

(A 3.14) 


To obtain nfl we require the energy U = because, as in 7.12, 

using m = er. 


CO 


TO* = 






t 

// 


dr. 


(A 3.15) 


Now the external energy is given by —where R is the reaction 
field, as can be shown by the argument which led to 7.18. Therefore 
making use of 5.10 and the fact that 47 ra ^/3 = V, 




TO 


2 


•8 


a 




iS 


a® 2€5+1 a®2€,+ l 

Introducing U = into A 3.15 and integrating yields 

kT 2a®(2€j+l)cx 


(A 3.16) 


3 

m* == — 


2 aa(2€,+ l)-2(e,-l)a- 


(A 3.17) 


Inserting this expression into A 3.9 leads at once to equation A 3.7. 

In the above derivation the importance of the external energy 
inoluding the interaction with the reaction field should be stressed 
It is iateresting to compare this with the case of a non-polarizable 
molecule with dipole moment {jl (Onsager case), where the interaction 
with the reaction field is irrelevant. In this latter case the only variable 
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is the direction of (jl, but its amount, fi, is a constant. Thus since m ^ 
in this case, m* = = constant, so that according to A 3.16 is a 

constant as well. Hence in A 3.15 the terms containing JJ^ can be taken 
outside the integrals and, therefore, cancel. 


A 4. SHAPE OF ABSORPTION LINES (F9) 

We shall consider a set of linear harmonic oscillators of mass Mt 
charge e, proper frequency aio/27r. Their number N per unit volume is 
assumed to be so small that interaction between them can be neglected. 
Then in an external field E, which may depend on time, 

x= —J?oos0 (A 4.1) 


if a; is the displacement of the charge e, 0 the angle of E with the x-direc- 
tion, and a dot represents differentiation with respect to time U T^e 
oscillators are supposed to make frequent collisions with a medium which 
is in thermal equilibrium so that they also tend to reach equilibrium. To 
describe their behaviour we introduce a distribution function /(»,») bo 
that /(a?, x) dxdx represents per unit volume the number of oscillators 
with a displacement between x and x+dx, and a velocity between cb and 

x+dx. Hence 

J J dxdx = JV. (-A- 4.2) 

—00 


With the help of the distribution function the polarization P in the field 
direction is obtained as 

(A 4.3) 


CD 

P = 6OOS0 J J a;/(z,z) dxda. 

— 00 

Let us first consider the static ease for whi^ the field is independent of 
time, = J?o* If is the energy of an oscillator, i.e. , , , 

U = Di- eEx cos 9, (-A- *■*) 

Uo = iMwi **®) 

then, in equUibrium, according to the Boltzmann theorem, 

(A 4.6) 

where O is independent of x and x ^d^n determined with 
of equation A 4.2. Now using A 4.4 and o. 

( eJrzcoafl 1/e^cos^* I (A 4.7) 

g-mxAI^T = -fcji- kT / ^ / 

and hence we can write, considering A 4.5, 

* -nr 0 PA 


eEcosB ^/o 1 
f—fo— dx^”" 


(A 4.8) 


where the dots represent SterSed in linear 

In the calculation of the polarization P we snau 



174 


APPENDIX 


terms in SJ only. Therefore we need not consider these higher terms; to 
be exact we should also prove that the series for P, in terms of JE^ con¬ 
verges, but we shall omit this. In A 4.8 

f^=: Oe-^olkT (A 4.9) 

represents the eqxiilibrium distribution in the absence of a field. The fac¬ 
tor <7 is the same as in A 4.6 if in ceJculating the integral A 4.2 for the 
latter case only terms which are lineM in, or independent of, jE?© con¬ 
sidered. 

Now inserting A 4.8 into A 4.3 we find for the polarization in a con- 
■ Btant field 


— oo 

Here the first term vanishes because is an even function in a?. Inte¬ 
grating the second integral by parts we find, using A 4.2, 




e^Eoos^dN 

Moj^ 


(A 4.11) 


Averaging over all directions of the field leads, of course, to a replacement 
of oos*0 by 1/3. Hence the contribution Ae to the static dielectric con¬ 
stant is given by 


E SAftog* 


(A 4.12) 


Consideir now the case of a periodic field E described by 

E ^ E^e-^^K (A 4.13) 

Again we shall introduce the distribution function^ satisfying A 4.2, but 
•now dep endin g on time We shall calculate / by considering its rat© of 
change Qfjdt for a given veJue of x and x. This is composed of two terms, 
one due to collisions of the oscillators with the surrounding medium and. 
tibfi. other due to their motion. Denoting by (S//^)cou the rate of change of 
/• due to collisiozis the simplest assumption we can make is to put 




fwvi)r 


(A 4.14) 


in which the relaxation time t is assumed to be independent of x €U 3 id d>» 
This inoplies that equQibrium (/«iu) Is approached exponentially. 

Eor time-dependent fields we shall assume that A 4.14 still holds if 
/miu Is the equfiibrium distribution corresponding to the field at the time 
in. question. Thus making use of A 4.8 and A 4.13, 


Tetia — fi 


co^ % 
dx* 


Ma>% 


(A 4.16) 


For the distribution function f in this case we shall assume 

/(®, X, t) = /j{ar, x)+gix^ x)E^ cos 0 


(A 4.16) 
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where g{x ,») is a funotion of x and x, indepeadent of t. To satisfy A 4tJi 
this requires 


Jj g(Xtds)dxdA 0. 

—00 

Inserting A 4.16 and A 4.16 into A 4.14 now leads to 


(A 4.17) 


“ “^(‘'+3^ g)ai.^<oos8. 


(A 4.18) 


There is a second contribution to dJIdt due to the motion of the oscilla¬ 
tors. This contribution will be denoted by so that 

at 




(A 4.18) 


To find on expression for notice that all the oscillators contained, 

at a time in an interval the number of such oscillators being 

/(a?, ac, Aa?Aac, were at the time ^—Si in an interval of the same sizetAasAi 
with coordinates x—xBt provided the number has not been 

changed by collisions; all^(a!, jc, f—8i) Ag?Agft oscillators contained in thi s 
interval at the time t —Sf have meanwhile pctssed into another interval* 
Thus 


L{f(x—xBt, x—x 8<, *— 8t)—J{x ,», S*)} 

ot 


is lihe rsite of change of f due to motion. Hence with St 

dl 

BtJ 


0 , 


^Ix-^x 

dx Bx 


(A 4.20) 


or making use of A 4.1, 

\dtJ^ 


•XA-tOn ^ . 
dx ^ ^dx 


M 


dx' 


(A 4.21} 


Now from A 4.9 and A 4.6 it follows that 

flC+CuS 


I_ A^x = 0. 


(A 4.22} 


dx - '’Bx 

Thus inserting from A 4.16 into A 4.21 and neglecting terms proportkmal 
to E^, we obtain 




e""*"* COB 0. 


(A4J23) 


^ = —itx)gEo er-*"‘ cos 8. 
dt 


(A 4.24) 

Also from A 4.1t5, 

Inerting A 4.24. A 4.2S. »d A 4.18 into A 4.19 „ find ttot , 

satisfies the differential equation 

/ 1\ Bg . , .Bg^ (A 4.26) 

Md± Ma^rBx 

This can be proved with the help of A 4.1. 
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To solve it we put 


g ^ a 




(A 4.26) 


dx ^ dx 

where a and b are constant. 

This does not represent the most general solution of A 4.26> but it has 
been shown by Huby [IfS] that it is the only solution, satisfying all 
requirements. 

The polarization P according to A 4.3, A 4.16, and A 4.26 is thus given 
by 

^ 00 

P = r r ==" — e008*6fit- 

-o5 ^ (A 4.27) 

To find the constant a insert A 4.26 into A 4.25. M ak i n g use of 

(A 4.28) 


and 


we obtain 


Hence 


/ . 6 , « a Wo 

I dx'^ ^ dx/dds ~~ dx * 


(A 4.29) 
(A 4.30) 
(A 4.31) 


a = 


e a»oT*+l —icoT 


aioT*+(l—iair)* 

6 1/ 1— ioi^T , 14 ~^<t>oT \ 

2V1— i(ta-i-tOo)T 1 — i(at — tanW/' 


M.tiA 


(A 4.32) 


Together with A 4.27 and A 4.12 this leads to eqTiation. 13.9. 
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Magnosiuia oxide, 158« 
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MOodels, for crystaUine solids, 54 

— for Debye equations, 78 
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63 
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Molecules, 105 £E. paaavnt^ 

Neon, 111. 
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Onsager’s formula, 33 ft., 49 ft., 53, 
130, 171. 

Optic^ constants, 14, 162. 

—jxjlaiization, 105 ft., 14:9. 
Order-disorder trcmsitions in dipolar 
crystals, 53 ft., 132, 146. 

Oaqrgen, 108, 111. 

Paraffin, 108. 113 ff., 125 ft. 
Penta-methyl ohlor-benzene, 138. 
Pentane, 114. 

Phenolic resin, 135. 

Polar molecules, 26 ft., 105 ft. passim. 
Polarizability. 28, 105 ft. 

Polarization, 2, 64. 


Polarization waves, 149 ft. 

Power loss of dieleotrios in periodic 
field, 13, 73 ft., 92 ff„ 98 ff- 

Kate of unimoleoular reactions, 81. 
Keaction field, 25, 31, 34, 41, 165. 
Kefractive index, 28, 163. 

Kelaxation time, 73 ft., 8S, 01 £t«, 
121 ft., 128 ft. 

Kesonanoe absorption, 98 ft., 173 ft. 
Bubber-aulphur compounds, 135. 
Kubidium halide, 158. 

Self-energy, 144, 145, 155, 169. 
Sodium chloride, 109, 157. 

Strontium oxide, 158. 

Supeiposition principle, 6. 

Tallium chloride, 158. 
Temperature-dgpendenoe of cLieleotrio 
oonstant, see Dieleotrio constajnLt. 
Tertiary butyl chloride, 138. 

Titanium oxide, 109. 

Viscosity, 84 ft., 89, 123. 

Water, 107, 137 fit. 

Xenon, 111. 
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